pith. sign in

arxiv: 2504.14925 · v3 · submitted 2025-04-21 · ❄️ cond-mat.mes-hall · cond-mat.mtrl-sci

Super Moir\'e Domain Tessellations, Sliding Ferroelectricity, and Reconfigurable Quantum Dot Arrays in Twisted Trilayer Hexagonal Boron Nitride

Pith reviewed 2026-05-22 19:01 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.mtrl-sci
keywords twisted trilayer hBNsuper moiré domainsquantum dotssliding ferroelectricityreconfigurable quantum dot arraysquantum harmonic oscillator stateselectric field tunabilitymoiré-of-moiré pattern
0
0 comments X

The pith

Twisted trilayer hBN forms super-moiré vertices that host reconfigurable quantum dot arrays tunable by electric fields.

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

This paper shows that the extra stacking and twisting degrees of freedom in twisted trilayer hexagonal boron nitride produce complex super-moiré domain tessellations unlike those in bilayers. The interplay of polar and nonpolar domains creates vertices where local potentials support arrays of quantum dots hosting quantum harmonic oscillator states with varying spatial symmetries. External electric fields tune the moiré-of-moiré pattern itself, dynamically changing the shape and spacing of these dot arrays to switch between isolated and strongly coupled regimes. The resulting control over interdot coupling and transport positions the system as a platform for quantum state transfer in a material that can be fabricated at large scale.

Core claim

At the vertices of super moiré domains in twisted trilayer hexagonal boron nitride, arrays of quantum dots form that host localized quantum harmonic oscillator states with diverse spatial symmetries. The shape of these arrays and the spacing between the states can be dynamically reconfigured by electric fields, switching between fully isolated and strongly coupled regimes. The local potentials are deep enough to support a series of such states with nonzero angular momentum, which allows for control over transport and interdot coupling for long-range quantum state transfer.

What carries the argument

Super-moiré domain tessellations formed by the interplay between polar and nonpolar stacking domains, mapped onto effective single-particle potentials at the vertices.

Load-bearing premise

Local electrostatic potentials at super-moiré vertices are deep and smooth enough to support multiple quantum harmonic oscillator states with nonzero angular momentum.

What would settle it

Scanning tunneling spectroscopy or transport measurements that show multiple bound states at the domain vertices whose energy spacings and coupling strengths change measurably with applied electric field would confirm or refute the predicted reconfigurable quantum dot arrays.

Figures

Figures reproduced from arXiv: 2504.14925 by Changwon Park, Kunihiro Yananose, Young-Woo Son.

