{"id":"f8c8c99e-5e4c-4515-9bc5-1f6777961d94","arxiv_id":"1908.02332","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Spin-wave measurements show that MnBi2Te4 has frustrated intralayer magnetic exchange close to the classical limit for ferromagnetism.","lead":"Inelastic neutron scattering on the antiferromagnetic topological insulator MnBi2Te4 reveals frustrated magnetic interactions within its manganese layers, with competing ferromagnetic and antiferromagnetic exchange close to a stability boundary. The result suggests that chemically tuning this material could drive it into exotic spin textures such as skyrmions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The |J2/J1|=0.36 result depends on a uniform 0.85 meV broadening that the paper itself says cannot describe the full spectrum; until Q-dependent broadening is tested, the quantitative proximity to 1/3 is not secure.","rationale":"The reader's weakest-assumption analysis correctly identified the uniform, Q-independent broadening as the critical assumption behind |J2/J1| = 0.36. My reading confirms this concern and strengthens it by noting that the paper itself reports Model-c and Model-ab cannot be reconciled, that Model-c fixes SJ2 = 0 while Model-ab gives SJ2 = -0.10 meV, and that a single Heisenberg parameter set is explicitly insufficient. The broadening width (0.85 meV) exceeds the fitted nearest-neighbor exchange (0.28 meV), meaning the spin-wave peaks are heavily overlapping and the fit is not a clean measurement of the dispersion. Still, the qualitative evidence for frustration is credible: the data show broad, renormalized magnetic spectral weight and the paper is admirably transparent about its limitations. A single-crystal INS experiment or a Q-dependent broadening re-fit would settle whether |J2/J1| is really near 1/3 or merely consistent with it within large systematic uncertainty. Since the reader's verdict was already CONDITIONAL, my assessment does not change it; the paper should remain conditionally accepted pending sharper experimental constraints.","tokens_in":7922,"tokens_out":3802,"duration_ms":43090,"concrete_test":"Re-fit Model-ab to the same powder data with the broadening width allowed to vary with Q (e.g., piecewise-constant in two or three Q windows) or by separately fitting the low-Q (0.8-1.3 Å^-1) and high-Q (1.3-1.9 Å^-1) portions. If the best-fit |J2/J1| changes by more than about 0.1, or if the Q-dependent model substantially improves the residuals, then the uniform-broadening assumption is load-bearing and the central ratio is not robust until single-crystal INS resolves the intrinsic linewidth.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline ratio rests entirely on Model-ab, which fits the powder-averaged high-energy spectrum under the assumption of a single, Q-independent Lorentzian broadening of 0.85 meV. The paper explicitly flags an internal inconsistency: 'the sharpness of intralayer modes (at the gap edge) and the broad, high energy interlayer modes cannot be consistently modeled with a single set of Heisenberg parameters' (Model-c vs Model-ab). It later concedes 'clear evidence for strongly Q-dependent broadening,' and notes that low-Q data display a high-energy tail that Model-ab cannot reproduce. With FWHM = 0.85 meV larger than SJ1 = 0.28 meV and roughly three times the exchange energy scale, individual spin-wave branches are unresolved, so the extracted (SJ1, SJ2) depend sensitively on the assumed line shape, background, and fitted Q range. The problem is not just statistical: Model-c fixes SJ2 = 0 while Model-ab gives SJ2 = -0.10 meV, and the two models use different anisotropy (SD = 0.12 vs 0.55 meV). If the true broadening is Q- or energy-dependent, the inferred |J2/J1| = 0.36 could shift substantially; the claim that the system sits near the classical instability limit |J2/J1| = 1/3 is therefore quantitatively fragile, even though the qualitative case for frustration is plausible.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports powder inelastic neutron scattering (INS) measurements on the putative antiferromagnetic topological insulator MnBi2Te4. The authors model the spin dynamics with a Heisenberg Hamiltonian containing intralayer nearest- and next-nearest-neighbor exchanges J1 and J2, interlayer exchange Jc, and single-ion anisotropy D. Because one parameter set cannot simultaneously describe the low-energy gap edge and the high-energy intralayer modes, they construct two models: Model-c (low-energy gap edge, with resolution-limited broadening) and Model-ab (high-energy intralayer spectrum, with a uniform 0.85 meV Lorentzian broadening). The central claim is that |J2/J1| = 0.36, close to the classical instability limit |J2/J1| = 1/3 for ferromagnetic triangular layers. This is supported by DFT+U calculations at moderate U (about 2.7 eV) and by classical Monte Carlo simulations that show nearby spiral and skyrmion phases. The paper also emphasizes Ising-like anisotropy and metamagnetism.","tokens_in":8231,"tokens_out":5133,"duration_ms":53291,"significance":"If the extracted ratio |J2/J1| is robust, the result is significant: it would place MnBi2Te4 near a magnetic instability where small perturbations (chemical substitution, strain, field) could drive the system into helimagnetic or skyrmionic states, which would matter for quantum anomalous Hall and axion-insulator physics. The paper has clear strengths: it provides a quantitative two-model analysis of powder INS data, uses magnetization critical fields as independent constraints for the gap-edge model, performs DFT+U calculations over a range of U, and is unusually candid about inconsistencies between models. However, because the headline ratio rests entirely on a model whose broadening assumption the authors themselves acknowledge is inadequate, the numerical proximity to 1/3 is not yet established with the claimed precision.","major_comments":[{"comment":"The extraction of the headline ratio |J2/J1| = 0.36 depends exclusively on Model-ab, which convolutes calculated spectra with a single, Q-independent Lorentzian of FWHM = 0.85 meV. Because this width is larger than SJ1 = 0.28 meV, the individual spin-wave branches are not resolved in the powder-averaged data. The manuscript itself states that the gap-edge and high-energy modes cannot be modeled with a single set of Heisenberg parameters, and later concludes with 'clear evidence for strongly Q-dependent broadening.' Under these conditions the fitted (SJ1, SJ2) values are sensitive to the assumed line shape, energy-dependent background, and fitted Q range, so the numerical proximity to 1/3 is not secure. This is load-bearing because the paper's main claim is precisely that the system sits close to the classical instability limit. I request a robustness analysis: fits with Q- or energy-dependent broadening, alternative background models, and fits excluding the low-Q data that show the high-energy tail, to demonstrate that the ratio is stable or to quote a conservative range.","section":"Model-ab fitting paragraph and Table I"},{"comment":"The two independent parameter sets imply inconsistent single-ion anisotropies. Model-ab gives SD = 0.55 meV and an effective spin gap of 1.1 meV, while Model-c and the magnetization data give SD about 0.085–0.12 meV and the true gap of about 0.5 meV. The paper acknowledges that this discrepancy is unexplained. Nevertheless, the Monte Carlo phase diagram is computed at D/J1 = 0.4, obtained by combining SD from Model-c with SJ1 from Model-ab, and the text quotes 0.3 < D/J1 < 0.4 as an experimental finding. No single model in Table I yields D/J1 = 0.4; combining values from inconsistent models is an ad hoc procedure. The skyrmion phase in Fig. 3(a) appears only for |J2/J1| greater than about 0.5 at this anisotropy, whereas the fitted ratio is 0.36, so the proximity statement involves a substantial extrapolation. The authors should clearly label the Monte Carlo results as scenario calculations and present the phase diagram over the full range of D/J1 allowed by the data, including the Model-ab value D/J1 approximately equal to 2.","section":"Table I; 'Overall, our findings' paragraph"},{"comment":"The classical Monte Carlo phase diagram is computed using parameters obtained from the same INS data, so Fig. 3 is not an independent confirmation of skyrmion proximity; it is a scenario calculation whose input parameters are the same uncertain quantities that are being tested. This is not a statistical circularity in the strict sense, but the text's claim that skyrmion phases 'could be accessed in chemically tuned compounds' overstates the predictive content, especially given that the model places the skyrmion phase outside the fitted |J2/J1| range at D/J1 = 0.4. I recommend explicit language stating that the skyrmion and spiral phases occur for parameter ratios beyond those inferred from the current fits, and that the quoted proximity to the 1/3 instability is therefore a motivation for further study rather than a quantitative prediction.","section":"Monte Carlo simulation paragraph and Fig. 3"}],"minor_comments":[{"comment":"There are several typographical errors: 'teteradymite' should be 'tetradymite' (appearing twice), and 'Componds' should be 'Compounds' in the final section before the acknowledgments.","section":"Abstract and Introduction"},{"comment":"Reference [17] appears only as 'Supplementary Material reference'; the authors should include a complete citation to the supplementary material.","section":"References"},{"comment":"The comparison of Model-ab to the data at low Q is described as 'less satisfactory,' but the text does not quantify this disagreement or explain whether the fitting Q range (0.8–1.9 Å−1) was chosen to exclude the problematic low-Q region. A brief statement on how the Q range was selected would help the reader assess the fit's sensitivity.","section":"Fig. 1(c)"},{"comment":"For clarity, the sign convention and the parameterization in terms of S times the coupling constants should be stated explicitly in the text near Eq. (1); the table lists SJ1, SJ2, SJc, and SD, but the conversion from these products to the Hamiltonian parameters J1, J2, Jc, and D is not spelled out.","section":"Eq. (1)"},{"comment":"The sentence 'Recent first-principles electronic structure calculations with U = 5 eV predict that |J2/J1| < 0.03' would be easier to check if the authors explicitly contrasted this with their own DFT+U result at 5 eV in Table I, which gives |J2/J1| approximately 0.07; the difference likely arises from details of the method, and a one-sentence explanation would prevent confusion.","section":"DFT+U discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful contribution to the magnetic characterization of MnBi2Te4, and the candid discussion of the model inconsistency is commendable. However, the central quantitative claim (|J2/J1| = 0.36, proximity to 1/3) rests on a broadening model that the authors themselves call into question, and the Monte Carlo phase diagram is computed from an ad hoc combination of parameters from the two inconsistent models. I would support publication after the authors either provide robustness tests for the extracted ratio and anisotropy or substantially soften the claims, reframing the result as a qualitative indication of frustration. The paper is not fatally flawed and should not be rejected outright."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first inelastic neutron scattering study of spin waves in MnBi2Te4, and it gives a plausible qualitative case that the triangular Mn layers are close to magnetic frustration. The authors extract |J2/J1| = 0.36, near the classical instability at 1/3, and show that earlier DFT predictions of near-zero J2 are wrong. That is a real result with consequences for the quantum anomalous Hall and axion insulator discussions.\n\nThe paper does several things well. The powder INS data are carefully presented; the magnetization data on single crystals provide complementary constraints; and the authors are unusually candid about the limits of their model. They split the analysis into Model-c for the gap edge and Model-ab for the high-energy spectrum, and they note explicitly that a single Heisenberg parameter set cannot describe both. They also flag \"clear evidence for strongly Q-dependent broadening\" from powder data.\n\nThe soft spot is exactly that split. Model-ab, which determines J2/J1, assumes a uniform 0.85 meV Lorentzian broadening, larger than SJ1 = 0.28 meV. With individual branches unresolved, the fitted (SJ1, SJ2) depend on line shape, background, and Q range. The internal inconsistency is real: Model-c uses SD = 0.12 meV and SJ2 = 0, while Model-ab uses SD = 0.55 meV and SJ2 = -0.10 meV. The 0.85 meV broadening is the load-bearing assumption, and the paper admits it is an oversimplification. So the quantitative proximity to 1/3 is fragile; the qualitative frustration conclusion is probably right.\n\nTwo smaller issues. The DFT+U section tunes U to about 2.7 eV to match the data; that is consistency checking, not independent support. The classical Monte Carlo phase diagram uses the fitted parameters, so the skyrmion proximity statement is an extrapolation, not a prediction. The paper is honest about both, so I do not read them as deceptions.\n\nBottom line: this is a useful experimental contribution for anyone working on MnBi2Te4 or frustrated triangular magnets. It deserves a serious referee; I would send it to review. But the referee should demand a careful statement that the J2/J1 value is model-dependent, ideally with a single-crystal INS follow-up or at least an explicit test of Q-dependent broadening before it becomes a textbook number.","headline":"First INS spin-wave data on MnBi2Te4 make a credible case for frustrated intralayer exchange, but the headline |J2/J1|=0.36 rests on a uniform-broadening fit the paper itself contradicts.","tokens_in":8847,"tokens_out":2224,"would_cite":true,"duration_ms":22810,"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":"Inelastic neutron scattering on MnBi2Te4 finds intralayer frustration $|J_2/J_1| = 0.36$, placing the ferromagnetic triangular layers near the classical instability boundary at $1/3$ where helical and skyrmion phases compete.","keywords":["MnBi2Te4","antiferromagnetic topological insulator","inelastic neutron scattering","spin frustration","triangular lattice","J1-J2 Heisenberg model","spin waves","skyrmion phases"],"falsifier":"A single-crystal inelastic neutron scattering experiment that resolves the intralayer spin-wave branches along high-symmetry directions would settle it: if a direct fit gives $|J_2/J_1|$ clearly below $1/3$, or if the line shapes force a strongly momentum-dependent broadening instead of the uniform 0.85 meV assumed here, the proximity-to-instability claim fails.","tokens_in":7710,"feed_emoji":"🧲","tokens_out":18825,"duration_ms":174534,"temperature":0.7,"pith_summary":"The paper uses inelastic neutron scattering to map the spin waves of MnBi2Te4, a candidate antiferromagnetic topological insulator built from ferromagnetic Mn triangular layers with antiferromagnetic interlayer coupling. The measurements show that the intralayer magnetism is strongly frustrated: the next-nearest-neighbor exchange is antiferromagnetic, with $|J_2/J_1| = 0.36$, close to the classical threshold $1/3$ at which ferromagnetic ordering inside a triangular layer becomes unstable. The paper argues that this places the material near competing magnetic phases, so chemical doping, magnetic fields, or strain could plausibly drive it into helical, skyrmion, or other noncollinear states. That matters because the topological electronic states live in the adjacent Bi-Te layers, so switching the magnetic texture could switch topological transport behavior. DFT+U and classical Monte Carlo calculations are used to reinforce the same picture, with the measured exchange ratios reproduced at moderate correlation strength and spiral and skyrmion phases appearing nearby in parameter space.","feed_headline":"0.36 frustration ratio puts MnBi2Te4 at a magnetic tipping point","feed_subtitle":"At |J2/J1|=0.36, a modest field or chemical doping could push the layers into helical or skyrmion phases.","key_machinery":"The central object is a local-moment Heisenberg Hamiltonian on the stacked triangular lattice with intralayer nearest-neighbor exchange $J_1$, intralayer next-nearest-neighbor exchange $J_2$, interlayer exchange $J_c$, and single-ion anisotropy $D$ for the $S=5/2$ Mn moments. The load-bearing quantity is the dimensionless frustration ratio $|J_2/J_1|$ compared with the classical boundary $1/3$ between intralayer ferromagnetic and noncollinear order. The parameters are extracted by comparing powder-averaged linear spin-wave intensity calculations to the measured $S(Q,E)$, using separate fits for the low-energy gap edge (Model-c) and the full intralayer spectrum (Model-ab). Classical Monte Carlo simulations then map the fitted parameters into a field-temperature phase diagram that places the spiral and skyrmion phases nearby.","core_discovery":"The paper's central claim is that the magnetic interactions in MnBi2Te4 form a frustrated $J_1$-$J_2$ triangular-lattice system rather than a simple nearest-neighbor ferromagnet. From the powder-averaged spin-wave spectrum the authors extract $SJ_1 = 0.28(2)$ meV and $SJ_2 = -0.10(2)$ meV, giving $|J_2/J_1| = 0.36$, while the low-energy gap edge and magnetization data yield an interlayer coupling $SJ_c$ near $-0.055$ to $-0.085$ meV and a single-ion anisotropy $SD$ near 0.12 meV. The excitations show a spin gap of about 0.5 meV and an intrinsic linewidth that requires convolution by roughly 0.85 meV, far beyond the instrumental resolution. The high-energy intralayer modes need an effective anisotropy closer to $SD = 0.55$ meV, a discrepancy the paper attributes to momentum-dependent broadening that powder data cannot resolve. The major finding is that $|J_2/J_1|$ sits close to the classical instability limit $1/3$ for ferromagnetic triangular layers, so the ground state is nearly degenerate with helical and multi-$q$ states; classical Monte Carlo on the fitted parameters finds vertical spiral, skyrmion, and up-up-down-down phases at nearby points in the field-frustration phase diagram, and DFT+U reproduces the measured exchange ratios only near $U = 2.7$ eV.","pith_inferences":["The large gap between the low-energy anisotropy ($SD$ near 0.12 meV) and the high-energy effective value ($SD$ near 0.55 meV) suggests that the high-energy fit may be absorbing momentum-dependent broadening; resolving the linewidths on a single crystal would separate those effects.","If the frustration ratio is really this close to 1/3, strain, surface termination, or small off-stoichiometry should also be able to tip the surface layers into noncollinear textures, which could connect the bulk frustration to surface-sensitive measurements.","A direct extension would be to map the field-temperature phase diagram of Mn(Bi,Sb)2Te4, where the paper notes the anisotropy is smaller; skyrmion scattering should appear at accessible fields if the proximity argument is correct."],"forward_implications":["The measured frustration places MnBi2Te4 near the classical $|J_2/J_1|=1/3$ boundary, so small changes in exchange, anisotropy, or field could tip the triangular layers into a different magnetic state.","The strong intrinsic line broadening seen in the spin-wave spectrum points to magnon coupling to other degrees of freedom, which single-crystal neutron studies could directly characterize.","The DFT+U comparison fixes the effective correlation strength near $U=2.7$ eV, so the neutron data provide a quantitative anchor for electronic-structure calculations of this topological magnet.","The nearby spiral, skyrmion, and up-up-down-down phases in the classical phase diagram give concrete search targets for chemical substitution, particularly in Sb-substituted compounds where the anisotropy is smaller."],"supporting_citations":[{"why":"Provides the sample synthesis, the A-type antiferromagnetic order determination, and the magnetization critical fields that anchor the low-energy Model-c analysis.","marker":"[10]"},{"why":"Supplies the classical stability condition $|J_2/J_1| < 1/3$ for ferromagnetic triangular layers, the boundary the measured ratio is compared against.","marker":"[14]"},{"why":"Extends the triangular-lattice $J_1$-$J_2$ phase diagram and reinforces the same $1/3$ instability line used as the theoretical reference.","marker":"[15]"},{"why":"Supplies the numerical procedure for powder-averaged spin-wave intensity calculations used in both the gap-edge and intralayer model fits.","marker":"[19]"},{"why":"Provides single-crystal spin-wave results on hexagonal MnTe, whose similar Mn triangular layers give a comparable nearest-neighbor exchange used to benchmark $SJ_1$.","marker":"[20]"},{"why":"Earlier first-principles prediction of the magnetic and topological properties of MnBi2Te4 that frames the paper's experimental test of the exchange parameters.","marker":"[7]"},{"why":"Prior DFT+U calculation with U=5 eV predicting $|J_2/J_1|<0.03$; the neutron data contradict this and motivate the U-dependent re-analysis.","marker":"[9]"},{"why":"Supplies the four-state energy-mapping method used to extract Heisenberg exchange parameters from DFT+U total energies.","marker":"[22]"},{"why":"Supplies the DFT+U formulation used in the first-principles calculation of the magnetic interactions.","marker":"[23]"},{"why":"Shows that frustrated triangular magnets can host skyrmion phases, providing the motivation for expecting skyrmion states near the measured frustration.","marker":"[16]"}],"fun_headline_variants":["MnBi2Te4 magnetic interactions poised at frustration threshold","Spin waves reveal near-critical frustration in MnBi2Te4","MnBi2Te4 teeters on magnetism's edge with J2/J1=0.36","Frustrated spins in MnBi2Te4 close to magnetic tipping point"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fitted ratio $|J_2/J_1|=0.36$ assumes a single Heisenberg model with a uniform, momentum-independent spin-wave lifetime broadening of 0.85 meV; if the broadening varies with momentum or additional couplings exist, the ratio shifts.","fun_headline_variants_meta":{"raw":{"variants":["MnBi2Te4 magnetic interactions poised at frustration threshold","Spin waves reveal near-critical frustration in MnBi2Te4","MnBi2Te4 teeters on magnetism's edge with J2/J1=0.36","Frustrated spins in MnBi2Te4 close to magnetic tipping point"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000541,"raw_usage":{"total_tokens":2622,"prompt_tokens":1006,"completion_tokens":1616,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":1535}},"tokens_in":622,"tokens_out":1616,"duration_ms":10958,"temperature":1.0,"reasoning_tokens":1535,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:47:06.024062+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single-crystal inelastic neutron scattering experiment that resolves the intralayer spin-wave branches along high-symmetry directions would settle it: if a direct fit gives $|J_2/J_1|$ clearly below $1/3$, or if the line shapes force a strongly momentum-dependent broadening instead of the uniform 0.85 meV assumed here, the proximity-to-instability claim fails.","supporting_citations":[{"cited_title":"Tanaka and N","cited_arxiv_id":null,"evidence_quote":"Extends the triangular-lattice $J_1$-$J_2$ phase diagram and reinforces the same $1/3$ instability line used as the theoretical reference."},{"cited_title":"Vaknin, D","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical procedure for powder-averaged spin-wave intensity calculations used in both the gap-edge and intralayer model fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides single-crystal spin-wave results on hexagonal MnTe, whose similar Mn triangular layers give a comparable nearest-neighbor exchange used to benchmark $SJ_1$."},{"cited_title":"Xiang, C","cited_arxiv_id":null,"evidence_quote":"Supplies the DFT+U formulation used in the first-principles calculation of the magnetic interactions."},{"cited_title":"Murao, F","cited_arxiv_id":null,"evidence_quote":"Shows that frustrated triangular magnets can host skyrmion phases, providing the motivation for expecting skyrmion states near the measured frustration."}],"review_version":1}