Pith. sign in

REVIEW 4 major objections 5 minor 99 references

Cell Natural Orbitals in Quantum Materials

T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper introduces cell natural orbitals, eigenstates of a unit-cell density matrix, as a way to identify the few local orbitals that best capture chosen bands, yielding short-ranged Wannier models for twisted WSe2.

desk verdict A genuinely new wavefunction-derived route to local-orbital selection, with an overreaching 'optimal' claim and some reproducibility gaps, but the core construction is sound and worth engaging. read the letter →

arxiv 2608.10059 v1 pith:M64PE7DZ submitted 2026-08-10 cond-mat.mtrl-sci cond-mat.mes-hallcond-mat.str-el

Harshitra Mahalingam , Nishchhal Verma , Daniel Muñoz-Segovia , Raquel Queiroz This is my paper · ORCID
classification cond-mat.mtrl-scicond-mat.mes-hallcond-mat.str-el
keywords cellnaturalorbitalsunit-cellone-particlereduceddensitymatrixWannierizationquantumgeometrytwistedbilayerWSe2moirématerialstopologicalbandsentanglement
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper aims to make the choice of local orbitals for a band-structure model something the wavefunctions decide, rather than a guess. The authors define the unit-cell one-particle reduced density matrix, obtained by restricting the projector onto a chosen set of bands to a single unit cell, and show that its eigenstates—the cell natural orbitals—form a local, symmetry-labelled basis whose eigenvalues rank how much each orbital participates in those bands. They argue that the subset of these orbitals with maximal total spectral weight, subject to spanning the target bands at every momentum, provides optimal trial states for building tight-binding models, even when the bands are topologically nontrivial. They demonstrate the construction on twisted bilayer WSe2, producing local models with exponentially decaying hoppings whose Berry curvature tracks the continuum model across twist angles. If correct, the construction turns a routine object—the band projector—into a systematic diagnostic of multi-orbital character and a practical tool for correlated-materials modeling.

What carries the argument

The central object is the unit-cell one-particle reduced density matrix $L = \Pi P \Pi$, where $P$ projects onto the target Bloch bands and $\Pi$ restricts to a single unit cell; its eigenstates are the cell natural orbitals and its eigenvalues are their occupations. Equivalently, through the singular value decomposition of $U = \Pi P$, each occupation $\lambda_a$ is the squared singular value connecting a CNO $|\tau_a\rangle$ to a band-subspace state $|s_a\rangle$ whose rescaled Bloch coefficients $s_{a,n,k} = \langle u^{\mathrm{per}}_{kn}|\tau_a\rangle/\sqrt{\lambda_a}$ form the envelope functions that carry the zeros diagnosing topological or symmetry obstructions. The construction also supplies the completeness criterion $\det S_k \neq 0$, with $S_k$ the overlap Gram matrix between CNO Bloch states and target bands, and the ranking functional $\lambda_A$ for choosing among candidate CNO sets. These ingredients together convert the band projector into trial states for the single-shot Wannierization procedure used to obtain explicit tight-binding models.

What would settle it

Vary the disentanglement weighting $\eta_{km}$ and the plane-wave cutoff $N_G$ in the twisted WSe2 construction and check whether the exponential decay length of the hoppings or the Berry curvature agreement changes materially; alternatively, solve the interacting problem in the CNO-derived model and see whether a lower-$\lambda_A$ set with more local interactions reproduces the continuum many-body spectrum better than the maximal-$\lambda_A$ set.

Watch

Extended reading notes

Core claim

The central claim is that the eigenstates $|\tau_a\rangle$ of $L = \Pi P \Pi$, called cell natural orbitals, are the local degrees of freedom intrinsically selected by a target band manifold: they are uniquely determined by the Bloch wavefunctions and the real-space partition, carry site-symmetry irreps and Wyckoff positions, and their eigenvalues $\lambda_a$ measure the occupation of each orbital in the manifold. For Chern bands the projected envelope functions $s_{a,n,k}$ necessarily contain zeros, so no single compact orbital can span the manifold; these zeros are irremovable, while symmetry-enforced zeros can be lifted by moving the partition to the appropriate Wyckoff position. A CNO subset that spans the target bands (det $S_k \neq 0$ at every $k$) and maximizes the total spectral weight $\lambda_A = \mathrm{Tr}(P P_A)/N_k$ is proposed as the optimal trial-state set for single-shot Wannierization. Applied to twisted bilayer WSe2 at 3°, both an MM-centered three-orbital set and a multi-centered set produce exponentially localized Wannier functions and Wannier models whose Berry curvature reproduces the continuum result, with the multi-centered set carrying the higher total spectral weight (0.899 versus 0.755).

Load-bearing premise

The claim that a CNO set is optimal rests on equating 'best' with maximal total spectral weight plus symmetry compatibility, a proxy the paper never compares against an independent measure of whether a local model is good for correlated calculations.

Editorial extensions