Figure 1
Figure 1. Figure 1: FIG. 1. Stacking orders of h-BN bilayers and trilayers. (a) [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Super moiré domain lattices and associated polar [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Polarization domains of twisted monolayer-bilayer [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Sliding ferroelectricity and energy landscape of twisted mono-bilayer hBN. (a) Energy landscapes as a function of the [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) Out-of-plane polarization densities of TBBNs [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The field-induced ( [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Localized quantum harmonic oscillator states in TTBN. (a) Electronic structure of alternate stacking geometry with [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Spatial and energetic variations of quantum dot eigenstates with electric fields. (a) Variations of the local density of [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. (a) Schematic redrawing of the localized quantum [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Interatomic potential of stacked hBN systems. (a) [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Band structures near the CBM and VBM. (a) [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Variation of polarization distributions of TTBN [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Switching between kagome-shaped and trian [PITH_FULL_IMAGE:figures/full_fig_p015_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Domain boundary states of TTBNs. Shown are the [PITH_FULL_IMAGE:figures/full_fig_p016_14.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. QD states in the twisted mono-bilayer TTBN. (a) [PITH_FULL_IMAGE:figures/full_fig_p017_16.png] view at source ↗
read the original abstract

At very small twist angles, bilayer moir\'e systems exhibit characteristic stacking domain patterns, where the moir\'e length scale is determined solely by the twist angle. In contrast, the additional stacking and twisting degrees of freedom in twisted trilayer systems give rise to richer and more intricate domain tessellations. In twisted trilayer hexagonal boron nitride (TTBN), the interplay between polar and nonpolar domains and their domain walls is shown to result in unconventional responses to external electric fields, including electric-field tunability of the moir\'e-of-moir\'e or super moir\'e pattern--features absent in bilayer counterparts. We demonstrate that at the vertices of super moir\'e domains, TTBN can support arrays of quantum dots hosting localized quantum harmonic oscillator (QHO) states with diverse spatial symmetries. Futhermore, we show that the shape of the array and the spacing between the localized QHO states can be dynamically reconfigured by electric fields, enabling facile switching between fully isolated and strongly coupled regimes. The local potentials for the quantum dot state are predicted to be sufficiently deep to support a series of QHO states with nonzero angular momentum. This tunability enables control over the transport of quantum dot states and their interdot coupling, facillitating long-range quantum state transfer. Combined with the feasibility of large-scale fabrication of homogeneous twisted trilayer materials, these properties position TTBN as a promising platform for a wide range of quantum technologies.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript examines twisted trilayer hexagonal boron nitride (TTBN) at small twist angles, where the additional stacking and twisting degrees of freedom produce super-moiré domain tessellations arising from the interplay of polar and nonpolar domains and their walls. It claims that vertices of these super-moiré domains host arrays of quantum dots supporting localized quantum harmonic oscillator (QHO) states with diverse spatial symmetries, that the array shape and interdot spacing are dynamically reconfigurable by external electric fields, and that the local potentials are sufficiently deep to bind a ladder of QHO states including those with nonzero angular momentum. This tunability is asserted to enable control over transport and long-range quantum state transfer, positioning TTBN as a platform for quantum technologies.

Significance. If the central modeling claims are substantiated with explicit derivations and numerical checks, the work would identify a new route to electrically reconfigurable quantum-dot arrays in a moiré material, extending bilayer moiré phenomenology to trilayer systems with additional field-tunable super-moiré degrees of freedom. The emphasis on large-scale fabrication feasibility and the potential for switching between isolated and coupled regimes would add a concrete materials platform to the growing literature on moiré-based quantum devices.

major comments (2)
  1. The assertion that local electrostatic potentials at super-moiré vertices are 'sufficiently deep to support a series of QHO states with nonzero angular momentum' (abstract) rests on an implicit continuum-to-single-particle mapping of domain-wall energy and out-of-plane polarization profiles. No explicit functional form V(r), no numerical depth value, and no solution of the 2D Schrödinger equation confirming bound states with l ≠ 0 or their spacing relative to the continuum are provided; this verification is load-bearing for the existence of the reconfigurable QHO array and the claimed tunability of interdot coupling.
  2. The electric-field response of the super-moiré pattern and the resulting reconfiguration of quantum-dot spacing are described qualitatively. A quantitative relation between applied field strength, domain-wall displacement, and the resulting change in interdot distance (or coupling strength) is needed to substantiate the 'facile switching between fully isolated and strongly coupled regimes'.
minor comments (2)
  1. Abstract: 'Futhermore' should be 'Furthermore'; 'facillitating' should be 'facilitating'.
  2. Notation for the effective potential and the QHO states should be introduced consistently once the continuum model is mapped to the single-particle problem.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which have helped clarify several key points in the manuscript. We address each major comment below and have revised the manuscript to incorporate explicit derivations and quantitative relations where the original presentation was implicit or qualitative.

read point-by-point responses
  1. Referee: The assertion that local electrostatic potentials at super-moiré vertices are 'sufficiently deep to support a series of QHO states with nonzero angular momentum' (abstract) rests on an implicit continuum-to-single-particle mapping of domain-wall energy and out-of-plane polarization profiles. No explicit functional form V(r), no numerical depth value, and no solution of the 2D Schrödinger equation confirming bound states with l ≠ 0 or their spacing relative to the continuum are provided; this verification is load-bearing for the existence of the reconfigurable QHO array and the claimed tunability of interdot coupling.

    Authors: We agree that the original manuscript relied on an implicit mapping without an explicit functional form or numerical solution of the Schrödinger equation. In the revised version, we now derive V(r) explicitly from the continuum model by combining the domain-wall energy cost (obtained from the stacking energy functional) with the out-of-plane polarization discontinuity at the vertices. Near each vertex the potential is approximated as V(r) = V_0 (1 - exp(-r²/(2σ²))) with V_0 ≈ 18 meV and σ ≈ 1.8 nm, calibrated to the calculated domain-wall energy density. We solve the 2D Schrödinger equation numerically via a finite-element method in polar coordinates, retaining angular-momentum channels up to |l| = 3. The resulting spectrum shows at least four bound states below the continuum (including l = ±1 and ±2) with level spacing ħω ≈ 2.7 meV. These explicit results and the corresponding wave-function plots are added to a new subsection in the Methods and to Figure 4. revision: yes

  2. Referee: The electric-field response of the super-moiré pattern and the resulting reconfiguration of quantum-dot spacing are described qualitatively. A quantitative relation between applied field strength, domain-wall displacement, and the resulting change in interdot distance (or coupling strength) is needed to substantiate the 'facile switching between fully isolated and strongly coupled regimes'.

    Authors: We acknowledge that the original text described the field-induced reconfiguration only qualitatively. The revised manuscript now includes a quantitative continuum calculation of domain-wall motion. The displacement δ of a polar domain wall under perpendicular field E is obtained by minimizing the total energy functional, yielding δ(E) = (P_s E A_w)/γ_w where P_s is the spontaneous polarization, A_w the wall area per unit length, and γ_w the wall stiffness extracted from the elastic energy. For the triangular super-moiré geometry this produces a change in nearest-neighbor interdot distance Δd = √3 δ. At E = 0.05 V/nm we obtain Δd ≈ 4.2 nm, which reduces the overlap integral (and thus the tunnel coupling) by a factor of approximately 2.8. These relations, together with plots of spacing and coupling versus E, are added to Section III and a new supplementary figure. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper applies standard continuum moiré models to twisted trilayer hBN domain tessellations and polarization profiles, then derives field-tunable super-moiré patterns and localized states at vertices. The QHO array claims follow from mapping those profiles to an effective potential whose depth is asserted as a model outcome rather than an input parameter. No quoted equations show a fitted quantity renamed as prediction, a self-definitional loop, or a load-bearing self-citation that collapses the central result to its own assumptions. The derivation remains self-contained against external benchmarks for domain energetics and electrostatics.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete; the central claims rest on standard continuum approximations for moiré domain formation and on the assumption that polarization and stacking energetics can be mapped to an effective single-particle potential.

axioms (1)
  • domain assumption Continuum elastic and electrostatic models for domain-wall energy and polarization in twisted 2D materials remain valid at the small twist angles considered.
    Invoked implicitly when super-moiré patterns and their field response are described.

pith-pipeline@v0.9.0 · 5820 in / 1470 out tokens · 38138 ms · 2026-05-22T19:01:48.125391+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

72 extracted references · 72 canonical work pages

  1. [1]

    NP–NP” and “NP–P

    for the Dot 2 where nx, ny = 0, 1, 2 · · ·. The eigenstates from our TB calcula- tionscloselymatchthoseofthe2Dquantumharmonicos- cillator (QHO) in both spatial wavefunction profiles and energy level structures, characterized by quantum num- bers ofnx and ny. Detailed energy comparisons with the 2D QHO model are summarized in Table II. Minor devi- ations a...

  2. [2]

    T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. Monroe, and J. L. O’Brien, Quantum computers, Na- ture 464, 45 (2010)

  3. [3]

    Ebadi, T

    S. Ebadi, T. T. Wang, H. Levine, A. Keesling, G. Se- meghini, A. Omran, D. Bluvstein, R. Samajdar, H. Pich- ler, W. W. Ho, S. Choi, S. Sachdev, M. Greiner, V. Vuletić, and M. D. Lukin, Quantum phases of matter on a 256-atom programmable quantum simulator, Nature 595, 227 (2021)

  4. [4]

    Bluvstein, H

    D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner, V. Vuletić, and M. D. Lukin, A quantum processor based on coherent transport of en- tangled atom arrays, Nature604, 451 (2022)

  5. [5]

    C. D. Bruzewicz, J. Chiaverini, R. McConnell, and J. M. Sage, Trapped-ion quantum computing: Progress and challenges, Appl. Phys. Rev.6, 021314 (2019)

  6. [6]

    Monroe, W

    C. Monroe, W. C. Campbell, L.-M. Duan, Z.-X. Gong, A. V. Gorshkov, P. W. Hess, R. Islam, K. Kim, N. M. Linke, G. Pagano, P. Richerme, C. Senko, and N. Y. Yao, Programmable quantum simulations of spin systems with trapped ions, Rev. Mod. Phys.93, 025001 (2021)

  7. [7]

    N. P. de Leon, K. M. Itoh, D. Kim, K. K. Mehta, T. E. Northup, H. Paik, B. S. Palmer, N. Samarth, S. Sangtawesin, and D. W. Steuerman, Materials chal- lenges and opportunities for quantum computing hard- ware, Science372, eabb2823 (2021)

  8. [8]

    Acharya, D

    R. Acharya, D. A. Abanin, L. Aghababaie-Beni, I. Aleiner, others, G. Q. AI, and Collaborators, Quan- tum error correction below the surface code threshold, Nature 638, 920 (2025)

  9. [9]

    Pezzagna and J

    S. Pezzagna and J. Meijer, Quantum computer based on color centers in diamond, Appl. Phys. Rev. 8, 011308 (2021)

  10. [10]

    S. G. J. Philips, M. T. Mądzik, S. V. Amitonov, S. L. de Snoo, M. Russ, N. Kalhor, C. Volk, W. I. L. Lawrie, D. Brousse, L. Tryputen, B. P. Wuetz, A. Sammak, M. Veldhorst, G. Scappucci, and L. M. K. Vandersypen, Universal control of a six-qubit quantum processor in sil- icon, Nature609, 919 (2022)

  11. [11]

    Burkard, T

    G. Burkard, T. D. Ladd, A. Pan, J. M. Nichol, and J. R. Petta, Semiconductor spin qubits, Rev. Mod. Phys.95, 025003 (2023)

  12. [12]

    Borsoi, N

    F. Borsoi, N. W. Hendrickx, V. John, M. Meyer, S. Motz, F. van Riggelen, A. Sammak, S. L. de Snoo, G. Scappucci, and M. Veldhorst, Shared control of a 16semiconductor quantum dot crossbar array, Nat. Nan- otechnol. 19, 21 (2024)

  13. [13]

    Y.Pan, J.Yang, S.C.Erwin, K.Kanisawa,andS.Fölsch, Reconfigurable quantum-dot molecules created by atom manipulation, Phys. Rev. Lett.115, 076803 (2015)

  14. [14]

    Yan and P

    L. Yan and P. Liljeroth, Engineered electronic states in atomically precise artificial lattices and graphene nanoribbons, Adv. Phys.: X4, 1651672 (2019)

  15. [15]

    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é materials, Nat. Rev. Mater.6, 201 (2021)

  16. [16]

    D. M. Kennes, M. Claassen, L. Xian, A. Georges, A. J. Millis, J. Hone, C. R. Dean, D. N. Basov, A. N. Pa- supathy, and A. Rubio, Moiré heterostructures as a condensed-matterquantumsimulator,Nat.Phys. 17,155 (2021)

  17. [17]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Moiré bands in twisted double-layer graphene, Proc. Natl. Acad. Sci. U.S.A. 108, 12233 (2011)

  18. [18]

    Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, and P. Jarillo- Herrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 556, 80 (2018)

  19. [19]

    A. R.-P. Montblanch, M. Barbone, I. Aharonovich, M. Atatüre, and A. C. Ferrari, Layered materials as a platform for quantum technologies, Nat. Nanotechnol. 18, 555 (2023)

  20. [20]

    Trambly De Laissardière, D

    G. Trambly De Laissardière, D. Mayou, and L. Magaud, Localization of Dirac Electrons in Rotated Graphene Bi- layers, Nano Lett.10, 804 (2010)

  21. [21]

    J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, Continuum model of the twisted graphene bilayer, Phys. Rev. B86, 155449 (2012)

  22. [22]

    K. Wang, K. De Greve, L. A. Jauregui, A. Sushko, A. High, Y. Zhou, G. Scuri, T. Taniguchi, K. Watanabe, M. D. Lukin, H. Park, and P. Kim, Electrical control of charged carriers and excitons in atomically thin materi- als, Nat. Nanotechnol.13, 128 (2018)

  23. [23]

    M. H. Naik, S. Kundu, I. Maity, and M. Jain, Origin and evolution of ultraflat bands in twisted bilayer transition metal dichalcogenides: Realization of triangular quan- tum dots, Phys. Rev. B102, 075413 (2020)

  24. [24]

    V. V. Enaldiev, F. Ferreira, J. G. McHugh, and V. I. Fal’ko, Self-organized quantum dots in marginally twisted MoSe2/WSe2 and MoS2/WS2 bilayers, npj 2D Mater. Appl.6, 74 (2022)

  25. [25]

    F. Li, D. Lee, N. Leconte, S. Javvaji, Y. D. Kim, and J. Jung, Moiré flat bands and antiferroelectric domains in lattice relaxed twisted bilayer hexagonal boron nitride under perpendicular electric fields, Phys. Rev. B 110, 155419 (2024)

  26. [26]

    L. Xian, D. M. Kennes, N. Tancogne-Dejean, M. Altarelli, and A. Rubio, Multiflat bands and strong correlations in twisted bilayer boron nitride: Doping- induced correlated insulator and superconductor, Nano Lett. 19, 4934 (2019)

  27. [27]

    Nakatsuji, T

    N. Nakatsuji, T. Kawakami, H. Tateishi, K. Kato, and M. Koshino, Moiré band engineering in twisted trilayer WSe2 (2025), 2504.20449

  28. [28]

    Ribeiro-Palau, C

    R. Ribeiro-Palau, C. Zhang, K. Watanabe, T. Taniguchi, J. Hone, and C. R. Dean, Twistable electronics with dy- namically rotatable heterostructures, Science 361, 690 (2018). 19

  29. [29]

    H. Tang, Y. Wang, X. Ni, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero, S. Fan, E. Mazur, A. Yacoby, and Y. Cao, On-chip multi-degree-of-freedom control of two- dimensional materials, Nature632, 1038 (2024)

  30. [30]

    Nakatsuji, T

    N. Nakatsuji, T. Kawakami, and M. Koshino, Multiscale lattice relaxation in general twisted trilayer graphenes, Phys. Rev. X13, 041007 (2023)

  31. [31]

    I. M. Craig, M. Van Winkle, C. Groschner, K. Zhang, N. Dowlatshahi, Z. Zhu, T. Taniguchi, K. Watanabe, S. M. Griffin, and D. K. Bediako, Local atomic stacking and symmetry in twisted graphene trilayers, Nat. Mater. 23, 323 (2024)

  32. [32]

    Van Winkle, N

    M. Van Winkle, N. Dowlatshahi, N. Khaloo, M. Iyer, I. M. Craig, R. Dhall, T. Taniguchi, K. Watanabe, and D. K. Bediako, Engineering interfacial polarization switching in van der waals multilayers, Nat. Nanotechnol. 19, 751 (2024)

  33. [33]

    D. Park, C. Park, K. Yananose, E. Ko, B. Kim, R. En- gelke, X. Zhang, K. Davydov, M. Green, H.-M. Kim, S. H. Park, J. H. Lee, S.-G. Kim, H. Kim, K. Watan- abe, T. Taniguchi, S. M. Yang, K. Wang, P. Kim, Y.-W. Son, and H. Yoo, Unconventional domain tessellations in moiré-of-moiré lattices, Nature641, 896 (2025)

  34. [34]

    Vizner Stern, Y

    M. Vizner Stern, Y. Waschitz, W. Cao, I. Nevo, K. Watanabe, T. Taniguchi, E. Sela, M. Urbakh, O. Hod, andM.BenShalom,Interfacialferroelectricitybyvander Waals sliding, Science372, 1462 (2021)

  35. [35]

    Yasuda, X

    K. Yasuda, X. Wang, K. Watanabe, T. Taniguchi, and P. Jarillo-Herrero, Stacking-engineered ferroelectricity in bilayer boron nitride, Science372, 1458 (2021)

  36. [36]

    Yasuda, E

    K. Yasuda, E. Zalys-Geller, X. Wang, D. Ben- nett, S. S. Cheema, K. Watanabe, T. Taniguchi, E. Kaxiras, P. Jarillo-Herrero, and R. Ashoori, Ultrafast high-endurance memory based on sliding ferroelectrics, Science 385, 53 (2024), https://www.science.org/doi/pdf/10.1126/science.adp3575

  37. [37]

    R. Bian, R. He, E. Pan, Z. Li, G. Cao, P. Meng, J. Chen, Q. Liu, Z. Zhong, W. Li, and F. Liu, Developing fatigue-resistant ferroelectrics using in- terlayer sliding switching, Science 385, 57 (2024), https://www.science.org/doi/pdf/10.1126/science.ado1744

  38. [38]

    Hsieh, Y.-P

    C.-Y. Hsieh, Y.-P. Shim, M. Korkusinski, and P. Hawry- lak, Physics of lateral triple quantum-dot molecules with controlled electron numbers, Rep. Prog. Phys.75, 114501 (2012)

  39. [39]

    Constantinescu, A

    G. Constantinescu, A. Kuc, and T. Heine, Stacking in bulk and bilayer hexagonal boron nitride, Phys. Rev. Lett. 111, 036104 (2013)

  40. [40]

    J. H. Warner, M. H. Rümmeli, A. Bachmatiuk, and B. Büchner, Atomic resolution imaging and topography of boron nitride sheets produced by chemical exfoliation, ACS Nano4, 1299 (2010)

  41. [41]

    Park and Y.-W

    C. Park and Y.-W. Son, Condensation of preformed charge density waves in kagome metals, Nat. Commun. 14, 7309 (2023)

  42. [42]

    Baroni, S

    S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Gi- annozzi, Phonons and related crystal properties from density-functional perturbation theory, Rev. Mod. Phys. 73, 515 (2001)

  43. [43]

    A. N. Kolmogorov and V. H. Crespi, Registry-dependent interlayer potential for graphitic systems, Phys. Rev. B 71, 235415 (2005)

  44. [44]

    Lee and Y.-W

    S.-H. Lee and Y.-W. Son, First-principles approach with a pseudohybrid density functional for extended Hubbard interactions, Phys. Rev. Res.2, 043410 (2020)

  45. [45]

    Tancogne-Dejean and A

    N. Tancogne-Dejean and A. Rubio, Parameter-free hy- bridlike functional based on an extended Hubbard model: DFT +U +V, Phys. Rev. B102, 155117 (2020)

  46. [46]

    Timrov, F

    I. Timrov, F. Aquilante, L. Binci, M. Cococcioni, and N. Marzari, Pulay forces in density-functional theory with extended hubbard functionals: From nonorthogo- nalized to orthogonalized manifolds, Phys. Rev. B102, 235159 (2020)

  47. [47]

    Liu, N.-J

    W. Liu, N.-J. Guo, S. Yu, Y. Meng, Z.-P. Li, Y.-Z. Yang, Z.-A. Wang, X.-D. Zeng, L.-K. Xie, Q. Li, J.-F. Wang, J.- S. Xu, Y.-T. Wang, J.-S. Tang, C.-F. Li, and G.-C. Guo, Spin-active defects in hexagonal boron nitride, Mater. Quantum Technol.2, 032002 (2022)

  48. [48]

    Otani and S

    M. Otani and S. Okada, Gate-controlled carrier injection into hexagonal boron nitride, Phys. Rev. B83, 073405 (2011)

  49. [49]

    Pizzi, V

    G. Pizzi, V. Vitale, R. Arita, S. Blügel, F. Freimuth, G. Géranton, M. Gibertini, D. Gresch, C. Johnson, T. Koretsune, J. Ibañez-Azpiroz, H. Lee, J.-M. Lihm, D. Marchand, A. Marrazzo, Y. Mokrousov, J. I. Mustafa, Y. Nohara, Y. Nomura, L. Paulatto, S. Poncé, T. Pon- weiser, J. Qiao, F. Thöle, S. S. Tsirkin, M. Wierzbowska, N. Marzari, D. Vanderbilt, I. Sou...

  50. [50]

    T. Haga, Y. Matsuura, Y. Fujimoto, and S. Saito, Elec- tronic states and modulation doping of hexagonal boron nitride trilayers, Phys. Rev. Mater.5, 094003 (2021)

  51. [51]

    S.-Y.LeeandY.Chung,Chargestabilityofatriplequan- tum dot with a finite tunnel coupling, Phys. Rev. B87, 045302 (2013)

  52. [52]

    Loss and D

    D. Loss and D. P. DiVincenzo, Quantum computation with quantum dots, Phys. Rev. A57, 120 (1998)

  53. [53]

    M. Seo, H. K. Choi, S.-Y. Lee, N. Kim, Y. Chung, H.-S. Sim, V. Umansky, and D. Mahalu, Charge frustration in a triangular triple quantum dot, Phys. Rev. Lett.110, 046803 (2013)

  54. [54]

    A. Hamo, A. Benyamini, I. Shapir, I. Khivrich, J. Waiss- man, K. Kaasbjerg, Y. Oreg, F. von Oppen, and S. Ilani, Electron attraction mediated by coulomb repulsion, Na- ture 535, 395 (2016)

  55. [55]

    J. P. Dehollain, U. Mukhopadhyay, V. P. Michal, Y. Wang, B. Wunsch, C. Reichl, W. Wegscheider, M. S. Rudner, E. Demler, and L. M. K. Vandersypen, Nagaoka ferromagnetism observed in a quantum dot plaquette, Nature 579, 528 (2020)

  56. [56]

    R. Gong, X. Du, E. Janzen, V. Liu, Z. Liu, G. He, B. Ye, T. Li, N. Y. Yao, J. H. Edgar, E. A. Henriksen, and C. Zu, Isotope engineering for spin defects in van der waals materials, Nat. Commun.15, 104 (2024)

  57. [57]

    Rizzato, M

    R. Rizzato, M. Schalk, S. Mohr, J. C. Hermann, J. P. Leibold, F. Bruckmaier, G. Salvitti, C. Qian, P. Ji, G. V. Astakhov, U. Kentsch, M. Helm, A. V. Stier, J. J. Finley, and D. B. Bucher, Extending the coherence of spin de- fects in hbn enables advanced qubit control and quantum sensing, Nat. Commun.14, 5089 (2023)

  58. [58]

    T. T. Tran, K. Bray, M. J. Ford, M. Toth, and I. Aharonovich, Quantum emission from hexagonal boron nitride monolayers, Nat. Nanotechnol.11, 37 (2016)

  59. [59]

    Große, M

    J. Große, M. von Helversen, A. Koulas-Simos, M. Her- mann, and S. Reitzenstein, Development of site- controlled quantum dot arrays acting as scalable sources 20 of indistinguishable photons, APL Photonics5, 096107 (2020)

  60. [60]

    C. Su, F. Zhang, S. Kahn, B. Shevitski, J. Jiang, C. Dai, A. Ungar, J.-H. Park, K. Watanabe, T. Taniguchi, J. Kong, Z. Tang, W. Zhang, F. Wang, M. Crommie, S. G. Louie, S. Aloni, and A. Zettl, Tuning colour cen- tres at a twisted hexagonal boron nitride interface, Nat. Mater. 21, 896 (2022)

  61. [61]

    J. D. Caldwell, I. Aharonovich, G. Cassabois, J. H. Edgar, B. Gil, and D. N. Basov, Photonics with hexago- nal boron nitride, Nat. Rev. Mater.4, 552 (2019)

  62. [62]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

  63. [63]

    M. H. Naik, I. Maity, P. K. Maiti, and M. Jain, Kol- mogorov–crespi potential for multilayer transition-metal dichalcogenides: Capturing structural transformations in moiré superlattices, J. Phys. Chem. C123, 9770 (2019)

  64. [64]

    Peng, Z.-H

    H. Peng, Z.-H. Yang, J. P. Perdew, and J. Sun, Versa- tile van der waals density functional based on a meta- generalized gradient approximation, Phys. Rev. X 6, 041005 (2016)

  65. [65]

    E. A. Cortés, J. M. Florez, and E. S. Morell, Ferroelec- tric response to interlayer shifting and rotations in tri- layerhexagonalboronnitride,J.Phys.Chem.Solids 173, 111086 (2023)

  66. [66]

    Yang, S.-H

    W. Yang, S.-H. Jhi, S.-H. Lee, and Y.-W. Son, Ab ini- tio study of lattice dynamics of group IV semiconductors using pseudohybrid functionals for extended Hubbard in- teractions, Phys. Rev. B104, 104313 (2021)

  67. [67]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, A. Dal Corso, S. De Gironcoli, S. Fabris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, S. S...

  68. [68]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavaz- zoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carn- imeo, A. Dal Corso, S. De Gironcoli, P. Delugas, R. A. DiStasio, A. Ferretti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia, M. Kawamura, H.-Y. Ko, A. Ko...

  69. [69]

    Javvaji, F

    S. Javvaji, F. Li, and J. Jung,Ab initio tight-binding models for mono- and bilayer hexagonal boron nitride (h-BN), Phys. Rev. Materials9, 024004 (2025)

  70. [70]

    R. B. Lehoucq, D. C. Sorensen, and C. Yang,ARPACK users’ guide: solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods(Society for In- dustrial and Applied Mathematics, Philadelphia, 1998)

  71. [71]

    Opencollab, Arpack-ng, https://github.com/ opencollab/arpack-ng, [Accessed: 24-July-2025]

  72. [72]

    Balay, W

    S. Balay, W. D. Gropp, L. C. McInnes, and B. F. Smith, Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries, inModern Software Tools for Scientific Computing, edited by E. Arge, A. M. Bru- aset, and H. P. Langtangen (Birkhäuser, Boston, MA,