{"id":"3ae0078b-7db5-415b-9acd-cf9ef822d697","arxiv_id":"2601.06778","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In doped twisted bilayer graphene, moiré phonons—particularly the phason—become infrared active, with the phason carrying exactly one electron per moiré cell of doping and giving a Drude-like THz response.","lead":"Adding electrons to twisted bilayer graphene makes its slow moiré vibrations—especially the sliding phason mode—carry charge and absorb light, creating sharp THz resonances inside the electronic gap. The paper derives these resonances from a microscopic model and predicts the phason's effective mass and topological charge, offering a measurable probe of lattice relaxation and electron-phonon coupling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gap persistence along the sliding path is unverified; a closure would break the quantized phason charge and Drude formula.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the quantization of the sliding Chern number and the resulting Drude formula require an adiabatic gap for every intermediate stacking. The paper gives a detailed derivation, analytical proofs of phason gaplessness and the f-sum rule, and numerical convergence checks, but it does not verify the gap along the sliding path. This is an addressable numerical check rather than a fatal flaw, so the verdict remains CONDITIONAL (UNCHANGED).","tokens_in":36672,"tokens_out":16115,"duration_ms":178946,"concrete_test":"For theta=1.05 deg (and 1.54 deg, 2 deg), diagonalize H[phi0(r-zeta)] on a grid of zeta spanning one moire cell, using the same 10-star basis and full u0(r). At each zeta, compute the direct gap between occupied and empty bands. If the gap remains open, evaluate C(zeta)=2 sum_occ int_mBZ d^2k/(2pi)^2 Omega_{kx,zeta_x}(k) and verify it is zeta-independent and equals -n; if the gap closes at any zeta, repeat with larger truncations to determine whether the closing is a finite-basis artifact. This directly tests the quantization underpinning Eqs. 11-15.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the phason carries a quantized charge C=-n (Eq. 12) and yields the Drude-like conductivity (Eq. 14) depends on the electron system remaining gapped at every intermediate stacking configuration zeta along the adiabatic sliding path. Eq. 12 is only a well-defined integer when the occupied manifold is separated by a gap for all zeta; otherwise the mixed Berry curvature integral is not quantized and C can take non-integer values. The paper asserts in the Discussion that the argument applies \"as long as the electron bands remain gapped,\" but no numerical check of the gap along the path is shown for any twist angle. The convergence data in Fig. 4 themselves show that reaching C=-n requires including many bands and the full relaxation field u0(r) explicitly in the electronic Hamiltonian (Supplement S2.B), so the quantization is not robust to common approximations. If the gap closes during sliding, the effective mass m* (Eq. 15) and the Drude weight in Eq. 14 would change accordingly. This is the load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that moiré phonons in twisted bilayer graphene, in particular the acoustic phason, acquire an electric dipole moment when the flat bands are doped away from charge neutrality. The mechanism is the interband electron-phonon coupling that transfers spectral weight from the electron-hole continuum to the phonon resonances. The central formal result is that the phason carries a topologically quantized charge equal to the doping n per moiré cell, expressed through a sliding Chern number C=-n, and that this charged phason gives a Drude-like low-frequency optical conductivity with an effective mass that increases with lattice relaxation. The numerical part computes the optical absorption at three twist angles, including the magic angle, with and without relaxation, and checks convergence up to 60 bands.","tokens_in":36947,"tokens_out":11394,"duration_ms":109108,"significance":"If the central claim holds, the paper identifies a generic and experimentally accessible phenomenon: infrared-active moiré phonons and a quantized sliding phason charge in doped twisted bilayer graphene. The THz response predicted here would provide a direct probe of electron-phonon coupling and disorder in moiré systems. The manuscript has notable strengths: the supplementary information contains explicit derivations of the electron-phonon matrix elements, a proof that the phason remains gapless in the model, a formal proof that the computed diagrams exhaust the optical f-sum rule, and numerical convergence tests up to 60 bands. The numerical confirmation of the analytical Drude formula in Fig. 4, when it works, is a non-trivial check. The main risk is that the topological quantization C=-n and the resulting Drude formula require the electron system to remain gapped along the entire adiabatic sliding path, and this condition is asserted but not verified.","major_comments":[{"comment":"The quantization C=-n and the Drude-like formula Eq. (14) presuppose that the occupied manifold remains separated by a gap for every intermediate stacking configuration along the sliding path. The manuscript only states this condition in the Discussion (\"as long as the electron bands remain gapped\") and gives no calculation of the miniband gap as a function of the sliding coordinate ζ for any twist angle. In addition, the claim that the phason charge equals the number of added/removed electrons is made for arbitrary doping, but at partial fillings n=±1,±2,±3 the noninteracting band structure is metallic; an interaction-induced gap would be needed, and such a many-body gap is not analyzed. This is a load-bearing assumption: if the gap closes during sliding, C in Eq. (12) is not an integer and Eq. (14) does not follow. Please provide a gap-vs-ζ calculation for the angles and fillings used,","section":"Charged phason, Eq. (12); Discussion"},{"comment":"The prefactor in Eq. (11) contains \\hbar^2 e^2, but the low-frequency derivation in SI S3.B, Eq. (S47b), yields a phason conductivity without any factor of \\hbar^2. With physical units restored, Eq. (11) is dimensionally inconsistent and does not reduce to the Drude form Eq. (14) with m* given by Eq. (15) unless \\hbar=1 is assumed without being stated. Relatedly, with n dimensionless (electrons per moiré cell), m* = \\theta^2\\varrho/(2\\eta^2 n) has units of mass per area rather than mass; the Drude formula then uses the dimensionless filling n, but the text calls m* an \"effective mass\" and does not specify which density convention is used. Please correct the prefactor and state the units of m* explicitly.","section":"Eq. (11) and SI S3.B"},{"comment":"The supplementary information itself states that reproducing the electron-phason coupling requires the relaxation field u0(r) to appear explicitly in the electronic Hamiltonian, not just as renormalized model parameters, and that this is decisive for capturing the quantized response. This is an internal limitation of the central numerical claim: the quantization is not robust to the common approximation of absorbing relaxation into model parameters. The authors should state this limitation in the main text. In addition, Fig. 4 shows convergence of the conductivity to the analytical curve, but not the value of C itself as a function of band truncation; reporting C directly for increasing band number would make the quantization claim more transparent.","section":"SI S2.B and Fig. 4"}],"minor_comments":[{"comment":"The RPA resummation is justified by \"the large (N=4) number of fermion species\"; four is not large, and the statement should be softened or supplemented by a more quantitative justification.","section":"Model, Eq. (8)"},{"comment":"The caption says the phason peaks are \"pinned to 0.01ω_m\" for visualization, while Fig. 4 is described as \"with no pinning\". This is understandable, but the distinction between a numerical pinning frequency and physical pinning by disorder should be clarified in the text.","section":"Fig. 2 caption"},{"comment":"The formula for η contains a typographical artifact \"1q\" and should read 1/sqrt(...). Please check the final typeset version.","section":"Eq. (10)"},{"comment":"The text refers to red and blue vertices in Fig. 3, but the printed figure may not show colors; please use shape-based labels or otherwise make the reference unambiguous.","section":"Optical conductivity, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is timely and the SI is detailed, but the central quantization claim rests on an unverified gap-persistence condition along the sliding path. This is the main risk. The spurious \\hbar^2 in Eq. (11) appears to be a fixable typo, but it is currently load-bearing for the relation between the central equations and should be corrected. The numerical convergence checks are a genuine strength, and with the requested gap calculations and unit clarifications the paper could become a strong contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth reading. It computes the optical conductivity of twisted bilayer graphene with E1 moiré phonon vertices and finds infrared-active phonons inside the single-electron gap when the flat band is filled. The new piece is the direct microscopic calculation of the charged-phonon response and the phason Drude formula with an effective mass that grows with relaxation. Individual ingredients — charged phonons, moiré phonons, sliding charge pumping — were known; the combination and the specific predictions are not. The SI is thorough: explicit EPC matrix elements, an f-sum rule proof, a phason stability proof via the diamagnetic/paramagnetic cancellation, and convergence checks up to 60 bands. That is real work and it holds together.\n\nThe soft spots are proportionate. The phonon damping is fitted from transport with theta^-3 scaling — fine for an estimate, not a prediction. The deformation potentials and tunneling amplitudes are model inputs. More substantively, the quantized phason charge C=-n and the Drude formula assume the electron system stays gapped along the adiabatic sliding path. The paper states this in the Discussion but never checks it numerically for all intermediate stackings. The authors do show that reaching C=-n requires many bands and the full relaxation field in the Hamiltonian, so the quantization is not an artifact of a crude model; but the gap-persistence premise is a real unverified step. It is a condition for the central claim, not a contradiction of it.\n\nWho is this for? Experimentalists doing THz spectroscopy on tBG and theorists working on moiré mechanics. It is a solid within-subfield advance. It deserves a serious referee. I would send it to review, asking specifically for a gap check along the sliding path and a brief comment on how sensitive the Drude weight is to the fitted gamma.","headline":"A direct calculation of charged moiré phonon optics in tBG with real substance; the main open point is the assumed gap persistence along the sliding path.","tokens_in":37416,"tokens_out":1496,"would_cite":true,"duration_ms":17514,"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":"Doped twisted bilayer graphene makes its phason — a sliding moiré phonon — carry a quantized charge equal to the doping, and absorb light in the terahertz range.","keywords":["twisted bilayer graphene","moiré phonons","phasons","sliding Chern number","optical conductivity","electron-phonon coupling","terahertz spectroscopy","lattice relaxation"],"falsifier":"Measure the low-frequency optical conductivity of magic-angle twisted bilayer graphene at full filling of the flat band. If there is no narrow Drude-like resonance within the single-electron gap whose weight scales with the square of the doping and whose mass grows with relaxation (e.g., under pressure), the central claim is falsified.","tokens_in":36575,"feed_emoji":"⚡","tokens_out":3446,"duration_ms":34321,"temperature":0.7,"pith_summary":"The paper tries to show that in twisted bilayer graphene, the ultra-low-energy collective vibrations of the moiré pattern — in particular the sliding mode called the phason — acquire an electric dipole moment as soon as electrons are added or removed from the flat bands. Because the charge redistribution follows the stacking pattern, a sliding of the honeycomb layers acts like an adiabatic pump that drags electrons with it. The paper argues that the charge of the phason is topologically quantized and equals the doping level measured from neutrality, and that this charged phason produces a measurable low-frequency (terahertz) optical absorption inside the single-electron gap, even though the system is a band insulator at full filling. If true, this turns optical absorption into a direct probe of the electron–phonon coupling and of disorder in moiré materials.","feed_headline":"Sliding phonons in doped twisted graphene carry one electron of charge","feed_subtitle":"Theory: each added electron per moiré cell charges the phason and creates a THz Drude peak inside the gap.","key_machinery":"The identifying object is the sliding Chern number C: the integral of the sliding Berry curvature over the moiré Brillouin zone, summed over occupied bands. It is quantized when the electron system is gapped, equals minus the doping, and enters the phason conductivity formula. The phason normal mode is represented by a collective coordinate that translates the relaxed stacking configuration, with a geometric factor that makes the effective mass depend on lattice relaxation.","core_discovery":"The paper claims that moiré phonons of E1 symmetry, although back-folded acoustic modes of the layers, become infrared active when the flat bands are filled or emptied. The phason in particular acquires an electric charge equal to the number of electrons per moiré cell added or removed relative to neutrality. This charge is a topological invariant — a sliding Chern number built from a mixed Berry curvature that pairs the electronic crystal momentum with the coordinate that slides the relaxed moiré pattern. As a consequence, the low-frequency optical conductivity acquires a Drude-like term with an effective mass that increases with lattice relaxation. The numerical calculation shows convergen","pith_inferences":["Because the phason carries charge, an electric field should exert a shear force between the layers, suggesting an electromechanical way to drive or detect sliding of the moiré lattice in a device — a consequence the paper hints at but does not develop.","The same sliding-Chern-number argument would apply to any moiré homobilayer with gapped bands, so charged phasons may be generic in twisted van der Waals stacks, not just in twisted bilayer graphene.","The quantized charge implies a topological contribution to the dielectric or second-order nonlinear optical response near the phason frequency, which could be tested by harmonic generation experiments.","Pinning by contacts or disorder converts the Drude peak into a finite-frequency resonance; measuring its position as a function of sample dimensions could separate contact pinning from bulk disorder."],"forward_implications":["At full flat-band filling, the system should display a DC-like conductivity from sliding even though it is a band insulator; the phason resonance is the optical signature.","The phason peak position and width should be controlled by disorder (pinning) and by pressure, which tunes layer adhesion and thus the effective mass.","The oscillator strength of the moiré phonon modes is a direct measure of electron-phonon coupling strength in moiré materials.","The effect is not restricted to the magic angle: peak frequencies scale with twist angle, so THz spectroscopy across twist angles can test the mechanism.","The phason remains gapless in the clean limit and the optical sum rule is preserved, so the charged-phonon spectral weight is borrowed from the electron-hole continuum."],"fun_headline_variants":["Moiré phasons acquire electron charge in doped twisted bilayer graphene","Sliding phonons in tBG gain charge per moiré cell when doped","Charged phasons in twisted bilayer graphene give Drude-like THz peak","Doped twisted graphene: phason charge equals added electrons per moiré cell","Infrared-active phasons in tBG carry one electron charge"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantization of the phason charge and the resulting Drude formula assume the electron system remains gapped for every intermediate stacking configuration along the adiabatic sliding path, and that the full relaxation field is included explicitly in the electronic Hamiltonian; if the gap closes or relaxation is only treated perturbatively, the equality C=-n and the effective-mass formula break down.","fun_headline_variants_meta":{"raw":{"variants":["Moiré phasons acquire electron charge in doped twisted bilayer graphene","Sliding phonons in tBG gain charge per moiré cell when doped","Charged phasons in twisted bilayer graphene give Drude-like THz peak","Doped twisted graphene: phason charge equals added electrons per moiré cell","Infrared-active phasons in tBG carry one electron charge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1272,"prompt_tokens":741,"completion_tokens":531,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":485,"tokens_out":531,"duration_ms":6399,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:17:26.586715+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the low-frequency optical conductivity of magic-angle twisted bilayer graphene at full filling of the flat band. If there is no narrow Drude-like resonance within the single-electron gap whose weight scales with the square of the doping and whose mass grows with relaxation (e.g., under pressure), the central claim is falsified.","supporting_citations":[],"review_version":1}