If this is right

  • CNO occupations give a quantitative, projector-only measure of multi-orbital character for any set of bands, applicable to both lattice and continuum models.
  • For Chern bands, no single compact local orbital can span the manifold: the leading CNO envelope necessarily vanishes somewhere in the Brillouin zone, a real-space signature of the Wannier obstruction.
  • Symmetry-enforced zeros in CNO envelopes are removable by recentering the unit-cell partition at the appropriate Wyckoff position, which can reduce the number of orbitals needed to describe an obstructed atomic insulator.
  • In twisted bilayer WSe2 at 3°, the CNO trial states produce exponentially decaying Wannier functions and short-ranged tight-binding models whose Berry curvature matches the continuum result for the top band.
  • Because the only input is a projector, the same construction carries over to first-principles wavefunctions, so the occupation spectrum could become a routine companion to a band-structure calculation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the paper's optimality metric is not benchmarked; a natural test is to compare CNO-selected models against many-body spectra or the locality of interaction matrix elements, not only band dispersion and Berry curvature.
  • Editorial inference: the final localization depends on the manually chosen disentanglement function $\eta_{km}$ for bands outside the target manifold and on the plane-wave cutoff $N_G$; CNOs inform but do not fix these choices, so the practical 'single-shot' claim is weaker than the construction itself.
  • Editorial inference: because the construction takes only the projector as input, it could be added as a routine diagnostic in automated Wannierization pipelines, flagging band inversions through the momentum positions of envelope-function zeros.
  • Editorial inference: the paper's closing speculation that the dominant CNO favors charge-density-wave order while subdominant modes stabilize fractional Chern states is testable by solving the interacting problem in the CNO basis, which would directly tie the occupation hierarchy to interaction scales.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript introduces the unit-cell one-particle reduced density matrix L = ΠPΠ for a set of target Bloch bands and defines its eigenstates as cell natural orbitals (CNOs). It shows that CNO occupations λ_a quantify overlap with the target manifold, relate to unit-cell entanglement entropy via the singular value decomposition of U = ΠP, and carry Wyckoff-position and site-symmetry labels when the partition respects symmetry. The paper proves that Chern bands force zeros in CNO envelope functions and illustrates symmetry-enforced zeros in a monolayer TMD model, then proposes a construction of tight-binding models by selecting CNO subsets that span the target bands and maximize the spectral weight λ_A (Eq. 28), followed by single-shot Wannierization. The method is applied to continuum twisted bilayer WSe2 at θ = 3°, yielding two three-orbital local models whose hoppings decay exponentially and whose Berry curvature approximately reproduces the continuum result.

Significance. The construction is elegant and largely rigorous: L is Hermitian positive semidefinite, the SVD relations in Sec. II B are clean, and the symmetry arguments in Sec. III provide a genuine real-space diagnostic for Wannier obstructions. The explicit BHZ and monolayer TMD examples are convincing, and the tWSe2 demonstration addresses a timely problem. If the optimality claim can be properly benchmarked, the method would give a useful, parameter-light alternative to manual Wannier orbital selection. The main limitation is that the 'optimal' claim rests on a spectral-weight metric that is not validated against independent measures of model quality, and the numerical demonstration lacks stated convergence parameters.

major comments (4)
  1. [Sec. IV A / Abstract] The paper equates the 'best' local orbital set with a set maximizing the total spectral weight λ_A (Eq. 28) plus symmetry compatibility, but this quality metric is never benchmarked against an independent measure of model quality such as locality of interaction matrix elements, many-body spectra, or convergence of observables with respect to retained orbitals. As written, the abstract's 'optimal trial states' claim is therefore an assertion rather than a demonstrated result. I request either a relaxation of the wording ('well-suited' or 'systematically ranked') or an additional numerical test that compares the λ_A ranking with an independent criterion.
  2. [Sec. V B / Fig. 7] The comparison between the MM-centered and multi-centered CNO sets does not isolate the effect of λ_A. The two constructions use different disentanglement weightings η_km outside P — the MM set needs η_km 'slow enough to draw weight from the third and fourth bands,' while the multi-centered set uses a different slow decay — and the multi-centered set places orbitals on three Wyckoff sites, which trivially reduces the real-space spread. I ask for a controlled comparison: either use the same disentanglement function for both sets, or report results for several η_km choices and show that the relative ranking and spreads are stable.
  3. [Sec. IV C / Sec. V] The numerical results are not reproducible without the plane-wave cutoff N_G. Section IV C states that N_G 'must be converged' but gives no value and no convergence study, while the exponential decay of hoppings (Fig. 7e) and the Berry curvature comparison (Fig. 7d,f) depend on the real-space grid resolution. Please state the N_G used and provide a convergence test showing that the reported spreads, hopping decay constants, and Berry curvature profiles are converged.
  4. [Sec. IV B / Eq. (31)] Because single-shot Wannierization with N_τ = N_P retains the target bands exactly, reproducing the continuum dispersion (Fig. 7c) is a consistency check and not evidence that the CNO-selected model is better than an arbitrary trial set. The paper acknowledges that dispersion is a weaker test, but the Berry curvature comparison shows a 'slight' deviation in the second band that is not quantified. I recommend reporting a quantitative deviation (e.g., relative L2 error of Ω(k) or of the quantum metric) for both models, so the reader can evaluate the fidelity claim.
minor comments (5)
  1. [Sec. I] There is a typo in the Introduction: 'translational symetry' should be 'translational symmetry'.
  2. [Sec. II A] Equation (9) uses the notation ⟨N_cell⟩ but this quantity is not defined before use; please define it explicitly (the text says it is the total electron number per unit cell in the target bands, but the notation is introduced only in the equation).
  3. [Sec. IV B] In the text following Eq. (32), η_km is called 'the projector' in one place ('the projector ηkm must extend beyond P'); since η_km is a weighting factor, not a projection operator, the phrasing should be corrected.
  4. [Fig. 7 caption] The figure caption appears to contain duplicated or rearranged panel labels (c) and e) appear twice, and some labels are inconsistent between the top and bottom panels); please make the panel callouts precise.
  5. [Sec. VI] The conclusion states that the leading CNO 'defines an impurity space fixed by the wavefunctions' and that 'the occupation hierarchy orders the interaction scales,' but these are expectations rather than results shown here; I suggest marking them clearly as outlook items or citing the companion work [49] with a concrete summary.

Circularity Check

1 steps flagged · score 4.0 of 10

CNO construction is self-contained; the only circular element is using the exact-retention property of single-shot Wannierization as evidence that the model reproduces the continuum dispersion.

  1. fitted input called prediction [Sec. IV B (Single-Shot Wannierization); Sec. V B, Fig. 7c]
    "When the selected CNOs span the target bands at every k, they can serve as trial states for Wannierization [19], which retains the target bands exactly and uses the CNOs only to fix the gauge. ... Reproducing the dispersion is a weaker test than reproducing the wavefunctions"

    Eq. (31) defines the Wannier functions as a unitary rotation of the target Bloch states, and the paper states that this retains the target bands exactly. The Wannier-basis Hamiltonian is therefore unitarily equivalent to the continuum Hamiltonian on the target subspace, so the model dispersion in Fig. 7c coincides with the input band structure by construction for any spanning trial set. Presenting this overlay as evidence for the CNO-based construction is circular: it validates only the SVD gauge-fixing, not the orbital selection. The paper itself downgrades the dispersion test as a 'weaker test', but the conclusion still lists reproducing the dispersion among the model's achievements.

