{"id":"63b15bf5-4e2e-48c0-8e21-e4affb059109","arxiv_id":"2608.09932","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Cell natural orbitals, defined from the unit-cell reduced density matrix, give an optimally truncated local-orbital representation of topological bands and a hierarchy of projected interaction strengths, demonstrated in chiral twisted bilayer graphene.","lead":"The paper introduces a decomposition of band wavefunctions into 'cell natural orbitals' ranked by occupation, and applies it to twisted bilayer graphene to show that a few local orbitals capture the dominant interactions. The framework offers a systematic way to truncate interacting topological band models, but rests on an unpublished companion paper and several unproven shortcuts.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central interaction-hierarchy claim depends on an unproven transfer from the periodic-gauge CNO kernel Lambda (Eq. 12) to the physical, non-periodic form factor Lambda_phys (Eq.","rationale":"The reader's verdict is CONDITIONAL and the reader's weakest assumption is the transfer from Lambda to Lambda_phys; my analysis agrees. The paper's rigorous components are real: the Eckart-Young argument for Lambda is correct, the topological-zero theorem for a fixed CNO is sound, and the chiral TBG numerics (exponential decay after the second CNO, vanishing Im F_11(q), AA-site localization) are a genuine demonstration. The weakness is that the central practical claim, that CNO occupations rank interaction channels and determine the minimal local-orbital set, is not derived from the proved statements. The physical form factor Lambda_phys is a different object, and the CNO decomposition is explicitly not its SVD; the hierarchy in interactions is therefore an empirical observation from one model. This is not internal inconsistency, so the paper is not rejected; it is a missing generalization that a second model test or an analytic bound would settle. I also note the separated-zero two-CNO argument is an existence claim that does not follow for the CNO basis itself. These observations do not change the reader's verdict: CONDITIONAL remains appropriate.","tokens_in":24489,"tokens_out":11362,"duration_ms":104342,"concrete_test":"Take a non-TBG Chern lattice model with C=1 (Haldane) and C=2 (e.g., the spin-S multifold lattice model of Eq. 39 with S=1). Compute the CNOs from L and then evaluate the exact physical form-factor matrix M(q)_{k,k'} = Lambda_phys_{k,k+q} with the same screened Coulomb V_q used in the paper. Compare (i) the relative Hilbert-Schmidt error ||M(q) - M_CNO^{(N_tau)}(q)|| / ||M(q)|| for N_tau=1, 2, 3, and (ii) the on-site U_a from Eq. 30 sorted by lambda_a. If U_a is not monotonically ordered by lambda_a, or if the rank-2 reconstruction error of M(q) is not small (say >10%) while the Lambda-truncation error sum_{a>2} lambda_a^2 is small, then the interaction hierarchy seen in chiral TBG is model-specific and the central claim fails as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest claim is that the CNO occupation spectrum sets a hierarchy in projected interactions. The rigorous parts are the SVD and truncation-error statement for the overlap kernel Lambda (Eqs. 12-15 and Appendix E) and the topological-zero theorem for s_{1,k}. Neither of these, however, controls the physical form factor Lambda_phys entering the projected density operator (Eq. 24). The authors explicitly state in Sec. III A that Lambda is not the exact SVD of the physical form-factor matrix, and in Sec. IV B they concede that the origin of this hierarchy is not obvious a priori and that it works well for magic-angle twisted bilayer graphene in the chiral limit. The on-site interaction strengths U_a (Eq. 30) and the exponential locality after N_tau=2 are computed from F_ab(q) in Eq. 29, which mixes Lambda_phys with CNO envelope functions; nothing in the Eckart-Young optimality for Lambda forces this combination to be ranked by lambda_a. The claim therefore rests on a numerical observation in a single, highly special model (chiral TBG, with AA-site charge concentration and chiral symmetry). The paper also asserts without proof that two CNOs suffice for arbitrary Chern number when zeros are separated, via a well-conditioned transformation to a single-zero representation; because CNOs are fixed by L, such a transformation would change the orbitals, so the original CNO basis is not shown to supply the needed second mode. These gaps make the general framework conditional rather than established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces Cell Natural Orbitals (CNOs) as eigenstates of the unit-cell reduced one-particle density matrix built from the band projector, and develops a singular-value decomposition of the band-projected overlap kernel. The formal results include: (i) the Eckart–Young optimality of the CNO truncation with Hilbert–Schmidt error equal to the sum of discarded squared occupations (Sec. III C and Appendix E); (ii) a topological obstruction theorem stating that the leading CNO envelope must vanish on a Chern band (Sec. III B); and (iii) a claimed hierarchy in band-projected interactions, with numerical evidence in chiral twisted bilayer graphene at the magic angle (Sec. IV B). The paper also studies unit-cell entanglement, complexity measures, and momentum-dependent self-energies in a Hubbard-I approximation. The central physical claim is that the CNO occupation spectrum sets a systematically improvable local-orbital hierarchy for interacting topological bands.","tokens_in":24788,"tokens_out":3436,"duration_ms":32929,"significance":"If the central claim holds, CNOs provide a valuable, parameter-free construction of local orbitals for topological flat bands, with a natural ranking of interaction channels. The formal parts of the paper are genuinely rigorous: the truncation-error proof (Appendix E) is self-contained, the line-bundle argument for the protected zero of the leading envelope is concise and correct, and the Peschel relation is used appropriately. The chiral TBG numerics (Figs. 1, 5, 6) are concrete and will be useful to the community, and the explicit comparison of power-law versus exponential interaction tails as a function of the Gram-matrix singularity (Sec. V, Fig. 9) is instructive. The main weakness is that the decisive transfer from the overlap-kernel SVD to the physical form factor is asserted rather than proved, and the general multi-zero argument for arbitrary Chern number is not established.","major_comments":[{"comment":"The central claim that the CNO occupation spectrum sets a hierarchy in projected interactions is not derived from the SVD optimality of the overlap kernel. The on-site form factor F_ab(q) in Eq. (29) contains the physical form factor Λ_phys and the CNO envelopes divided by sqrt(S_k S_{k+q}); Eckart–Young optimality for Λ does not constrain this combination. The paper explicitly concedes in Sec. III A that \"our decomposition is not the exact SVD of the physical form-factor matrix\" and in Sec. IV B that \"the origin of this hierarchy is not obvious a priori.\" The hierarchy observed in Fig. 6 is therefore a numerical observation in one model (chiral TBG), not a demonstrated general property. The third central result in Sec. VI should either be proved or restated as a conjecture with a precise condition under which it holds.","section":"Sec. IV B, Eq. (29)"},{"comment":"The assertion that \"since all bands with the same Chern number belong to the same topological class, there should exist a well-conditioned transformation mapping a band with multiple zeros to one with a single zero of higher vorticity\" is not proven, and it is not sufficient for the stated conclusion. Even if such a transformation exists, it acts on the band wavefunctions, not on the CNOs, which are fixed by the unit-cell reduced density matrix L. Consequently the original CNO basis is not shown to supply the needed second mode when zeros are separated. The proof only covers the coincident-zero case via completeness; the general claim that Nτ=2 regularizes S_k for arbitrary Chern number remains unsubstantiated.","section":"Sec. IV A, paragraph on separated zeros"},{"comment":"The hierarchy U_f > U_cf > U_c follows exactly in Appendix G only for the special all-to-all unit-cell-density interaction of Eq. (G7), where the interaction scales are squared CNO occupations. The paper then states in the same appendix that \"the hierarchy persists for realistic Coulomb interactions\" without a derivation. Since the realistic case is precisely the one used for the main physical claim in Sec. IV B, the relation between the exactly solvable limiting case and the generic case should be made explicit, and the unexplained step should be flagged as an assumption.","section":"Sec. VI vs Appendix G"}],"minor_comments":[{"comment":"There is a typo: \"We refer the read to Ref. [1]\" should be \"We refer the reader to Ref. [1].\"","section":"Sec. II, first paragraph"},{"comment":"The word \"wavefucntions\" in the sentence following Eq. (4) is misspelled; it should be \"wavefunctions.\"","section":"Sec. II, Eq. (4) text"},{"comment":"References [57] and [58] appear to be the same paper (B. Mera and T. Ozawa, Phys. Rev. B 106, 245134 (2022)). The duplicate should be removed and the citation numbers adjusted.","section":"References [57,58]"},{"comment":"The caption for Fig. 7 says \"with (c) its imaginary part\" but panels (a)–(d) are not individually labeled in the caption text; please align the caption references with the figure panel labels for readability.","section":"Sec. IV B, Fig. 7 caption"},{"comment":"The statement that equality in Eq. (13) holds only if the band is a single k-independent orbital is correct for a single band, but for a multi-band projector the equality condition should be stated for each band separately; a one-sentence clarification would avoid ambiguity.","section":"Sec. III B, Eq. (13)"},{"comment":"The term \"Mott Semimetal phase\" appears with inconsistent capitalization; please harmonize with \"Mott semimetal\" used elsewhere.","section":"Appendix H, last paragraph"}],"recommendation":"major_revision","confidential_remarks":"The paper depends heavily on a companion paper (Ref. [1], \"to appear\") for the foundational derivation of CNOs and the periodic embedding in continuum models. Since the present manuscript's central physical claim is explicitly conditional on a numerical observation in one model, I would recommend that the editor request the companion paper (or a preprint version) for the review process, and that the authors strengthen the generality of the interaction-hierarchy claim or narrow it to a clearly stated conjecture. The formal core is sound and the paper is likely to be influential, so the main revision effort should focus on the unproven transfer from the overlap kernel to the physical form factor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper. The formal core holds up. The SVD of the band-projected overlap kernel via the unit-cell RDM is correct, the Eckart-Young truncation error equals the sum of discarded squared occupations, and the protected zero of the leading envelope is a clean restatement of the known bundle obstruction. The Peschel relation and the cell-entropy bound are fine. Those parts are rigorous and give the community something genuinely useful: an ansatz-free, prescribed way to choose local orbitals and rank them by occupancy.\n\nThe new and useful bit is the demonstration that in chiral TBG, two CNOs regularize the Gram matrix, kill the power-law tails in the projected interaction, and give an on-site hierarchy with the AA-centered mode dominant. That is concrete and believable, and the figures show it. The Hubbard-I spectral function with the gapless crossing at Gamma is a nice illustration, and the authors are honest that it is a single-band toy and that the two-sector crossing is missed.\n\nSoft spots, in order of real weight. First, the conclusions' claim that \"the occupation spectrum sets a hierarchy in projected interactions\" is not proven. The rigorous statements control the kernel Lambda of Eq. (12), not the physical form factor Lambda_phys entering the projected density in Eq. (24). The authors state in Sec. III A that the decomposition is not the exact SVD of the physical form-factor matrix, and in Sec. IV B they concede that the origin of the hierarchy is not obvious a priori. So the hierarchy is a numerical observation in one highly special model, with AA-site charge concentration and chiral symmetry doing much of the work. That is acceptable as a conjecture, but the central-claims section should not present it as an established result of the framework. Second, the two-CNO-sufficiency claim for arbitrary Chern number is hand-waved: the well-conditioned transformation argument would change the orbitals, while CNOs are fixed by L, so it does not show the original CNO set supplies the needed second mode. This is a smaller gap but should be flagged. Third, the foundational CNO paper is cited as \"to appear\"; referees will want to see it or, failing that, a fuller derivation here.\n\nThis paper deserves a serious referee. The formal core is sound, the TBG application is a useful model-level result, and the gaps are addressable rather than fatal. The revision should mostly recalibrate claims, separating theorem from observation, and either prove or soften the transfer from Lambda to Lambda_phys. For people building local-orbital models of interacting topological or moir\\'e bands, and for DMFT or ED practitioners wanting a principled impurity-orbital choice, this is worth engaging with.","headline":"Solid methods paper with a clean SVD core and a useful chiral-TBG demonstration; the interaction-hierarchy claim is real but established only numerically for one model, so the conclusions need recalibration.","tokens_in":25324,"tokens_out":2824,"would_cite":true,"duration_ms":27411,"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":"Cell natural orbitals give every topological band a ranked set of local orbitals, with Chern topology forcing the leading orbital to vanish somewhere on the Brillouin zone.","keywords":["cell natural orbitals","topological bands","quantum geometry","density form factors","unit-cell reduced density matrix","topological zero","moiré materials","twisted bilayer graphene"],"falsifier":"Compute the full physical form-factor matrix $\\Lambda^{\\mathrm{phys}}_{k,k+q}$ and the resulting onsite interaction matrix for a Chern band whose leading CNO envelope has zeros at separated momenta rather than one coincident zero, using a screened Coulomb interaction; if a single extra CNO no longer removes power-law interaction tails, or if the ordering of interaction strengths contradicts the ordering of $\\lambda_a$, the central hierarchy claim fails.","tokens_in":24285,"feed_emoji":"","tokens_out":8093,"duration_ms":70439,"temperature":0.7,"pith_summary":"Topological bands resist exponentially localized symmetric Wannier functions, so representing their interactions in real space is usually a compromise. This paper introduces cell natural orbitals (CNOs), obtained by diagonalizing the unit-cell reduced density matrix of the band projector, to make that compromise systematic. The CNO occupations $\\lambda_a$ rank how much each local orbital contributes to the band's overlaps; truncating to the largest occupations is the optimal low-rank approximation, with error equal to the discarded squared occupations. The paper proves that the leading CNO envelope must vanish on a Chern band, forcing at least two local orbitals, and shows numerically that in magic-angle chiral twisted bilayer graphene the dominant CNO is centered at the AA site and successive channels carry weaker interactions. If correct, the framework gives a controlled, improvable local-orbital basis for correlated topological materials and explains when power-law interaction tails appear.","feed_headline":"One local orbital can never capture a Chern band","feed_subtitle":"Cell Natural Orbitals rank interaction channels, so weaker channels can be treated more cheaply in twisted bilayer graphene.","key_machinery":"The central object is the unit-cell reduced density matrix $L=\\Pi P\\Pi$ (Eq. 1), the Brillouin-zone average of the band projector $P$ restricted to one unit cell by $\\Pi$. Its eigenstates are the CNOs and its eigenvalues $\\lambda_a$ are occupations. The same nonzero spectrum appears in $\\Lambda=U^\\dagger U$ through the SVD of $U=\\Pi P$, so diagonalizing the small $N_\\alpha\\times N_\\alpha$ matrix $L$ gives the singular vectors of the momentum-space overlap kernel. The argument is carried by two identities: the Eckart-Young optimality of the rank-$N_\\tau$ truncation, with error $\\sum_{a>N_\\tau}\\lambda_a^2$, and the Poincare-Hopf count that the leading envelope's zeros have total vorticity equal to the Chern number. The Gram matrix $S_k=\\sum_{a\\le N_\\tau}\\lambda_a|s_{a,k}|^2$ controls locality: where $S_k$ vanishes, the orthonormalization factor $S_k^{-1/2}$ is non-analytic and interactions develop power-law tails; adding a CNO that is finite at the zero restores exponential decay.","core_discovery":"The paper's central discovery is that the unit-cell reduced density matrix $L_{\\alpha\\beta} = \\frac{1}{N_k}\\sum_{k\\in BZ}\\sum_{n\\in P} u_{n,\\alpha}(k) u^*_{n,\\beta}(k)$ defines a set of local orbitals, the cell natural orbitals, that carry a singular-value decomposition of the band-projected overlap kernel $\\Lambda_{km,k'n} = \\frac{1}{N_k}\\langle u_{k,m}|u_{k',n}\\rangle$. Because $L$ and $\\Lambda$ are the two adjoint products of the single rectangular operator $U=\\Pi P$, they share the nonzero spectrum $\\{\\lambda_a\\}$, and the right singular vectors $s_{a,k} = \\langle u_k|\\tau_a\\rangle/\\sqrt{\\lambda_a}$ are envelope functions on the Brillouin zone. Three claims follow. First, keeping the $N_\\tau$ largest occupations is the optimal rank-$N_\\tau$ approximation to $\\Lambda$, with exact truncation error $\\sum_{a>N_\\tau}\\lambda_a^2$. Second, for a band with Chern number $C$, the leading envelope $s_{1,k}$ is a global section of a nontrivial line bundle and must vanish at points whose vorticities sum to $2\\pi C$, so $\\lambda_1<1$, $S_{\\mathrm{cell}}>0$, and no single local orbital can represent a Chern band. Third, the occupations set a hierarchy in projected interactions: in chiral twisted bilayer graphene at the magic angle, the dominant CNO sits at the AA site, a second, delocalized CNO regularizes the Gram matrix at the protected zero, power-law interaction tails disappear, and the on-site interaction strengths are ordered across CNO channels so that subdominant channels can be treated as static mean fields while the leading channel requires dynamical self-energy. The paper notes that the kernel $\\Lambda$ is not the exact SVD of the physical form-factor matrix $\\Lambda^{\\mathrm{phys}}_{k,k+q}$ entering projected interactions; the hierarchy in that physical matrix is found numerically in chiral TBG.","pith_inferences":["The truncation-error bound suggests a practical diagnostic for embedding calculations: the minimal number of impurity orbitals for a correlated band could be chosen directly from the CNO occupations, with the sum of discarded $\\lambda_a^2$ as a controlled error, rather than from an ad hoc Wannierization.","If the observed transfer of the hierarchy to the physical form factor holds beyond chiral TBG, then CNO truncation could be tested on simpler Chern bands, such as the Haldane or multifold models, by computing the onsite structure factor $F_{ab}(q)$ at rank 1 versus rank 2; a mismatch between the $\\lambda_a$ ordering and the interaction ordering would falsify the general claim.","The paper's mechanism implies that local orbitals chosen purely for real-space localization are not enough: what matters for short-ranged interactions is whether the truncated set keeps the Gram matrix bounded away from zero, reframing Wannierizability as a quantitative coverage condition on the Brillouin zone rather than exponential decay.","The unit-cell entanglement entropy and the CNO complexity measures may serve as cheap pre-screening observables for when multi-orbital physics will be needed, for example across topological phase transitions where the entropy derivative peaks."],"forward_implications":["Any band with nonzero Chern number requires at least two CNOs: the leading envelope must vanish with total vorticity equal to $2\\pi C$, so a single local orbital is structurally incapable of representing the band and the unit-cell entanglement entropy is strictly positive.","Retaining the $N_\\tau$ largest occupations is the optimal rank-$N_\\tau$ approximation to the overlap kernel, with exact error $\\sum_{a>N_\\tau}\\lambda_a^2$, giving a quantitative stopping rule for building local-orbital models of projected interactions.","In chiral twisted bilayer graphene at the magic angle, including the second CNO regularizes the Gram matrix at the protected zero, turning power-law interaction tails into exponential decay while preserving the hierarchy of on-site interaction strengths.","The single-particle spectral function inherits momentum dependence from the CNO band weights $\\lambda_a|s_{a,k}|^2$: at the leading envelope's zero the correlation scale is set by the subleading channel, producing a gapless crossing within a single chiral sector whose fate is governed by coupling between the two chiral sectors."],"supporting_citations":[{"why":"Defines cell natural orbitals as eigenstates of the unit-cell reduced density matrix and supplies the construction this paper builds on.","marker":"[1]"},{"why":"Provides the heavy-fermion model of twisted bilayer graphene whose f-orbital the dominant CNO is compared to.","marker":"[2]"},{"why":"Establishes the obstruction to exponentially localized symmetric Wannier functions in topological bands, motivating a multi-orbital local representation.","marker":"[3]"},{"why":"Gives the mathematical Wannier-obstruction result used to argue the leading CNO envelope cannot be nowhere vanishing.","marker":"[4]"},{"why":"Identifies the relation between the single-particle correlation matrix and the fermionic reduced density matrix, connecting CNO occupations to unit-cell entanglement entropy.","marker":"[51]"},{"why":"Supports the statement that zeros of overlaps with a fixed trial orbital carry winding fixed by the Chern number.","marker":"[59]"},{"why":"The TKNN result fixing the Chern number as a topological invariant, used for the vorticity sum $2\\pi C$ of the envelope zeros.","marker":"[62]"},{"why":"Provides the standard projected-interaction formalism with the physical form factors used to define CNO interaction channels.","marker":"[66]"}],"fun_headline_variants":["Chern bands resist single-orbital local representation","Cell Natural Orbitals decode interaction hierarchy in bands","Quantum geometry sets orbital requirements for Chern bands","CNO decomposition reveals optimal local basis for topology","Twisted bilayer graphene interactions ordered by cell orbitals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the hierarchy found in the periodic-embedding overlap kernel survives in the physical form factor entering projected interactions, a property the paper verifies numerically for chiral twisted bilayer graphene but does not prove for general topological bands.","fun_headline_variants_meta":{"raw":{"variants":["Chern bands resist single-orbital local representation","Cell Natural Orbitals decode interaction hierarchy in bands","Quantum geometry sets orbital requirements for Chern bands","CNO decomposition reveals optimal local basis for topology","Twisted bilayer graphene interactions ordered by cell orbitals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000431,"raw_usage":{"total_tokens":2362,"prompt_tokens":1271,"completion_tokens":1091,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":887,"completion_tokens_details":{"reasoning_tokens":1020}},"tokens_in":887,"tokens_out":1091,"duration_ms":10295,"temperature":1.0,"reasoning_tokens":1020,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:11:59.290675+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full physical form-factor matrix $\\Lambda^{\\mathrm{phys}}_{k,k+q}$ and the resulting onsite interaction matrix for a Chern band whose leading CNO envelope has zeros at separated momenta rather than one coincident zero, using a screened Coulomb interaction; if a single extra CNO no longer removes power-law interaction tails, or if the ordering of interaction strengths contradicts the ordering of $\\lambda_a$, the central hierarchy claim fails.","supporting_citations":[{"cited_title":"Thonhauser and D","cited_arxiv_id":null,"evidence_quote":"The TKNN result fixing the Chern number as a topological invariant, used for the vorticity sum $2\\pi C$ of the envelope zeros."},{"cited_title":"Paul, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the standard projected-interaction formalism with the physical form factors used to define CNO interaction channels."}],"review_version":1}