{"id":"12349ab1-fdff-4fd3-85dd-9d77398c1688","arxiv_id":"2507.18904","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"DFT analysis of GdRu2X2 (X = Si, Ge, Sn) links increased destabilizing band energy and competing exchange interactions to easier skyrmion formation, predicting GdRu2Sn2 as a high-temperature, low-field skyrmion host.","lead":"Using density functional theory calculations, this paper compares chemical bonding and electronic instabilities in three related magnetic compounds, GdRu2X2 with X = Si, Ge, and Sn. It proposes that stronger destabilizing energy and competing magnetic interactions predict skyrmion formation at higher temperatures and lower magnetic fields.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eband-based skyrmion descriptor is underdetermined: Fig. 11 has only two experimental anchors, Eband and competing-exchange metrics co-vary, and the Sn point is hypothetical, so the causal role of Eband is not established.","rationale":"The manuscript is a careful DFT study of bonding, Fermi-surface nesting, and exchange interactions; the individual calculations (COHP/COBI, LRF, J-couplings) are standard and internally consistent. The stress-test concern sits at the inference level of Fig. 11: with two experimental points plus a hypothetical Sn phase, the Eband trend is underdetermined. The reader's weakest assumption names this same issue, and their CONDITIONAL verdict matches my assessment. I add two technical refinements: (i) several monotonic descriptors (J1, MAD-1, lattice volume, FSN multiplicity) co-vary with Eband, so the specific controlling variable is not identified; and (ii) Eband is an absolute integrated quantity without per-electron normalization or error estimates, making cross-compound comparison fragile. These considerations do not invalidate the paper's computational results, but they do mean the abstract's phrasing that the trend is 'correlated' or 'revealed' should be softened to a hypothesis. No change from CONDITIONAL is needed.","tokens_in":16063,"tokens_out":7659,"duration_ms":81364,"concrete_test":"Apply the identical QE+LOBSTER pipeline to two external benchmarks: one experimentally established centrosymmetric skyrmion host with known skyrmion T and H (e.g., Gd2PdSi3) and one non-skyrmion isostructural control (e.g., LaRu2Si2). Compute Eband from Eq. (2) with per-atom and per-valence-electron normalization, using the same k-mesh and projection settings as in the paper. If the external skyrmion host does not fall on the Fig. 11 trend, or if the non-skyrmion control falls inside the GdRu2Si2-GdRu2Ge2 window, then Eband is not a sufficient, general predictive descriptor and the central claim reduces to a two-point observation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the integrated density-of-energy band energy Eband in Eq. (2) is a predictive descriptor: a more destabilizing Eband implies skyrmion formation at higher T and lower H. The evidence in Fig. 11 is two experimental compounds (GdRu2Si2, GdRu2Ge2) plus a hypothetical, variable-cell-optimized GdRu2Sn2 (Table S1) that has not been synthesized. With only two experimental points, any quantity that changes monotonically across the series will reproduce the trend. The paper itself reports several such co-varying quantities: the a-axis expands 4.164 to 4.232 to 4.311 Å (Table S1), J1 changes from -0.7 K to 71.1/76.2 K while J2 stays near -153 K (Table 1), MAD-1 increases 0.015 to 0.019 to 0.020, and the FSN pattern gains additional Q vectors. Fig. 11 even plots MAD-1 as shading, so the attribution to Eband rather than to frustration metrics is not identified. In addition, Eband is an integrated band-energy-like sum obtained from a LOBSTER projection; its absolute scale depends on the PAW/pseudopotential and projection settings, and no per-atom, per-electron, or reference-normalized values or error estimates are reported. A more negative Eband may therefore be a projection artifact rather than a physical instability. Because the third anchor is hypothetical, the correlation is a two-point interpolation, not a validated predictive law. The CONDITIONAL verdict is correct, but the abstract should not state the trend as established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the centrosymmetric skyrmion hosts GdRu2X2 (X = Si, Ge, Sn) using DFT-based chemical bonding analysis (COHP, COBI, MO diagrams), Fermi surface nesting through the Lindhard response function, density-of-energy (DOE) analysis, and energy-mapping of exchange interactions J1-J3. It proposes that a more destabilizing integrated band energy Eband (Eq. 2) and stronger competition between ferromagnetic and antiferromagnetic exchange interactions correlate with skyrmion formation at higher temperature and lower magnetic field, as summarized in Fig. 11. The experimental anchors are GdRu2Si2 and GdRu2Ge2, while GdRu2Sn2 is a hypothetical, variable-cell-optimized compound used as a predicted third point. The authors explicitly acknowledge that the Sn compound has not been synthesized and that the proposed trend needs experimental verification.","tokens_in":16501,"tokens_out":2903,"duration_ms":32205,"significance":"If validated, the proposed Eband-based descriptor would offer a computationally accessible chemical guideline for designing centrosymmetric skyrmion materials, a goal of current interest. The work is systematic in combining bonding indicators, Fermi-surface nesting, and exchange-parameter trends across a chemically coherent series, and it is transparent about the hypothetical nature of the Sn member. The main value is the formulation of a falsifiable design hypothesis rather than an established predictive law; in its present form the correlation rests on only two experimental compounds and one computed point, and the causal role of Eband is not isolated from other co-varying quantities.","major_comments":[{"comment":"The central descriptor Eband is an integrated density-of-energy value obtained from LOBSTER projections, but no numerical Eband values, no per-atom or per-electron normalization, and no uncertainty estimates are reported. Since the magnitude of such integrated quantities depends on the choice of projection basis and PAW potentials, the reader cannot assess whether the Si-Ge-Sn ordering is robust or a projection artifact. Please report the actual Eband numbers, test sensitivity to the projection settings, and indicate how the integration was truncated or normalized.","section":"Eq. (2), Fig. 9, and Fig. 11"},{"comment":"The proposed correlation is underdetermined because Eband co-varies with several other monotonic quantities across the series, including the a-axis length (4.164 to 4.232 to 4.311 Å), J1 (-0.7 to 71.1 to 76.2 K), MAD-1 (0.015 to 0.019 to 0.020), and the number of Fermi-surface nesting vectors. With only two experimental anchors, any monotonic descriptor would reproduce the Si-Ge trend, so the specific attribution to Eband is not demonstrated. The manuscript should discuss this co-variance explicitly and, ideally, provide a test that separates Eband from lattice-expansion or frustration metrics, or acknowledge in the main text that Eband is one of several interchangeable indicators.","section":"Fig. 11 and Table 1/Table S1"},{"comment":"The third point on the trend, GdRu2Sn2, is a variable-cell-optimized hypothetical structure with negative formation energy but no experimental realization, as stated in the text. Including this point in the same correlation plot as the two experimental compounds conflates a prediction with validation. The plot and the surrounding discussion should visually and textually distinguish experimental from predicted data, and the abstract/conclusions should not imply an established three-point law. The authors' own caveat that 'it warrants rigorous, experimental evidence for the Sn material' should be reflected more strongly in the presentation of Fig. 11.","section":"Fig. 11 and Table S1"}],"minor_comments":[{"comment":"The phrase 'Yoshimoshi et al.' in Section 3 should be 'Yoshimochi et al.' to match reference 40.","section":"General"},{"comment":"The abstract says 'These results reveal some correlation,' which is appropriately hedged, but the conclusions state the proposed trend more strongly. Please align the strength of the claims in the conclusions with the two-point nature of the evidence.","section":"Abstract and Conclusions"},{"comment":"The sentence preceding Eq. (3) says 'the following equations:5' with the reference number appearing after a colon; this is an odd citation placement and should be corrected.","section":"Eq. (3)"},{"comment":"Figure 11 would benefit from a caption specifying which points are experimental and which are calculated, and from error bars or shaded regions reflecting the uncertainty in Eband and the experimental skyrmion boundaries.","section":"Figure captions"},{"comment":"The SI uses 'pseudo potential' and 'diffused' inconsistently; these should be 'pseudopotential' and 'diffuse' throughout.","section":"Supporting Information"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the DFT workflow is standard. The main risk is that the central Eband descriptor is not yet convincingly separated from other monotonic structure-property trends. The authors are candid about the hypothetical Sn point, which is a strength. I would encourage the editor to request a revision that either provides sensitivity analysis for Eband and a clear experimental/predicted distinction in the key figure, or explicitly reframes the paper as proposing a hypothesis rather than demonstrating the correlation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent, internally consistent DFT/bonding study of GdRu2X2. The genuinely new content is the systematic COHP/COBI/DOE comparison and the specific proposal that Eband tracks the T and H conditions for skyrmion formation. That proposal is a hypothesis, not a demonstrated law. There are only two experimental skyrmion hosts in the dataset (Si and Ge); the third point is a variable-cell-optimized, unsynthesized GdRu2Sn2. Any monotonic descriptor would fit the two-point trend, and several descriptors actually co-vary here: the a-axis, J1, MAD-1, and the number of FSN vectors. So Fig. 11 cannot identify Eband as the controlling variable.\n\nWhat is good: the workflow is standard and carefully described; the FSN/Lindhard analysis is a fair reproduction of known results for Si and Ge, with the Sn extension new; the exchange-constant mapping is clearly presented; and the paper is transparent, at the end and in the conclusions, that the Sn prediction needs experimental verification. The self-citation to ref. 41 is appropriate rather than inflationary. This is not a circularity problem—no fitted constants are fed back into the trend—it is an underdetermination problem.\n\nThe soft spots are in proportion: Eband is an integrated LOBSTER quantity with no error bars, no per-electron normalization, and no test of sensitivity to PAW/projection settings. A more negative number may partly be a projection artifact. The Hubbard U (6.7 eV) is fixed; that is minor because the trends are robust across standard settings, but a U = 6.0/7.4 eV check would cost little. The abstract and the phrase 'demonstrates' overstate what two points plus a prediction can support.\n\nI'd send this to peer review, not desk reject. A serious referee can ask the authors to (i) temper the abstract and conclusions, (ii) add explicit discussion of the co-varying descriptors, and (iii) ideally include at least one more experimental anchor or a projection-sensitivity analysis. The paper is useful for computational materials chemists and for experimentalists looking for a cheap screening descriptor; they should read the trend as a target to test on GdRu2Sn2 and related compounds.\n\nRecommendation: engage. Worth referee time, likely publishable after revision.","headline":"A solid DFT/bonding study whose central predictive trend is a two-point correlation; worth publishing but with the Eband claim dialed back from 'demonstrated' to 'hypothesis'.","tokens_in":16990,"tokens_out":3233,"would_cite":false,"duration_ms":36020,"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 computed band-energy integral predicts the temperature and field at which skyrmions appear in centrosymmetric GdRu2X2 magnets.","keywords":["magnetic skyrmions","GdRu2X2","chemical bonding","Fermi surface nesting","density of energy","RKKY exchange interactions","centrosymmetric magnets","electronic instability"],"falsifier":"Synthesize GdRu2Sn2, determine its crystal structure, and measure its magnetic phase diagram: if its skyrmion pocket appears at a lower temperature or a higher field than GdRu2Ge2, or if the compound forms in a different structure, the proposed $E_{\\mathrm{band}}$ ordering is disproved.","tokens_in":15877,"feed_emoji":"🧲","tokens_out":9958,"duration_ms":94300,"temperature":0.7,"pith_summary":"This paper tries to locate the electronic origin of skyrmion formation conditions in centrosymmetric (inversion-symmetric) magnets, taking the isostructural series GdRu2X2 (X = Si, Ge, Sn) as a model. It argues that as the X-p orbitals become more extended from Si-3p to Ge-4p to Sn-5p, the coupling between the localized Gd spins and the itinerant [Ru2X2] conduction layer improves, the integrated destabilizing band energy $E_{\\mathrm{band}}$ from a density-of-energy analysis grows, and ferromagnetic versus antiferromagnetic exchange competition in the Gd square lattice strengthens. The paper then proposes that these computed quantities set the skyrmion conditions: the larger $E_{\\mathrm{band}}$ and the stronger the competing interactions, the higher the temperature and the lower the magnetic field at which skyrmions appear. If that correlation holds, a standard electronic-structure calculation becomes a predictive screen for designing skyrmion materials.","feed_headline":"One band-energy number orders skyrmion temperatures and fields","feed_subtitle":"More extended p orbitals raise that energy, pushing skyrmions to higher temperature and lower field.","key_machinery":"The load-bearing object is the density-of-energy (DOE) function, whose integral up to the Fermi level defines the band energy $E_{\\mathrm{band}}=\\int_{-\\infty}^{E_F}\\mathrm{DOE}(E)\\,dE$; this single number converts the whole band structure's destabilizing antibonding character into one correlation variable. The Lindhard response function supplies the Fermi-surface nesting vectors that mark the underlying electronic instability, and an energy-mapping analysis over six spin-ordered supercells supplies the exchange couplings $J_1$, $J_2$, and $J_3$ whose FM/AFM competition is the second correlation variable. Together they turn chemical orbital extension into a quantitative prediction about skyrmion phase boundaries.","core_discovery":"On the paper's own terms, chemical substitution in GdRu2X2 tunes the electronic instability that precedes skyrmion formation. The projected bonding analysis finds antibonding Ru–Ru and Ru–X states near the Fermi level, and these destabilizing contributions accumulate in the density-of-energy integral $E_{\\mathrm{band}}=\\int_{-\\infty}^{E_F}\\mathrm{DOE}(E)\\,dE$, which increases from Si to Ge to Sn. Fermi-surface nesting evolves in step: GdRu2Si2 shows one vector $Q=(q,0,0)$, GdRu2Ge2 shows two inequivalent vectors $Q=(q,0,0)$ and $Q_A=(q,q,0)$, and GdRu2Sn2 shows multiple nesting vectors. Exchange-energy mapping gives a nearly unchanged antiferromagnetic $J_2$ near $-153$ K while $J_1$ changes from weakly negative to strongly ferromagnetic, so the FM/AFM competition intensifies along the series. The paper concludes that larger destabilizing energy and stronger competing exchange make skyrmion formation more accessible, matching the measured Si and Ge phase diagrams and predicting that the unsynthesized GdRu2Sn2 would host skyrmions at higher temperature and lower field than its siblings.","pith_inferences":["The proposed $E_{\\mathrm{band}}$ descriptor is not tested beyond this one family; a natural extension would be to compute the same density-of-energy integral for other centrosymmetric skyrmion hosts and see whether their measured skyrmion temperatures and fields follow the same ordering.","The GdRu2Sn2 leg of the trend rests on a hypothetical variable-cell-optimized structure, so the decisive test is experimental: a real GdRu2Sn2 that adopts a different structure, or that fails to show the predicted easier skyrmion window, would reduce the correlation to a two-point observation.","One could also test the correlation in reverse within this family by computing $E_{\\mathrm{band}}$ for additional X-site substitutions (for example, mixed or doped compositions) and checking whether the skyrmion window moves monotonically with the computed energy."],"forward_implications":["The measured shift of GdRu2Ge2's skyrmion pocket to higher temperature and lower field relative to GdRu2Si2 is attributed to its larger integrated destabilizing energy and stronger FM/AFM competition.","GdRu2Sn2 is predicted to be the easiest skyrmion host of the three, with the largest $E_{\\mathrm{band}}$ and multiple Fermi-surface nesting vectors.","The number and type of nesting vectors extracted from the Lindhard response provide a computational fingerprint that can be checked before a candidate material is synthesized.","Because $J_1$ changes strongly while $J_2$ stays nearly constant, chemical substitution on the X site offers a way to tune the competing-interaction balance without destroying the dominant antiferromagnetic coupling."],"supporting_citations":[{"why":"Provides the experimental skyrmion temperature-field windows for GdRu2Ge2 that anchor the comparison.","marker":"[41]"},{"why":"Reports the square skyrmion lattice and phase boundaries in GdRu2Si2.","marker":"[42]"},{"why":"Documents competing exchange interactions at inequivalent wavevectors in GdRu2Ge2, grounding the multi-Q nesting picture.","marker":"[40]"},{"why":"Links a nested Fermi-surface band to skyrmion formation in GdRu2Si2.","marker":"[51]"},{"why":"Introduces the density-of-energy function whose integral defines $E_{\\mathrm{band}}$.","marker":"[68]"},{"why":"Supplies the plane-wave density functional theory method for the electronic structure calculations.","marker":"[52]"},{"why":"Provides the local-orbital projection procedure behind the chemical bonding and density-of-energy analyses.","marker":"[54]"},{"why":"Underwrites the energy-mapping formalism used to extract the exchange couplings.","marker":"[69]"}],"fun_headline_variants":["Chemical bonding reshapes skyrmion phase space","Longer p orbitals push skyrmions to higher T","Antibonding states enable skyrmion tuning","GdRu2X2 links bonding to skyrmion fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole trend assumes that the computed integrated band energy $E_{\\mathrm{band}}$ accurately controls how easily skyrmions appear, and that the predicted crystal structure of the still-unsynthesized GdRu2Sn2 is the structure that would actually form.","fun_headline_variants_meta":{"raw":{"variants":["Chemical bonding reshapes skyrmion phase space","Longer p orbitals push skyrmions to higher T","Antibonding states enable skyrmion tuning","GdRu2X2 links bonding to skyrmion fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000339,"raw_usage":{"total_tokens":1965,"prompt_tokens":1134,"completion_tokens":831,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":750,"completion_tokens_details":{"reasoning_tokens":765}},"tokens_in":750,"tokens_out":831,"duration_ms":8410,"temperature":1.0,"reasoning_tokens":765,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:05:36.800533+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Synthesize GdRu2Sn2, determine its crystal structure, and measure its magnetic phase diagram: if its skyrmion pocket appears at a lower temperature or a higher field than GdRu2Ge2, or if the compound forms in a different structure, the proposed $E_{\\mathrm{band}}$ ordering is disproved.","supporting_citations":[],"review_version":1}