full rationale

The construction itself is not circular: L = ΠPΠ in Eq. (4) is built purely from the target-band projector and a spatial partition, and Eqs. (5)-(7) define CNOs and occupations with no fitted constants. The symmetry and topology statements (Secs. III B-III C) rest on external theorems rather than on the authors' own prior results. Self-citations such as [49] are not load-bearing for the derivation. The 'optimal trial states' criterion in Eq. (28) is an unbenchmarked proposal, which is a correctness risk rather than a circularity, and the body in fact ranks candidates by both spectral weight and Wannier spread. The one genuine circular element is the tWSe2 validation: the single-shot procedure is stated to retain the target bands exactly, so the dispersion overlay in Fig. 7c is a mathematical identity, not independent evidence. Because the localization and Berry-curvature comparisons are non-identical outputs and depend on the CNO trial states, the central claim retains independent content, giving a score of 4 rather than higher.

Assumptions & free parameters 3 free parameters · 5 assumptions · 2 invented entities

The central construction rests on standard linear algebra and free-fermion results (the UC-1pRDM is the Peschel correlation matrix restricted to a cell). No fitted constants enter the definition of CNOs themselves. The paper's leverage comes from four non-trivial inputs: (1) the Gaussian-state and Wick-theorem identification of the λ_a spectrum with the unit-cell entanglement spectrum (standard, cited); (2) the Chern-section theorem used for topology-enforced zeros (standard, cited); (3) the physical-embedding assumption needed for exponential localization (stated but not proven); (4) the continuum model of tWSe2 with Ref. [93] parameters, which all application results inherit. Two hand-tuned elements enter the application: the disentanglement weight η_km, chosen by inspection, and the real-space grid cutoff N_G, unreported. The spectral-weight optimality metric (Eq. 28) is an unvalidated design choice. No new physical entity is postulated; CNOs are a new basis, not a new particle or force, and they carry no falsifiable handle outside the paper.

free parameters (3)
  • Disentanglement weighting η_km for bands outside P = Not given numerically; 'an exponentially decaying function in band energies', decay set by inspection.
    In the topological case N_τ > N_P (Sec. IV B), single-shot Wannierization requires η_km for bands outside the target manifold. Its decay rate is chosen by hand ('slow enough to allow for weight from both the third and fourth bands', Sec. V B) and directly affects the Wannier spread and Berry curvature agreement reported in Fig. 7.
  • Plane-wave cutoff / real-space grid size N_G = Not reported in the text.
    For continuum models, Sec. IV C defines CNOs on a real-space grid derived from the reciprocal cutoff and states that N_G 'must be converged'; no value or convergence study is given for the tWSe2 results.
  • Unit-cell partition center (Wyckoff position) = 1a and 1c for monolayer WSe2; MM, MX and XM for tWSe2; swept r_c for the parabolic model, all hand-chosen.
    The CNO spectrum, envelope zeros, and local model all depend on the partition center (Secs. II C, III B, IV D). The paper treats the center as a variational knob (minimizing S_cell) and the tWSe2 comparison between the MM-centered and multi-centered models depends on this choice.
assumptions (5)
  • standard math The target manifold is a Fermi sea of Bloch states and the unit-cell reduced density matrix factorizes over CNO occupation modes (Gaussian-state and Wick's theorem structure).
    Used in Sec. II C (Eqs. 19-21) to identify λ_a with the unit-cell entanglement spectrum; standard result for free fermions (Peschel, Refs. [47,48]).
  • standard math A nonzero Chern number implies no smooth nowhere-vanishing section of the Bloch line bundle over the BZ torus, so any fixed-vector overlap ⟨u_k^per|τ_a⟩ must vanish somewhere.
    Invoked in Sec. III C for the topology-enforced zeros; cited to Refs. [60-63]. Standard theorem.
  • domain assumption The physical embedding with orbital positions r_α and the embedding matrix V_k captures the locality needed for exponential localization of CNOs.
    Eqs. (5)-(6): the paper states that omitting V_k would give eigenstates that are not exponentially localized (Sec. II A). The physical-embedding choice is an input assumption about how orbitals sit in the cell.
  • domain assumption The continuum model for tWSe2 with parameters of Ref. [93] correctly describes the bands of interest.
    All tWSe2 results (CNO spectra, local models, Berry curvatures) inherit the Devakul-Crépel-Zhang-Fu continuum model; the CNO tool only re-represents whatever bands this model produces.
  • domain assumption Non-singular overlap matrix det S_k ≠ 0 at all k, plus the chosen η_km, yields exponentially localized Wannier functions in the over-complete (N_τ > N_P) case.
    Secs. IV B and V B assert this ('a finite det S_k guarantees an exponentially localized set of Wannier functions'); the proof is standard only for the square N_τ = N_P case, and the over-complete disentanglement case additionally requires smoothness of the selected sub-bundle, which is not analyzed.
invented entities (2)
  • Cell natural orbitals (CNOs)
    purpose: Local orbital basis selected by the target-band projector; used as ranked trial states for Wannierization and as the eigenbasis of the unit-cell correlation matrix.
    The paper's central new object. Its usefulness is demonstrated only through internal consistency checks: the Wannier models reproduce the input bands' dispersion by construction, and the Berry curvature agreement depends on the tuned η_km. There is no falsifiable prediction that could confirm or refute the CNO basis from outside the paper.
  • Envelope functions s_{a,n,k}
    purpose: Momentum-resolved overlap of each CNO with the target bands; their zeros diagnose topological and symmetry obstructions.
    A derived quantity (Eq. 16), not an independent observable. The claim that its zeros are irremovable for Chern bands follows from a standard theorem (Sec. III C), but the diagnostic value is shown only on the paper's own examples (BHZ, WSe2, tWSe2), with no experimental or external ab initio confirmation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Cell Natural Orbitals in Quantum Materials." pith.science (2026). https://pith.science/paper/M64PE7DZ

@misc{pith2026260810059,
  author       = {Pith},
  title        = {Pith review of: Cell Natural Orbitals in Quantum Materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/M64PE7DZ}},
  note         = {Machine review of arXiv:2608.10059}
}
abstract

