{"id":"5b8a8605-37d8-4078-808c-10c701107638","arxiv_id":"1908.05417","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A microscopic calculation shows that the quantum anomalous Hall ferromagnet in twisted bilayer graphene is stable against spin and valley magnons, and that valley wave fluctuations limit the ordering temperature.","lead":"This paper computes the spin and valley magnon spectra of quantum anomalous Hall ferromagnets in twisted bilayer graphene using a microscopic Hartree-Fock plus Bethe-Salpeter approach. It finds the ferromagnet is stable, estimates the skyrmion pair energy, and argues that thermal valley fluctuations set the ordering temperature below the mean-field value, which may explain the small measured Curie temperature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The estimate TV≈ΔV/kB treats thermal excitation of a gapped valley wave as if it set the 2D Ising transition; it is not derived from the nonlinear sigma model and is the main soft spot.","rationale":"The reader's CONDITIONAL verdict identifies the T_V estimate as the weak point, and I agree that this is where the manuscript is least secure. My concern is slightly sharper: the issue is not only the unverified λ>a_M assumption, but the conceptual step that identifies the ordering temperature with the valley-wave gap. Thermal occupation of a gapped mode is a fluctuation correction to the order parameter, not the definition of the critical temperature; in the 2D Ising universality class the transition is governed by the balance of domain-wall free energy, which the paper itself estimates through T_DW. The paper's own Eq. (9) implies T_DW>2.62ΔV/k_B, so taking T_V=ΔV/k_B requires a model for how nonlinear interactions among valley waves renormalize the stiffness and lower the transition, which is not provided. The Monte Carlo test of the effective NLσM would settle whether 11 K is a real prediction of the model or an order-of-magnitude estimate. The stability of the QAHF against spin and valley particle-hole excitations is better supported: the Bethe-Salpeter spectra are nonnegative, the spin-wave Goldstone mode is correctly identified, and the SM provides a second parameter set (Δb=Δt=30 meV) where the single-band projection is cleaner and the stability conclusions survive. The single free parameter ε is fit to the experimental gap, which makes the quantitative 11 K number conditional, but that was already reflected in the reader's verdict. I am therefore not moving the verdict; 'no change needed' is the honest recommendation. If the proposed Monte Carlo test shows a materially different T_c, the T_V prediction would need to be downgraded to a qualitative statement.","tokens_in":13489,"tokens_out":6729,"duration_ms":68314,"concrete_test":"Run a classical Monte Carlo simulation of the lattice discretization of the anisotropic nonlinear sigma model in Eq. (8), using the microscopic values of ΔV and ρ⊥ extracted from the valley-wave spectrum at ε=30 and a ρz estimated from the microscopic model. Locate the Ising transition temperature from the Binder cumulant or from the peak of the specific heat and compare it with ΔV/k_B. If T_c differs from 11 K by more than about 30%, the paper's T_V estimate is not supported by its own model. A cheaper partial test is to determine λ microscopically; if λ/a_M is not significantly larger than 1, the ordering is controlled by domain walls and T_V≈ΔV/k_B fails on the paper's own criterion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative finite-temperature claim, T_V ≈ ΔV/k_B, is not actually derived. In the 'Valley wave' section the authors first estimate the domain-wall proliferation temperature from Eq. (9), k_B T_DW ≈ 2.62(λ/a_M)^2 ΔV, and argue λ>a_M so T_DW>2.62ΔV. They then state that valley waves are 'already thermally excited when k_B T exceeds ΔV' and conclude that T_V is limited by valley waves. This conflates having a thermally populated gapped mode with the loss of Ising order. In a 2D system with a discrete symmetry, the phase transition temperature is set by the free energy of topological defects and finite-temperature fluctuations, not by the single-particle gap; long-range order can survive well above a single-magnon gap. Indeed, if their own T_DW estimate applies, it gives the ordering scale and exceeds ΔV/k_B by at least a factor of 2.62 for λ>a_M. Thus the specific prediction T_V≈11 K rests on an unjustified proportionality, and the paper's own caveat that the agreement 'might be a coincidence' is apt. This is the load-bearing weakness because the paper's new contribution relative to concurrent work is precisely the valley-wave-limited T_V; the spin and valley stability part, supported by nonnegative Bethe-Salpeter spectra, is not endangered by this concern.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a microscopic theory of collective excitations of the quantum anomalous Hall ferromagnet (QAHF) at integer filling ν=3 in twisted bilayer graphene aligned with hBN. Using a continuum moiré model and Hartree-Fock mean-field states, the authors solve Bethe-Salpeter equations for spin-flip and valley-flip particle-hole excitations. They find nonnegative spin and valley excitation spectra, an exact gapless Goldstone spin wave at Q=0, and a gapped valley wave. From the spin-wave dispersion they extract the spin stiffness and estimate the skyrmion-antiskyrmion pair energy, which is comparable to the Hartree-Fock gap. From the valley-wave gap they estimate a valley ordering temperature T_V that is reduced from the mean-field transition temperature and is reported to be consistent with the experimental Curie temperature.","tokens_in":13654,"tokens_out":4360,"duration_ms":45668,"significance":"If the results hold, the paper provides a useful microscopic stability analysis of the QAHF state in TBG and a concrete mechanism (valley-wave fluctuations) for suppressing the ordering temperature below the mean-field value. The central stability claim is supported by the exact zero-energy spin wave at Q=0 required by SU(2) symmetry and by the nonnegative numerical spectra over the explored parameter range; this part is technically solid and will be of value to the moiré materials community. The quantitative finite-temperature prediction, however, rests on an ad hoc identification and needs to be underpinned by a proper finite-temperature calculation.","major_comments":[{"comment":"","section":"Valley wave, after Eq. (9)"},{"comment":"The quantitative agreement between the computed T_V (about 11 K) and the experimental Curie temperature is weakened by the fact that the dielectric constant ε is fitted to the experimental ν=3 charge gap of about 2 meV (and hence ΔV is not an independent prediction). This is not circular in the strict sense, but it reduces the independence of the finite-temperature comparison. Please state explicitly in the text that the T_V estimate is not a parameter-free prediction and, if possible, show the sensitivity of T_V to ε over a wider range to indicate how robust the 11 K value is.","section":"Ferromagnetism, fitting ε"}],"minor_comments":[{"comment":"The inset showing σ_{yx} as a function of temperature is informative, but the jump at T_MF would be easier to see if the figure included a vertical dashed line at the transition temperature in the inset as well as in the main panel.","section":"Fig. 2(b)"},{"comment":"The notation ρ_z and ρ_⊥ for the anisotropic valley stiffness is clear, but the text never defines n_0 explicitly in the main body; please state that n_0 is the density of one electron per moiré unit cell (it is defined in the Introduction in a slightly different context).","section":"Valley wave, Eq. (8)"},{"comment":"The sentence 'we addressed valley ordering temperature limited by valley wave excitations, which has not been studied previously in TBG' would be more precise as 'we address ...', since the paper is being submitted for publication rather than retrospectively announcing completed work.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":"The microscopic stability analysis is sound and publishable on its own. The main weakness is the derivation of T_V ≈ ΔV/k_B, which is not justified and is the paper's claimed new contribution. I would encourage the editor to require either a proper finite-temperature derivation from the NLSM or a clear reframing of T_V as a heuristic estimate before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read the Wu–Das Sarma paper on collective excitations of the TBG quantum anomalous Hall ferromagnet. The core new thing here is the calculation of both spin and valley magnon spectra from Bethe–Salpeter equations, plus the attempt to use the valley wave gap to estimate the ordering temperature. The stability part of the paper is well done: the spin wave is gapless at Q=0 as required by Goldstone's theorem, the spectra are nonnegative over the parameter range explored, and the identification of the gapped valley wave as the reason the valley-polarized state beats the valley-coherent state is sensible. The skyrmion-pair energy being comparable to the Hartree–Fock gap is a useful observation, and the authors are properly cautious about which one is the lowest charged excitation.\n\nThe paper also earns credit for being transparent. They openly note the three concurrent preprints and state that the valley-wave-limited ordering temperature analysis is new relative to those. They also flag, without being asked, that the quantitative agreement with the experimental Curie temperature \"might be a coincidence.\" That honesty matters.\n\nThe main soft spot is the finite-temperature claim, T_V ≈ ΔV/k_B. As far as I can see, this is not derived. The argument is that valley waves are already thermally excited when k_B T exceeds ΔV, so they should limit the ordering. But thermal occupation of a gapped mode does not by itself destroy Ising order; in two dimensions the transition is set by topological defect proliferation, and the single-particle gap can be well below the ordering scale. The paper's own domain-wall estimate gives T_DW > 2.62 ΔV/k_B, which would be a higher temperature, so the claim that valley waves dominate relies on an inequality that is plausible but not proven. This does not sink the stability analysis, but it means the specific prediction T_V ≈ 11 K should be read as a heuristic estimate, not a computed consequence.\n\nThere is a second, lesser concern: the single fitting parameter ε is chosen to match the experimental charge gap, and then the computed T_V is compared with the Curie temperature. That is not circular, since different observables are involved, but it does reduce the independence of the test. Also, at the fitted ε the interaction scale is comparable to the charge neutrality gap, so the single-band projection is only marginally justified. The authors address this in the supplemental material with the Δ_b = Δ_t case, which is good, but the main-text numbers sit at the edge of the model's validity.\n\nOverall, the central claim — that the ν=3 QAHF is robust against spin and valley particle-hole fluctuations — is likely correct and is supported by the computed spectra. The valley ordering temperature is the one place where the theory is under-derived. This paper is for people working on moiré magnets and topological flat bands; it deserves a serious referee and, if the heuristic nature of T_V is made clearer, publication. I would cite it for the magnon spectra and stability analysis.","headline":"Solid microscopic theory for spin and valley magnons in TBG QAHF; the stability analysis is robust, but the valley-ordering temperature estimate is a heuristic argument rather than a derivation.","tokens_in":14290,"tokens_out":1950,"would_cite":true,"duration_ms":21738,"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":"The ν=3 quantum anomalous Hall ferromagnet in twisted bilayer graphene is stable against spin and valley magnons, and its valley ordering temperature is set by the gapped valley-wave mode.","keywords":["twisted bilayer graphene","quantum anomalous Hall ferromagnet","magnons","Bethe-Salpeter equation","valley Ising order","skyrmions","nonlinear sigma model","moiré materials"],"falsifier":"Measure the ν=3 anomalous Hall resistance as a function of temperature in devices with several dielectric environments and compare the ordering temperature with the zero-temperature valley-wave gap: if $T_c$ tracks the mean-field gap rather than $\\Delta_V/k_B$, the valley-wave-limited-ordering claim fails. Alternatively, a calculation including remote moiré bands that finds a negative valley-wave energy at any momentum $Q$ would disprove the stability conclusion.","tokens_in":13152,"feed_emoji":"🧲","tokens_out":13886,"duration_ms":129221,"temperature":0.7,"pith_summary":"This paper asks whether the quantum anomalous Hall ferromagnet at filling ν=3 in twisted bilayer graphene—an ordered state with quantized Hall resistance at zero magnetic field—is stable, and what controls its ferromagnetic transition temperature. Solving Bethe-Salpeter equations for the two collective excitations, waves with one spin flipped and waves with one valley flipped, the authors find both spectra nonnegative, so the state survives small particle-hole fluctuations. The spin wave is gapless, as Goldstone's theorem requires, and its stiffness implies that a skyrmion-antiskyrmion pair costs energy comparable to the Hartree-Fock gap $\\Delta_{HF}$. The valley wave is gapped, favoring valley-polarized over valley-coherent order, and this gap limits the ordering temperature to roughly $k_B T_V\\approx\\Delta_V$, far below the mean-field value. If correct, this explains why the measured transport gap of the ν=3 state is several meV while its Curie temperature is only about 9 K.","feed_headline":"11 K: valley waves set twisted graphene's Curie temperature","feed_subtitle":"Thermal valley magnons push the transition well below mean field and explain the gap-versus-temperature mismatch.","key_machinery":"The load-bearing machinery is a ladder sum of repeated electron-hole interactions: a magnon is written as a particle-hole pair with variational coefficients $z_{k,Q}$, and minimizing its energy yields the Bethe-Salpeter equations $E(Q)z_{k,Q}=\\sum_{k'}K^{(Q)}_{kk'}z_{k',Q}$, whose kernel is the quasiparticle energy cost of the flip minus the attractive electron-hole interaction. A closed-loop gauge-invariance argument makes the eigenvalues well defined despite the phase ambiguity of Bloch wave functions. The long-wavelength results are captured by two nonlinear $\\sigma$ models: an O(3) model for the spin with stiffness $\\rho_s$, and a valley model with Ising anisotropy $u$ and anisotropic stiffnesses $\\rho_z,\\rho_\\perp$. The domain-wall ansatz $(\\pi_x,\\pi_y,\\pi_z)=(\\operatorname{sech}(x/\\lambda),0,\\tanh(x/\\lambda))$ supplies the domain-wall energy $J=4u\\lambda$ used to compare valley-wave and domain-wall ordering temperatures.","core_discovery":"The central claim is that the ν=3 quantum anomalous Hall ferromagnet in hBN-aligned twisted bilayer graphene is a stable valley Ising ferromagnet whose two collective branches determine the transport and ordering properties. The spin magnon spectrum, from the Bethe-Salpeter equation $E_S(Q)z_{k,Q}=\\sum_{k'}H^{(Q)}_{kk'}z_{k',Q}$, has a gapless Goldstone mode at $Q=0$ with $E_{SW}=(2\\rho_s/n_0)Q^2$; its nonnegative spectrum verifies stability against spin fluctuations. The valley magnon spectrum from the analogous equation has a gapped lowest mode $E_{VW}=\\Delta_V+(2\\rho_\\perp/n_0)Q^2$, showing that valley-polarized order is preferred over valley-coherent order. From the spin stiffness $\\rho_s$ the paper obtains a skyrmion-antiskyrmion pair energy $\\Delta_{pair}=8\\pi\\rho_s$ comparable to the Hartree-Fock gap $\\Delta_{HF}$, and capable of being the lowest charged excitation. From the valley-wave gap $\\Delta_V$, a nonlinear $\\sigma$ model with Ising anisotropy yields $k_B T_V\\approx\\Delta_V$, reduced from the mean-field $T_{MF}$, and gives about 11 K for a 2 meV gap, close to the experimental Curie temperature.","pith_inferences":["A testable consequence of the paper's mechanism is that the ordering temperature should scale with the zero-temperature valley-wave gap as screening is changed; a device with a smaller $\\Delta_V$ should show a proportionally lower $T_c$.","In the parameter regime where $\\Delta_{pair}<\\Delta_{HF}$, the activated transport at low temperature may be carried by skyrmion-antiskyrmion pairs, which would give an activation energy and Hall signature distinct from single-particle excitations.","The gapped valley-wave result suggests that intervalley coherent states in other flat moiré bands may be generically disfavored whenever the valley wave is gapped, not only at ν=3.","Coupled spin-valley fluctuations, which the paper only sketches, could be what turns the transition first-order; a quantitative study of the coupled sigma models would test that scenario."],"forward_implications":["With a finite Hartree-Fock gap $\\Delta_{HF}$, the ν=3 state is stable against both spin-flip and valley-flip particle-hole excitations, so the observed quantized Hall effect is not on the verge of a collective instability.","The lowest charged excitation may be a skyrmion-antiskyrmion pair rather than a single-particle excitation when both hBN layers are aligned, changing how the measured transport gap should be interpreted.","The ferromagnetic ordering temperature is limited by thermal valley waves, so the transport gap and the Curie temperature need not be proportional; a 2 meV gap gives $T_V\\approx 11$ K, near the observed 9 K.","The model predicts a valley-polarized (valley Ising) ground state rather than a valley-coherent one, because the valley wave is gapped.","The same Bethe-Salpeter approach can be carried over to other broken-symmetry states in flat moiré bands."],"supporting_citations":[{"why":"Reports the experimental ν=3 quantum anomalous Hall effect in hBN-aligned twisted bilayer graphene, including the ~2 meV transport gap, ~9 K Curie temperature, and resistance jump that the paper's estimates target.","marker":"[8]"},{"why":"Argues for the valley-polarized versus valley-coherent energetics that the gapped valley-wave mode supports.","marker":"[53]"},{"why":"Provides the quantum Hall ferromagnet skyrmion formalism whose pair energy $\\Delta_{pair}=8\\pi\\rho_s$ and comparison with $\\Delta_{HF}$ is adopted here.","marker":"[64]"},{"why":"Supplies the continuum moiré Hamiltonian used to compute the flat bands and Berry curvature that seed the Hartree-Fock state.","marker":"[66]"},{"why":"Provides the experimental basis for the sublattice-potential parameters $(\\Delta_b,\\Delta_t)=(30,0)$ meV used in the main calculation.","marker":"[68]"},{"why":"Underwrites the gauge-invariance argument that makes the Bethe-Salpeter magnon eigenvalues well-defined.","marker":"[69]"},{"why":"Supplies the O(3) nonlinear sigma model from which the spin stiffness and skyrmion energy are extracted.","marker":"[70]"},{"why":"Gives the domain-wall proliferation estimate $k_B T_{DW}=2.62(\\lambda/a_M)^2\\Delta_V$ used to decide that valley waves, not domain walls, limit the ordering temperature.","marker":"[72]"}],"fun_headline_variants":["Valley waves set twisted graphene's magnetic order at 11 K","Skyrmion pair energy rivals gap in quantum anomalous Hall ferromagnet","Gapped valley magnons drive low Curie temperature in TBG","Twisted bilayer graphene: valley waves control ordering at 11 K"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The estimate $k_B T_V\\approx\\Delta_V$ assumes thermally excited valley waves destroy the valley Ising order before domain-wall proliferation does, which requires the domain-wall width $\\lambda$ to exceed the moiré period $a_M$; if instead $\\lambda\\le a_M$, domain walls could set the transition temperature.","fun_headline_variants_meta":{"raw":{"variants":["Valley waves set twisted graphene's magnetic order at 11 K","Skyrmion pair energy rivals gap in quantum anomalous Hall ferromagnet","Gapped valley magnons drive low Curie temperature in TBG","Twisted bilayer graphene: valley waves control ordering at 11 K"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1643,"prompt_tokens":966,"completion_tokens":677,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":601}},"tokens_in":582,"tokens_out":677,"duration_ms":6348,"temperature":1.0,"reasoning_tokens":601,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:14:12.500350+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ν=3 anomalous Hall resistance as a function of temperature in devices with several dielectric environments and compare the ordering temperature with the zero-temperature valley-wave gap: if $T_c$ tracks the mean-field gap rather than $\\Delta_V/k_B$, the valley-wave-limited-ordering claim fails. Alternatively, a calculation including remote moiré bands that finds a negative valley-wave energy at any momentum $Q$ would disprove the stability conclusion.","supporting_citations":[{"cited_title":"Anomalous Hall effect, magneto-optical properties, and nonlinear optical properties of twisted graphene systems","cited_arxiv_id":"1907.08932","evidence_quote":"Provides the quantum Hall ferromagnet skyrmion formalism whose pair energy $\\Delta_{pair}=8\\pi\\rho_s$ and comparison with $\\Delta_{HF}$ is adopted here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the continuum moiré Hamiltonian used to compute the flat bands and Berry curvature that seed the Hartree-Fock state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Underwrites the gauge-invariance argument that makes the Bethe-Salpeter magnon eigenvalues well-defined."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the domain-wall proliferation estimate $k_B T_{DW}=2.62(\\lambda/a_M)^2\\Delta_V$ used to decide that valley waves, not domain walls, limit the ordering temperature."}],"review_version":1}