{"id":"47bf1c7e-b1bd-458e-860d-40fefec8f95f","arxiv_id":"1908.04761","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"DFT calculations predict that compressive strain makes the FeRh/MgO easy axis in-plane while tensile strain makes it out-of-plane, with the switch at zero strain.","lead":"Using density functional theory calculations, the authors report that strain in an ultra-thin FeRh film on MgO changes the magnetocrystalline anisotropy and can flip the easy axis of the magnetization. The work points to strain as a possible control knob for antiferromagnetic spintronic memory applications.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The strain parameter co-strains the MgO substrate, so the predicted MCA switching may not be physically achievable for FeRh/MgO; the central claim is unsupported without separating film strain from substrate strain.","rationale":"The paper presents a plausible but incomplete computational study. The central quantitative evidence (Table I, Eq. (1)) is missing, and the abstract/results directional inconsistency confuses the reader. My principal concern is more specific: the way strain is applied in the slab calculation physically co-strains the MgO substrate. Because the in-plane lattice constant is the only lattice parameter in the supercell, varying s relative to bulk FeRh necessarily strains the MgO by about +0.6% to +1.1% tensile, which is not achievable for a thick rigid substrate. The MCA in Fe/MgO systems is known to be sensitive to the Fe-O distance, and here that distance changes because the O atoms belong to the strained MgO. The paper does not separate the FeRh strain effect from the substrate strain effect, so the predicted easy-axis switching could be an artifact of the model rather than a property of FeRh/MgO. A free-standing FeRh slab calculation under the same biaxial strains would settle this. The reader's weakest assumption about interface termination and thickness is related but not identical to this substrate-strain issue; I partially agree with the reader. Given these unresolved issues, the conditional verdict remains appropriate, but the conditions should explicitly include the strain-separation test and the restoration of the missing data.","tokens_in":4664,"tokens_out":18067,"duration_ms":179959,"concrete_test":"Compute the MCA of a free-standing 5-ML Fe-terminated FeRh(001) slab in the G-AFM phase under the same biaxial strains (s=-0.50%, 0, +0.50%) using the same PBE functional, cutoff, k-mesh, and relaxation protocol as in Section 2, but with no MgO substrate. If the MCA sign still changes from negative to positive over this range, the switching is intrinsic to FeRh under lateral strain; if it does not, the reported switching is an artifact of co-straining the MgO substrate and the FeRh/MgO claim fails. Additionally, recompute the FeRh/MgO slab with the in-plane lattice fixed to bulk MgO (4.21 Å) and report the MCA at that natural lattice-matched interface.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the slab model (Section 2), the in-plane lattice constant is a single parameter for the entire FeRh/MgO supercell. The strain s is defined relative to the bulk FeRh lattice constant (2.995 Å, with FeRh[110] aligned to MgO[100]). At s=0, the in-plane FeRh lattice constant is 2.995√2 = 4.235 Å, while the bulk MgO lattice constant is 4.21 Å; the MgO is then under about +0.6% tensile strain. At s=+0.5%, the MgO strain reaches roughly +1.1%. A thick MgO substrate cannot be strained this way, so the strain path studied does not correspond to a physically achievable epitaxial strain for FeRh on a thick MgO substrate. Because the O atoms at the interface belong to the strained MgO, the Fe-O distance and hybridization change with s, and the reported MCA sign change may be driven by substrate strain rather than by intrinsic FeRh d-band shifts. The natural FeRh/MgO mismatch (about -0.6% on FeRh) lies at the edge of the studied range, and the predicted spin reorientation near s=0 is not the natural interface condition. Furthermore, Table I and Eq. (1), which contain the MCA values and the perturbation-theory expression central to the claim, are missing from the manuscript. Without disentangling FeRh strain from substrate strain, the conclusion that MCA in FeRh/MgO is strain-controllable is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles DFT (VASP, PBE, PAW) calculations of the magnetocrystalline anisotropy (MCA) of a five-monolayer FeRh(001) film on an MgO(001) substrate, for biaxial strain s in the range -0.5% to +0.5%. The authors claim that the easy axis switches from in-plane under compressive strain to out-of-plane under tensile strain, and that the effect originates from strain-induced shifts of the interfacial Fe d-states (dx2-y2 and dxz,yz), analyzed via second-order perturbation theory and k-resolved band character. They also report that the G-type antiferromagnetic phase remains more stable than the ferromagnetic phase by about 20 meV/Fe over the strain range.","tokens_in":4920,"tokens_out":7908,"duration_ms":116513,"significance":"If the central claim holds, the paper would provide a concrete strain-engineering route for controlling the magnetic anisotropy of antiferromagnetic FeRh films, which is of interest for AFM spintronics. The study uses standard DFT methods, contains no fitted parameters, and offers a specific orbital-level mechanism. However, the manuscript as provided is incomplete: the numerical MCA data (Table I), the perturbation-theory expression (Eq. 1), and the supporting band-structure/PDOS figures (Figs. 2-3) are referenced but not present, and the slab model applies strain to the entire FeRh/MgO stack, which may not correspond to a physically achievable epitaxial condition on a thick MgO substrate. These issues prevent verification of the central claims in the current version.","major_comments":[{"comment":"The central quantitative claim — that MCA changes sign between s=-0.5% and s=+0.5% — is not verifiable from the provided text, because Table I (with the MCA values and moments) and Fig. 2 (the k-resolved orbital character) are referenced but their content is absent. The paper should include the actual MCA values (with units of erg/cm2) for each strain and each magnetic phase, and the figure should be present.","section":"Section 3, Table I and Fig. 2"},{"comment":"The slab model uses a single in-plane lattice constant for the entire FeRh/MgO supercell, and the strain s is applied to the whole stack. At s=0, the FeRh in-plane constant is 2.995√2 = 4.235 Å, while bulk MgO is ~4.21 Å, so the MgO substrate is already at +0.6% tensile strain; at s=+0.5%, the MgO strain is ~+1.1%. A thick MgO substrate would not be strained in this way in an epitaxial FeRh/MgO bilayer. The natural FeRh/MgO mismatch is about -0.6% (FeRh in-plane compressed), which lies at the edge of the studied range. The authors should disentangle film strain from substrate strain, e.g., by fixing the MgO at its bulk lattice constant and varying the FeRh in-plane constant, or by explicitly justifying the co-strained model. Without this, the predicted spin reorientation near s=0 is not clearly established for real FeRh/MgO heterostructures.","section":"Section 2"},{"comment":"There is a direct contradiction about the direction of the easy-axis switching: the abstract states the switching is 'from perpendicular to in-plane', while the introduction and Section 3 state that compressive strain gives an in-plane easy axis and tensile strain gives an out-of-plane easy axis. Next to this, the text 'gives the negative MCA of 0.47 erg/cm2 for s=0.50%' appears inconsistent with the claimed sign of the effect; this likely should read 's=-0.50%' and '-0.47 erg/cm2'. The sign convention and the direction of switching must be made consistent throughout.","section":"Abstract, Introduction, and Section 3"},{"comment":"The perturbation-theory expression for the MCA is missing from the manuscript; the text refers to 'Eq. (1)' but the equation itself does not appear. Since the orbital-resolved analysis is entirely based on this expression, the equation must be provided.","section":"Section 3, Eq. (1)"}],"minor_comments":[{"comment":"Typographical errors: 'usingab-initio' (Abstract), 'antiferromagetic' (Abstract), 'expitaxially' (Introduction), 'the the' (Section 3), and 'exchange correlation' (Section 2) should be corrected.","section":"Throughout"},{"comment":"The sentence 'Note because the two Fe atoms on each atomic plane in the G-AFM phase have opposite ms.' is incomplete; the following 'we only list its magnitude' is orphaned. This should be rephrased.","section":"Section 3"},{"comment":"The definition 'The MCA per interfacial area is determined by [E[100]-E[001]]' should clarify that the energy difference is divided by the interfacial area, and define the sign convention clearly.","section":"Section 2"},{"comment":"The Fig. 2 caption contains a duplicated sentence about k points with large negative or positive contributions. In the Conclusion, the notation 'Lyz,zx' in the matrix element is unclear and should be defined consistently.","section":"Figure 2 caption and Conclusion"},{"comment":"The calculation considers only one interface termination (Fe-terminated, O atop Fe) and one film thickness (5 ML). The authors should comment on the expected robustness of the strain-induced switching with respect to termination and thickness, or state that the prediction is specific to this model.","section":"Section 2"},{"comment":"The phrase 'the spin reorientation occurs around 0' should specify units; it should read 'around s=0' (i.e., the lattice-matched condition in the model).","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to be an early or incomplete draft: the central data table, equation, and figures are missing from the text, which is unusual for a submission. The co-straining of the substrate is a serious concern that the authors will need to address with new calculations or a careful discussion. If the missing content is supplied and the substrate-strain issue is resolved, the paper could be suitable for publication. I also noticed the internal contradiction about the switching direction; it suggests the manuscript has not been carefully proofread."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent but incomplete DFT paper. The new piece is the systematic strain-dependent MCA calculation for an FeRh/MgO bilayer, with a band-resolved explanation in terms of SOC between Fe d-states. That is a reasonable extension of prior work on structural stability under strain (refs 16,17). The methodology is standard (PBE, PAW, thin slab) and the perturbation-theory analysis is post-hoc interpretation, not fitting. Significance is incremental — strain-controlled MCA is established in other systems — but for FeRh/MgO this specific new application is useful.\n\nBut the central claim cannot be checked from the text. Table I, Fig. 2, Fig. 3, and Eq. (1) are all referenced but missing. The abstract says the easy axis switches from perpendicular to in-plane, while the results section says compressive strain gives in-plane and tensile gives out-of-plane. One of those is backwards. That kind of inconsistency matters in a paper whose only quantitative output is the sign of MCA.\n\nThe bigger problem is the strain definition. The supercell has one in-plane lattice constant for FeRh and MgO. At s=0 the FeRh is at its bulk value and the MgO is stretched by +0.6%; at s=+0.5% the MgO is stretched by over 1%. A thick MgO substrate cannot be strained that way, so the calculated 'epitaxial strain' path is not physically realizable for FeRh/MgO. Since the O atoms at the interface belong to the strained MgO, the Fe-O distance changes with s, and the reported MCA switching could be driven by substrate strain rather than intrinsic FeRh d-band shifts. The natural FeRh/MgO mismatch sits at about -0.6% strain on FeRh, at the edge of the range. So the central prediction is unsupported unless the two components' strains are decoupled.\n\nThe paper does have useful pieces. The observation that the G-AFM phase is stable across the strain range and the strain dependence of the orbital moment differences is a sensible direction. The mechanism section, once the figures are restored, could be informative.\n\nMy recommendation: this is worth a serious referee, but not in its current form. The authors need to supply the missing table and equation, fix the direction language, and redo or reframe the strain calculations so that the film strain is varied while the substrate remains at its own equilibrium lattice constant. If that changes the result, so be it; the paper would at least be consistent.","headline":"Plausible but unverifiable DFT prediction of strain-driven MCA switching in FeRh/MgO, undermined by missing data and a strain model that co-strains the substrate.","tokens_in":5488,"tokens_out":4149,"would_cite":false,"duration_ms":36841,"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","75.70.-i"],"model":"deepseek-v4-flash","headline":"A thin FeRh film on MgO can have its magnetic easy axis flipped by strain, because strain shifts the iron d-electron states that carry spin-orbit coupling.","keywords":["FeRh","magnetocrystalline anisotropy","epitaxial strain","spin-orbit coupling","antiferromagnetic spintronics","density functional theory","MgO substrate","easy axis reorientation"],"falsifier":"A direct DFT calculation of the same five-monolayer FeRh/MgO slab with the opposite (Rh-terminated) interface, or with full in-plane relaxation, that fails to show a magnetocrystalline anisotropy sign change between -0.5% and +0.5% strain would falsify the claim; so would an experimental torque or Kerr measurement on epitaxial FeRh/MgO under biaxial strain showing no easy-axis rotation in that strain window.","tokens_in":4418,"feed_emoji":"🧲","tokens_out":8320,"duration_ms":81561,"temperature":0.7,"pith_summary":"This paper uses first-principles electronic structure calculations to show that epitaxial strain can reorient the magnetization direction of an antiferromagnetic FeRh thin film grown on MgO. In the antiferromagnetic phase, compressing the film by half a percent makes the magnetic easy axis lie in the film plane, while stretching it by half a percent makes the easy axis stand out of plane, with the switch occurring near zero strain. The authors trace this to spin-orbit coupling between specific iron d-orbitals at the Fe/MgO interface, whose energies shift relative to the Fermi level as strain changes. If true, this gives a practical, substrate-driven knob for tuning magnetic anisotropy in antiferromagnetic spintronic devices.","feed_headline":"Strain switches FeRh magnetism from in-plane to out-of-plane","feed_subtitle":"Compressing or stretching the FeRh/MgO interface shifts iron orbital states and flips the easy axis.","key_machinery":"The load-bearing object is the second-order perturbation theory expression for the magnetocrystalline anisotropy, evaluated from the spin-orbit coupling between occupied and unoccupied states of the interfacial Fe atoms. In the slab geometry, the relevant channel is the coupling matrix element $\\langle d_{xz,yz}|L_{xy}|d_{x^2-y^2}\\rangle$ between iron $d$ orbitals; strain shifts these states relative to the Fermi level at particular $k$-points, converting a net negative anisotropy into a positive one. The slab supercell itself, with five FeRh monolayers on a MgO substrate and only atomic $z$-positions relaxed, is what makes the epitaxial strain well-defined.","core_discovery":"The central claim is that the magnetocrystalline anisotropy of an ultrathin FeRh/MgO bilayer changes sign around zero epitaxial strain, so that the magnetic easy axis rotates from in-plane under compression to out-of-plane under tension. The sign change is driven by strain-induced shifts of the spin-orbit-coupled $d$ states of the interfacial Fe atoms: at compressive strain the $\\langle d_{xz,yz}|L_{xy}|d_{x^2-y^2}\\rangle$ coupling dominates with a negative contribution to the anisotropy energy, while under tensile strain band shifts suppress that negative channel and enable a positive one, reversing the easy axis. The authors further find that the G-type antiferromagnetic phase remains more stable than the ferromagnetic phase by about 20 meV per Fe atom across the strain range considered, so the reorientation is an antiferromagnetic effect. The claim is specifically for the slab model used: five FeRh monolayers on MgO(001) with an Fe-terminated interface and oxygen atoms sitting atop iron.","pith_inferences":["The predicted sign of the strain effect likely depends on the Fe/O interfacial registry; a Rh-terminated interface or a different oxygen position could suppress or reverse the reorientation, which is a testable extension of the same calculation.","Because the reorientation is an interfacial effect, increasing film thickness should dilute the anisotropy change, so the practical strain window should shrink for films much thicker than five monolayers.","On piezoelectric substrates, dynamic strain could offer an in-operando switch of the easy axis, going beyond the static strains studied here."],"forward_implications":["In the G-AFM phase, the easy axis can be selected between in-plane and out-of-plane by choosing compressive or tensile epitaxial strain, without changing the magnetic order.","The interface Fe atoms, not the surface Fe atoms, dominate the anisotropy, so interface chemistry is the control parameter.","The G-AFM phase stays about 20 meV per Fe atom more stable than the FM phase over the +/-0.5% strain window, putting the switch in a robust antiferromagnetic state.","Strain engineering, including growth on lattice-mismatched or piezoelectric substrates, becomes a viable route to set magnetization direction in FeRh-based spintronic devices."],"supporting_citations":[{"why":"Shows FeRh films grown on MgO and piezoelectric substrates inherit lattice mismatch strain, the physical setting for the strain range studied.","marker":"[14,15]"},{"why":"Earlier ab initio calculations of the stability of FM and AFM FeRh structures under epitaxial strain, the baseline for the phase ordering found here.","marker":"[16,17]"},{"why":"Bulk FeRh under hydrostatic pressure, the bulk electronic-structure baseline the slab calculation extends to surfaces and interfaces.","marker":"[13]"},{"why":"Supplies the projector augmented wave method and the generalized-gradient approximation used in every calculation.","marker":"[18-22]"}],"fun_headline_variants":["Strain flips FeRh easy axis between in-plane and out-of-plane","Strain engineering reorients FeRh/MgO magnetism","Strain tunes FeRh anisotropy and reverses easy axis","Strain toggles FeRh/MgO magnetic easy axis","FeRh/MgO easy axis flips under strain"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole prediction rests on the slab model where FeRh is locked to the MgO lattice with an Fe-terminated interface and only atomic z-positions relax; if the real interface termination, film thickness, or strain relaxation differs, the strain-driven easy-axis switch could disappear or reverse.","fun_headline_variants_meta":{"raw":{"variants":["Strain flips FeRh easy axis between in-plane and out-of-plane","Strain engineering reorients FeRh/MgO magnetism","Strain tunes FeRh anisotropy and reverses easy axis","Strain toggles FeRh/MgO magnetic easy axis","FeRh/MgO easy axis flips under strain"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001114,"raw_usage":{"total_tokens":4619,"prompt_tokens":908,"completion_tokens":3711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":3624}},"tokens_in":524,"tokens_out":3711,"duration_ms":23446,"temperature":1.0,"reasoning_tokens":3624,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:56:43.765166+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct DFT calculation of the same five-monolayer FeRh/MgO slab with the opposite (Rh-terminated) interface, or with full in-plane relaxation, that fails to show a magnetocrystalline anisotropy sign change between -0.5% and +0.5% strain would falsify the claim; so would an experimental torque or Kerr measurement on epitaxial FeRh/MgO under biaxial strain showing no easy-axis rotation in that strain window.","supporting_citations":[],"review_version":1}