Understanding correlated quantum matter starts with an accurate model of the single-particle states that interact at low energies: their dispersion, band geometry, orbital content and charge density. In many cases, notably the topological bands of moire materials, it is not straightforward to find a real-space description with a few local orbitals that accomplishes this task. Here we provide a systematic way to identify the local degrees of freedom that best capture the band geometry and charge density of any chosen set of bands. We use the unit-cell one-particle reduced density matrix (UC-1pRDM), obtained by restricting the projector onto the target bands to a single unit cell. Its eigenstates, which we call cell natural orbitals (CNOs), form a local, symmetric basis uniquely determined by the Bloch wavefunctions and the choice of real-space partition. Their eigenvalues measure the occupation of each CNO in the target bands, quantifying entanglement across unit-cell boundaries and the importance of multi-orbital character. A set of CNOs that maximizes total spectral weight and reproduces the target band symmetries provides optimal trial states for Wannierization. We exemplify this by constructing a lattice model for twisted bilayer WSe$_2$ that tracks the orbital content across twist angles.

Figures

Figures reproduced from arXiv: 2608.10059 by the authors.

Figure 1
Figure 1. FIG. 1. Schematic overview of the cell natural orbital construction. Starting from a target band projector [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. a) Triangular-lattice three-orbital monolayer TMD [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. a) CNO weights [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. a) CNO occupation spectra for the lowest band of [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: a) and c) shows the CNO weights for the MM￾and MX-centered partitions respectively, across the twist angle range θ ∈ [1◦ , 5 ◦ ]. The line colors label each CNO across angles. The grey shadowed region in Fig. 5a) and c) indicates the topologically trivial regime in whi…
Figure 6
Figure 6. Figure 6: with the CNO projections ⟨u per kn |τa⟩ for the domi￾nant three CNOs in the MM-centered partition (Fig. 6a)) and in the MX-centered partition (Fig. 6c)). We first consider the MM-centered partition to con￾struct the UC-1pRDM of the top two bands, shown in Fig. 6a). The…
Figure 7
Figure 7. Figure 7: FIG. 7. Local real-space models obtained from the MM-centered CNO subset (top panel) and the multi-centered CNO subset [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

99 extracted references · 51 canonical work pages

  1. [1]

    T¨ orm¨ a, Essay: Where can quantum geometry lead us?, Phys

    P. T¨ orm¨ a, Essay: Where can quantum geometry lead us?, Phys. Rev. Lett.131, 240001 (2023)

    2023

  2. [2]

    J. Yu, B. A. Bernevig, R. Queiroz, E. Rossi, P. T¨ orm¨ a, and B.-J. Yang, Quantum geometry in quantum materi- als, npj Quantum Materials10, 101 (2025)

    2025

  3. [3]

    Rossi, Quantum metric and correlated states in two- dimensional systems, Current Opinion in Solid State and Materials Science25, 100952 (2021)

    E. Rossi, Quantum metric and correlated states in two- dimensional systems, Current Opinion in Solid State and Materials Science25, 100952 (2021)

    2021

  4. [4]

    E. Y. Andrei, D. K. Efetov, P. Jarillo-Herrero, A. H. MacDonald, K. F. Mak, T. Senthil, E. Tutuc, A. Yazdani, and A. F. Young, The marvels of moir´ e materials, Nature Reviews Materials6, 201 (2021)

    2021

  5. [5]

    K. P. Nuckolls and A. Yazdani, A microscopic perspec- tive on moir´ e materials, Nature Reviews Materials9, 460 (2024)

    2024

  6. [6]

    Aoki, Flat bands in condensed-matter systems – per- spective for magnetism and superconductivity, Contem- porary Physics66, 1 (2025)

    H. Aoki, Flat bands in condensed-matter systems – per- spective for magnetism and superconductivity, Contem- porary Physics66, 1 (2025)

    2025

  7. [7]

    J. P. Provost and G. Vallee, Riemannian structure on manifolds of quantum states, Communications in Math- ematical Physics76, 289 (1980)

    1980

  8. [8]

    A. Gao, N. Nagaosa, N. Ni, and S.-Y. Xu, Quantum Ge- ometry Phenomena in Condensed Matter Systems, arXiv (2025), 2508.00469

    arXiv 2025

Show all 99 references
  1. [9]

    D. Vanderbilt, Berry Phases in Electronic Structure Theory: Electric Polarization, Orbital Magnetization and Topological Insulators, Cambridge University Press 10.1017/9781316662205 (2018)

    2018 doi

  2. [10]

    Ingham and R

    J. Ingham and R. Queiroz, Real-Space Imaging of Band Topology via Wavefunction Zeros, arxiv , 2607.21699 (2026)

    2026 arXiv

  3. [11]

    Verma and R

    N. Verma and R. Queiroz, Instantaneous response and quantum geometry of insulators, Proceed- ings of the National Academy of Sciences122, 10.1073/pnas.2405837122 (2025)

    2025 doi

  4. [12]

    Verma, P

    N. Verma, P. J. W. Moll, T. Holder, and R. Queiroz, Quantum geometry and the hidden scales in materials, Nature Reviews Physics8, 226 (2026)

    2026

  5. [13]

    Morales-Dur´ an, N

    N. Morales-Dur´ an, N. C. Hu, P. Potasz, and A. H. Mac- Donald, Nonlocal interactions in moir´ e hubbard systems, Phys. Rev. Lett.128, 217202 (2022)

    2022

  6. [14]

    F. Wu, T. Lovorn, E. Tutuc, I. Martin, and A. H. MacDonald, Topological insulators in twisted transition metal dichalcogenide homobilayers, Phys. Rev. Lett.122, 086402 (2019)

    2019

  7. [15]

    Cr´ epel and A

    V. Cr´ epel and A. Millis, Bridging the small and large in twisted transition metal dichalcogenide homobilayers: A tight binding model capturing orbital interference and topology across a wide range of twist angles, Physical Review Research6, 033127 (2024)

    2024

  8. [16]

    W.-X. Qiu, B. Li, X.-J. Luo, and F. Wu, Interaction- driven topological phase diagram of twisted bilayer mote2, Phys. Rev. X13, 041026 (2023)

    2023

  9. [17]

    J. Zang, J. Wang, J. Cano, and A. J. Millis, Hartree- fock study of the moir´ e hubbard model for twisted bi- layer transition metal dichalcogenides, Phys. Rev. B104, 075150 (2021)

    2021

  10. [18]

    J. D. Cloizeaux, Analytical Properties of n -Dimensional Energy Bands and Wannier Functions, Physical Review 135, A698 (1964)

    1964

  11. [19]

    Marzari, A

    N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Maximally localized wannier functions: Theory and applications, Reviews of Modern Physics84, 1419 (2012)

    2012

  12. [20]

    Souza, N

    I. Souza, N. Marzari, and D. Vanderbilt, Maximally local- ized wannier functions for entangled energy bands, Phys. Rev. B65, 035109 (2001)

    2001

  13. [21]

    Damle, L

    A. Damle, L. Lin, and L. Ying, Compressed represen- tation of Kohn–Sham orbitals via selected columns of the density matrix, J. Chem. Theory Comput.11, 1463 (2015)

    2015

  14. [22]

    Damle, L

    A. Damle, L. Lin, and L. Ying, SCDM-k: Localized or- bitals for solids via selected columns of the density ma- trix, J. Comput. Phys.334, 1 (2017)

    2017

  15. [23]

    Vitale, G

    V. Vitale, G. Pizzi, A. Marrazzo, J. R. Yates, N. Marzari, and A. A. Mostofi, Automated high-throughput Wan- nierisation, npj Computational Materials6, 66 (2020)

    2020

  16. [24]

    Bradlyn, L

    B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Topological quantum chemistry, Nature547, 298 (2017)

    2017

  17. [25]

    J. Cano, B. Bradlyn, Z. Wang, L. Elcoro, M. G. Vergniory, C. Felser, M. I. Aroyo, and B. A. Bernevig, 15 Building blocks of topological quantum chemistry: El- ementary band representations, Physical Review B97, 035139 (2018)

    2018

  18. [26]

    Zak, Band representations of space groups, Physical Review B26, 3010 (1982)

    J. Zak, Band representations of space groups, Physical Review B26, 3010 (1982)

    1982

  19. [27]

    L¨ owdin, Quantum theory of many-particle systems

    P.-O. L¨ owdin, Quantum theory of many-particle systems. i. physical interpretations by means of density matrices, natural spin-orbitals, and convergence problems in the method of configurational interaction, Phys. Rev.97, 1474 (1955)

    1955

  20. [28]

    Monsen and M

    B. Monsen and M. Claassen, Supercell wannier functions and emergent kondo lattices in topological bands, arXiv (2024), arXiv:2412.17190

    2024 arXiv

  21. [29]

    Cole and D

    T. Cole and D. Vanderbilt, Reduced wannier represen- tation for topological bands, Phys. Rev. B111, 205139 (2025)

    2025

  22. [30]

    Fischer, L

    A. Fischer, L. Klebl, D. M. Kennes, and T. O. Wehling, Supercell Wannier functions and a faithful low-energy model for Bernal bilayer graphene, Physical Review B 110, L201113 (2024)

    2024

  23. [31]

    Kang and O

    J. Kang and O. Vafek, Symmetry, Maximally Localized Wannier States, and a Low-Energy Model for Twisted Bilayer Graphene Narrow Bands, Physical Review X8, 031088 (2018)

    2018

  24. [32]

    J. Cao, M. Wang, S.-F. Qian, C.-C. Liu, and Y. Yao, Ab initio four-band wannier tight-binding model for generic twisted graphene systems, Phys. Rev. B104, L081403 (2021)

    2021

  25. [33]

    H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, Faithful tight-binding models and fragile topology of magic-angle bilayer graphene, Physical Review B99, 195455 (2019)

    2019

  26. [34]

    S. Carr, S. Fang, H. C. Po, A. Vishwanath, and E. Kaxi- ras, Derivation of Wannier orbitals and minimal-basis tight-binding Hamiltonians for twisted bilayer graphene: First-principles approach, Physical Review Research1, 033072 (2019)

    2019

  27. [35]

    M. J. Calder´ on and E. Bascones, Interactions in the 8- orbital model for twisted bilayer graphene, Physical Re- view B102, 155149 (2020)

    2020

  28. [36]

    Song and B

    Z.-D. Song and B. A. Bernevig, Magic-angle twisted bi- layer graphene as a topological heavy fermion problem, Phys. Rev. Lett.129, 047601 (2022)

    2022

  29. [37]

    C˘ alug˘ aru, M

    D. C˘ alug˘ aru, M. Borovkov, L. L. H. Lau, P. Coleman, Z.- D. Song, and B. A. Bernevig, Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters, Low Temperature Physics49, 640 (2023)

    2023

  30. [38]

    F. Xie, C. Li, J. Cano, and Q. Si, Kondo-lattice phenomenology of twisted bilayer wse 2 from compact molecular orbitals of topological bands, arxiv (2025), arXiv:2503.21769 [cond-mat.str-el]

    2025 arXiv

  31. [39]

    F. Xie, L. Chen, S. Sur, Y. Fang, J. Cano, and Q. Si, Su- perconductivity in twisted wse 2 from topology-induced quantum fluctuations, Phys. Rev. Lett.134, 136503 (2025)

    2025

  32. [40]

    F. Xie, L. Chen, Y. Fang, and Q. Si, Topological kondo semimetals emulated in heterobilayer transition metal dichalcogenides, Phys. Rev. Res.7, 033093 (2025)

    2025

  33. [41]

    Mahankali, F

    M. Mahankali, F. Xie, Y. Fang, L. Chen, S. Sur, S. Paschen, J. C. Souza, M. Haim, A. Gupta, N. Avra- ham, H. Beidenkopf, H. Tan, B. Yan, and Q. Si, Cor- related flat-band physics in a bilayer kagome metal based on compact molecular orbitals, arxiv (2025), arXiv:2503.09706 [con...

    2025 arXiv

  34. [42]

    C. Li, R. Dutta, F. Xie, J. R. Chelikowsky, J. Cano, M. H. Naik, and Q. Si, Topology and compact molec- ular orbitals in twisted bilayer wse 2, arxiv (2026), arXiv:2607.15611

    2026 arXiv

  35. [43]

    F. D. M. Haldane, Attachment of surface ”fermi arcs” to the bulk fermi surface: ”fermi-level plumbing” in topo- logical metals, arXiv (2014), 1401.0529

    2014 arXiv

  36. [44]

    L. K. Lim, J. N. Fuchs, and G. Montambaux, Geome- try of Bloch states probed by St¨ uckelberg interferome- try, Physical Review A - Atomic, Molecular, and Optical Physics92, 063627 (2015)

    2015

  37. [45]

    Herzog-Arbeitman, V

    J. Herzog-Arbeitman, V. Peri, F. Schindler, S. D. Hu- ber, and B. A. Bernevig, Superfluid Weight Bounds from Symmetry and Quantum Geometry in Flat Bands, Phys- ical Review Letters128, 087002 (2022)

    2022

  38. [46]

    K. E. Huhtinen, J. Herzog-Arbeitman, A. Chew, B. A. Bernevig, and P. T¨ orm¨ a, Revisiting flat band supercon- ductivity: Dependence on minimal quantum metric and band touchings, Physical Review B106, 014518 (2022)

    2022

  39. [47]

    Peschel, Calculation of reduced density matrices from correlation functions, J

    I. Peschel, Calculation of reduced density matrices from correlation functions, J. Phys. A: Math. Gen.36, L205 (2003)

    2003

  40. [48]

    Peschel and V

    I. Peschel and V. Eisler, Reduced density matrices and entanglement entropy in free lattice models, J. Phys. A: Math. Theor.42, 504003 (2009)

    2009

  41. [49]

    Verma, H

    N. Verma, H. Mahalingam, D. Mu˜ noz Segovia, and R. Queiroz, Cell natural orbitals in interacting topologi- cal bands, arxiv , to appear (2026)

    2026

  42. [50]

    Y. Xu, L. Elcoro, Z.-D. Song, B. J. Wieder, M. G. Vergniory, N. Regnault, Y. Chen, C. Felser, and B. A. Bernevig, High-throughput calculations of mag- netic topological materials, Nature586, 702 (2020)

    2020

  43. [51]

    B. J. Wieder, B. Bradlyn, J. Cano, Z. Wang, M. G. Vergniory, L. Elcoro, A. A. Soluyanov, C. Felser, T. Ne- upert, N. Regnault, and B. A. Bernevig, Topological ma- terials discovery from crystal symmetry, Nature Reviews Materials7, 196 (2022)

    2022

  44. [52]

    H. C. Po, A. Vishwanath, and H. Watanabe, Complete theory of symmetry-based indicators of band topology, Nature Communications8, 50 (2017)

    2017

  45. [53]

    J. Cano, L. Elcoro, M. I. Aroyo, B. A. Bernevig, and B. Bradlyn, Topology invisible to eigenvalues in ob- structed atomic insulators, Phys. Rev. B105, 125115 (2022)

    2022

  46. [54]

    Liu, W.-Y

    G.-B. Liu, W.-Y. Shan, Y. Yao, W. Yao, and D. Xiao, Three-band tight-binding model for monolayers of group- vib transition metal dichalcogenides, Phys. Rev. B88, 085433 (2013)

    2013

  47. [55]

    J. Zeng, H. Liu, H. Jiang, Q. F. Sun, and X. C. Xie, Multiorbital model reveals a second-order topological in- sulator in 1H transition metal dichalcogenides, Physical Review B104, L161108 (2021)

    2021

  48. [56]

    S. Qian, G. B. Liu, C. C. Liu, and Y. Yao, Cn -symmetric higher-order topological crystalline insulators in atomi- cally thin transition metal dichalcogenides, Physical Re- view B105, 045417 (2022)

    2022

  49. [57]

    Jung and Y

    J. Jung and Y. H. Kim, Hidden breathing kagome topol- ogy in hexagonal transition metal dichalcogenides, Phys- ical Review B105, 085138 (2022)

    2022

  50. [58]

    P. Eck, Y. Fang, D. Di Sante, G. Sangiovanni, and J. Cano, Recipe for higher order topology on the trian- gular lattice, Physical Review B107, 115130 (2023)

    2023

  51. [59]

    Holbrook, J

    M. Holbrook, J. Ingham, D. Kaplan, L. N. Holtzman, B. Bierman, B. Hou, N. Olsen, L. Nashabeh, Y. Li, S. Liu, 16 X. Zhu, D. Y. Qiu, D. Rhodes, K. Barmak, J. C. Hone, R. Queiroz, and A. N. Pasupathy, Real-space imaging of the band topology of transition metal dichalcogenides, Na...

    2026

  52. [60]

    Brouder, G

    C. Brouder, G. Panati, M. Calandra, C. Mourougane, and N. Marzari, Exponential localization of wannier func- tions in insulators, Phys. Rev. Lett.98, 046402 (2007)

    2007

  53. [61]

    Monaco, G

    D. Monaco, G. Panati, A. Pisante, and S. Teufel, Optimal decay of wannier functions in chern and quantum hall in- sulators, Communications in Mathematical Physics359, 61 (2018)

    2018

  54. [62]

    Mera and T

    B. Mera and T. Ozawa, Singular-connection approach to topological phases and resonant optical responses, Phys. Rev. B106, 245134 (2022)

    2022

  55. [63]

    Thonhauser and D

    T. Thonhauser and D. Vanderbilt, Insulator/chern- insulator transition in the haldane model, Physical Re- view B74, 235111 (2006)

    2006

  56. [64]

    Fischer, L

    A. Fischer, L. Klebl, V. Cr´ epel, S. Ryee, A. Rubio, L. Xian, T. O. Wehling, A. Georges, D. M. Kennes, and A. J. Millis, Theory of Intervalley-Coherent AFM Order and Topological Superconductivity in tWSe2, Physical Review X15, 041055 (2025)

    2025

  57. [65]

    Mu˜ noz-Segovia, V

    D. Mu˜ noz-Segovia, V. Cr´ epel, R. Queiroz, and A. J. Mil- lis, Twist-angle evolution of the intervalley-coherent an- tiferromagnet in twisted WSe2, Physical Review B112, 085111 (2025)

    2025

  58. [66]

    Wang, E.-M

    L. Wang, E.-M. Shih, A. Ghiotto, L. Xian, D. A. Rhodes, C. Tan, M. Claassen, D. M. Kennes, Y. Bai, B. Kim, K. Watanabe, T. Taniguchi, X. Zhu, J. Hone, A. Ru- bio, A. N. Pasupathy, and C. R. Dean, Correlated elec- tronic phases in twisted bilayer transition metal dichalco- geni...

    2020

  59. [67]

    Y. Tang, L. Li, T. Li, Y. Xu, S. Liu, K. Barmak, K. Watanabe, T. Taniguchi, A. H. MacDonald, J. Shan, and K. F. Mak, Simulation of Hubbard model physics in WSe2/WS2 moir´ e superlattices, Nature 2020 579:7799 579, 353 (2020)

    2020

  60. [68]

    Anderson, F

    E. Anderson, F. R. Fan, J. Cai, W. Holtzmann, T. Taniguchi, K. Watanabe, D. Xiao, W. Yao, and X. Xu, Programming correlated magnetic states with gate-controlled moir´ e geometry, Science381, 325 (2023)

    2023

  61. [69]

    Ghiotto, L

    A. Ghiotto, L. Wei, L. Song, J. Zang, A. B. Tazi, D. Os- trom, K. Watanabe, T. Taniguchi, J. C. Hone, D. A. Rhodes, A. J. Millis, C. R. Dean, L. Wang, and A. N. Pa- supathy, Stoner instabilities and ising excitonic states in twisted transition metal dichalcogenides, arXiv (202...

    2024 arXiv

  62. [70]

    Kn¨ uppel, J

    P. Kn¨ uppel, J. Zhu, Y. Xia, Z. Xia, Z. Han, Y. Zeng, K. Watanabe, T. Taniguchi, J. Shan, and K. F. Mak, Correlated states controlled by a tunable van Hove sin- gularity in moir´ e WSe2 bilayers, Nat Commun16, 1959 (2025)

    2025

  63. [71]

    Y. Xu, S. Liu, D. A. Rhodes, K. Watanabe, T. Taniguchi, J. Hone, V. Elser, K. F. Mak, and J. Shan, Correlated in- sulating states at fractional fillings of moir´ e superlattices, Nature587, 214 (2020)

    2020

  64. [72]

    E. C. Regan, D. Wang, C. Jin, M. I. Bakti Utama, B. Gao, X. Wei, S. Zhao, W. Zhao, Z. Zhang, K. Yu- migeta, M. Blei, J. D. Carlstr¨ om, K. Watanabe, T. Taniguchi, S. Tongay, M. Crommie, A. Zettl, and F. Wang, Mott and generalized Wigner crystal states in WSe2/WS2 moir´ e super...

    2020

  65. [73]

    Huang, T

    X. Huang, T. Wang, S. Miao, C. Wang, Z. Li, Z. Lian, T. Taniguchi, K. Watanabe, S. Okamoto, D. Xiao, S.- F. Shi, and Y.-T. Cui, Correlated insulating states at fractional fillings of the WS 2/WSe2 moir´ e lattice, Nat. Phys.17, 715 (2021)

    2021

  66. [74]

    H. Li, S. Li, E. C. Regan, D. Wang, W. Zhao, S. Kahn, K. Yumigeta, M. Blei, T. Taniguchi, K. Watanabe, S. Tongay, A. Zettl, M. F. Crommie, and F. Wang, Imag- ing two-dimensional generalized Wigner crystals, Nature 2021 597:7878597, 650 (2021)

    2021

  67. [75]

    T. Li, S. Jiang, B. Shen, Y. Zhang, L. Li, Z. Tao, T. De- vakul, K. Watanabe, T. Taniguchi, L. Fu, J. Shan, and K. F. Mak, Quantum anomalous Hall effect from in- tertwined moir´ e bands, Nature 2021 600:7890600, 641 (2021)

    2021

  68. [76]

    Z. Tao, B. Shen, S. Jiang, T. Li, L. Li, L. Ma, W. Zhao, J. Hu, K. Pistunova, K. Watanabe,et al., Valley-coherent quantum anomalous Hall state in AB-stacked MoTe2/ws e 2 bilayers, Physical Review X14, 011004 (2024)

    2024

  69. [77]

    B. A. Foutty, C. R. Kometter, T. Devakul, A. P. Reddy, K. Watanabe, T. Taniguchi, L. Fu, and B. E. Feld- man, Mapping twist-tuned multiband topology in bilayer WSe2, Science384, 343 (2024)

    2024

  70. [78]

    W. Zhao, B. Shen, Z. Tao, Z. Han, K. Kang, K. Watan- abe, T. Taniguchi, K. F. Mak, and J. Shan, Gate-tunable heavy fermions in a moir´ e kondo lattice, Nature616, 61–65 (2023)

    2023

  71. [79]

    W. Zhao, B. Shen, Z. Tao, S. Kim, P. Kn¨ uppel, Z. Han, Y. Zhang, K. Watanabe, T. Taniguchi, D. Chowdhury, J. Shan, and K. F. Mak, Emergence of ferromagnetism at the onset of moir´ e kondo breakdown, Nature Physics 20, 1772–1777 (2024)

    2024

  72. [80]

    Z. Han, Y. Xia, Z. Xia, W. Zhao, Y. Zhang, K. Watanabe, T. Taniguchi, J. Shan, and K. F. Mak, Topological kondo insulator in mote2/wse2 moir´ e bilayers, Nature Physics 22, 396–401 (2026)

    2026

  73. [81]

    W. Tian, B. Shen, L. Li, M. Zhang, F. Liu, C. Li, Y. Liu, F. Xu, K. Watanabe, T. Taniguchi, P. Li, L. Lu, Y. Xu, S. Jiang, T. Li, J. Shan, and K. F. Mak, Emergence of a symmetry-broken chern insulator near a moir´ e kondo breakdown, arxiv (2026), arXiv:2602.17349

    2026

  74. [82]

    Y. Xia, Z. Han, K. Watanabe, T. Taniguchi, J. Shan, and K. F. Mak, Superconductivity in twisted bilayer WSe2, Nature 2024 637:8047637, 833 (2024)

    2024

  75. [83]

    Y. Guo, J. Pack, J. Swann, L. Holtzman, M. Cothrine, K. Watanabe, T. Taniguchi, D. G. Mandrus, K. Barmak, J. Hone, A. J. Millis, A. Pasupathy, and C. R. Dean, Superconductivity in 5.0 ◦ twisted bilayer WSe2, Nature 2025 637:8047637, 839 (2025)

    2025

  76. [84]

    Y. Xia, Z. Han, J. Zhu, Y. Zhang, P. Kn¨ uppel, K. Watan- abe, T. Taniguchi, K. F. Mak, and J. Shan, Bandwidth- tuned Mott transition and superconductivity in moir´ e WSe2, Nature 2026 650:8102650, 585 (2026)

    2026

  77. [85]

    Y. Guo, J. Cenker, A. Fischer, D. Mu˜ noz-Segovia, J. Pack, L. Holtzman, L. Klebl, K. Watanabe, T. Taniguchi, K. Barmak, J. Hone, A. Rubio, D. M. Kennes, A. J. Millis, A. Pasupathy, and C. R. Dean, Angle evolution of the superconducting phase diagram in twisted bilayer WSe2, N...

    2026

  78. [86]

    Metzner and D

    W. Metzner and D. Vollhardt, Correlated lattice fermions ind=∞dimensions, Phys. Rev. Lett.62, 324 (1989)

    1989

  79. [87]

    Georges, G

    A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- berg, Dynamical mean-field theory of strongly correlated 17 fermion systems and the limit of infinite dimensions, Rev. Mod. Phys.68, 13 (1996)

    1996

  80. [88]

    Salmhofer and C

    M. Salmhofer and C. Honerkamp, Fermionic Renormal- ization Group Flows: Technique and Theory, Progress of Theoretical Physics105, 1 (2001)

    2001

  81. [89]

    Metzner, M

    W. Metzner, M. Salmhofer, C. Honerkamp, V. Meden, and K. Sch¨ onhammer, Functional renormalization group approach to correlated fermion systems, Reviews of Mod- ern Physics84, 299 (2012)

    2012

  82. [90]

    Platt, W

    C. Platt, W. Hanke, and R. Thomale, Functional renor- malization group for multi-orbital Fermi surface instabil- ities, Advances in Physics62, 453 (2013)

    2013

  83. [91]

    Dupuis, L

    N. Dupuis, L. Canet, A. Eichhorn, W. Metzner, J. M. Pawlowski, M. Tissier, and N. Wschebor, The nonper- turbative functional renormalization group and its appli- cations, Physics Reports910, 1 (2021)

    2021

  84. [92]

    Beyer, J

    J. Beyer, J. B. Hauck, and L. Klebl, Reference results for the momentum space functional renormalization group, The European Physical Journal B 2022 95:495, 65 (2022)

    2022

  85. [93]

    Devakul, V

    T. Devakul, V. Cr´ epel, Y. Zhang, and L. Fu, Magic in twisted transition metal dichalcogenide bilayers, Nature Communications12, 6730 (2021)

    2021

  86. [94]

    Tuo, M.-R

    C. Tuo, M.-R. Li, Z. Wu, W. Sun, and H. Yao, Theory of topological superconductivity and antiferromagnetic cor- related insulators in twisted bilayer wse2, Nature Com- munications16, 10.1038/s41467-025-64519-3 (2025)

    2025 doi

  87. [95]

    Knizia and G

    G. Knizia and G. K.-L. Chan, Density matrix embed- ding: A simple alternative to dynamical mean-field the- ory, Phys. Rev. Lett.109, 186404 (2012)

    2012

  88. [96]

    J. Karp, A. Hampel, and A. J. Millis, Dependence of DFT + DMFT results on the construction of the correlated orbitals, Physical Review B103, 195101 (2021)

    2021

  89. [97]

    S. H. Aronson, T. Han, Z. Lu, Y. Yao, J. P. Butler, K. Watanabe, T. Taniguchi, L. Ju, and R. C. Ashoori, Displacement field-controlled fractional chern insulators and charge density waves in a graphene/hbn moir´ e su- perlattice, Phys. Rev. X15, 031026 (2025)

    2025

  90. [98]

    J. P. Butler, T. Han, A. DiFabbio, Z. Hadjri, E. Aitken, K. Watanabe, T. Taniguchi, L. Ju, and R. C. Ashoori, 1/3 Fractional and Gapless Integer Quantum Anoma- lous Hall States in Rhombohedral Graphene, arxiv , 2606.06450 (2026)

    2026 arXiv

  91. [99]

    Z. Lu, T. Han, Y. Yao, A. P. Reddy, J. Yang, J. Seo, K. Watanabe, T. Taniguchi, L. Fu, and L. Ju, Fractional quantum anomalous Hall effect in multilayer graphene, Nature626, 759 (2024)

    2024

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.