arxiv: 2605.03459 · v1 · submitted 2026-05-05 · ❄️ cond-mat.mes-hall · cond-mat.other

Recognition: unknown

Cubic edge dispersion in a semi-Dirac Chern insulator

Marta Garc\'ia Olmos , David Mart\'in Tejedor , Mario Amado , Yuriko Baba , Rafael A. Molina

Authors on Pith no claims yet

Pith reviewed 2026-05-07 15:15 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.other

keywords semi-Dirac Chern insulatorchiral edge statescubic dispersionmass domain walltwo-band lattice modeltopological phase diagramtight-binding ribbon

0 comments

The pith

In a semi-Dirac Chern insulator the chiral edge states disperse as the cube of momentum.

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

The paper introduces a minimal two-band lattice model that realizes a Chern insulator with semi-Dirac bulk bands and shows that its chiral edge states obey a cubic dispersion relation at low energy. This result is obtained by solving the edge problem analytically through a mass-domain-wall construction in a semi-infinite geometry, then verified by direct diagonalization of a ribbon geometry. A reader would care because conventional Chern insulators produce linear edge dispersion inherited from Dirac cones, whereas the anisotropic semi-Dirac case produces a qualitatively different boundary spectrum that can alter transport and response properties.

Core claim

By constructing a minimal two-band lattice model for a semi-Dirac Chern insulator and applying a mass-domain-wall construction in semi-infinite geometry, the authors derive an explicit expression for the chiral edge states whose low-energy dispersion scales as E(k) ∝ k³, a result confirmed by numerical tight-binding calculations on a ribbon.

What carries the argument

The mass-domain-wall construction applied to the two-band semi-Dirac Hamiltonian, which produces an analytic solution for the edge-state wave function and dispersion.

If this is right

The linear dispersion relation of edge states is not universal in Chern insulators once the bulk bands become anisotropic semi-Dirac cones.
Cubic scaling changes the group velocity and density of states of the chiral modes relative to the linear case.
The same mass-domain-wall technique can be used to engineer other nonstandard edge dispersions by tuning bulk band anisotropy.
Topological phase diagrams of semi-Dirac models can be determined analytically before checking edge properties.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Cubic edge dispersion may produce distinct signatures in transport or Andreev reflection experiments compared with linear modes.
The result suggests that other anisotropic bulk dispersions could generate higher-order edge scalings that remain unexplored.
Realization in strained lattices or cold-atom systems would allow direct measurement of the k³ velocity scaling.

Load-bearing premise

The minimal two-band lattice model realizes a semi-Dirac Chern insulator whose topological properties and edge states are faithfully captured by the mass-domain-wall construction without higher-order corrections or lattice artifacts.

What would settle it

Numerical diagonalization of the tight-binding ribbon that shows the low-momentum edge dispersion deviating from cubic scaling would falsify the analytic prediction.

Figures

Figures reproduced from arXiv: 2605.03459 by David Mart\'in Tejedor, Mario Amado, Marta Garc\'ia Olmos, Rafael A. Molina, Yuriko Baba.

**Figure 1.** Figure 1: FIG. 1. (a) Bulk spectrum, (b) energy view at source ↗

**Figure 2.** Figure 2: FIG. 2. (a) Minimum energy gap view at source ↗

**Figure 3.** Figure 3: FIG. 3. Band dispersion and Zak phase for semi-infinite view at source ↗

**Figure 5.** Figure 5: FIG. 5. (a) Dispersion relation and (b) probability density of view at source ↗

**Figure 4.** Figure 4: FIG. 4. (a) Dispersion relation and (b) probability density of view at source ↗

read the original abstract

Topological edge states in Chern insulators are typically characterized by a linear dispersion relation inherited from the Dirac structure of the bulk Hamiltonian. Here we show that this paradigm can be fundamentally altered in systems with anisotropic semi-Dirac band structures. We introduce a minimal two-band lattice model realizing a semi-Dirac Chern insulator and determine its topological phase diagram analytically. Using a mass-domain-wall approach in a semi-infinite geometry, we derive an explicit expression for the chiral edge states and find that their low-energy dispersion scales cubically with momentum, $E(k)\propto k^3$. Numerical diagonalization of the corresponding tight-binding ribbon confirms the analytical prediction. Our results demonstrate that unconventional bulk band structures can produce qualitatively different boundary excitations, providing a route to engineering nonstandard chiral edge dynamics in topological materials and synthetic quantum systems.

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.

Desk Editor's Note private letter to a colleague

A minimal model shows cubic edge dispersion from semi-Dirac bulk in a Chern insulator.

read the letter

The key point from this paper is that they have a semi-Dirac Chern insulator where the edge states show cubic dispersion at low momentum, E(k) proportional to k cubed. They do a good job constructing a minimal two-band lattice model that achieves a semi-Dirac point while maintaining Chern insulator topology. They analytically determine the topological phase diagram and then use the mass-domain-wall approach in semi-infinite geometry to derive the explicit edge state wavefunction and its cubic dispersion. The numerical diagonalization of the tight-binding ribbon backs this up independently. That combination of analytic derivation and numerical check is what makes the result convincing. The potential soft spots are that this is a tailored minimal model, so the cubic behavior could be specific to their choice of hopping and mass terms rather than universal for semi-Dirac systems. The abstract doesn't include the Hamiltonian itself, which leaves some questions about the exact parameter values or possible lattice artifacts. Still, the mutual support between the two methods suggests the scaling is robust within the model. This work is for condensed matter physicists focused on topological edge states and how bulk band features influence them. Readers who work with lattice models or are looking for ways to get non-linear chiral modes will find it relevant. It deserves serious peer review because the claim is well-supported by the calculations described and opens a clear direction for further study. Recommendation: Yes, it should go to referees.

Referee Report

0 major / 3 minor

Summary. The manuscript introduces a minimal two-band lattice model realizing a semi-Dirac Chern insulator, analytically determines its topological phase diagram, and employs a mass-domain-wall construction in semi-infinite geometry to derive an explicit form for the chiral edge states whose low-energy dispersion scales as E(k) ∝ k³. Independent numerical diagonalization of the corresponding tight-binding ribbon is reported to confirm the analytic prediction.

Significance. If the central result holds, the work demonstrates that anisotropic semi-Dirac bulk dispersions can produce qualitatively non-linear (cubic) chiral edge modes, departing from the linear Dirac paradigm standard in Chern insulators. The combination of an explicit analytic edge-state wavefunction with direct lattice confirmation is a strength, as is the parameter-free character of the derived cubic scaling once the model is fixed. This opens a concrete route to engineering unconventional edge dynamics in both solid-state and synthetic topological platforms.

minor comments (3)

§2 (Model Hamiltonian): the explicit 2×2 lattice Hamiltonian and the precise values of the hopping amplitudes and mass term used for the phase diagram and edge-state calculation should be stated in a single equation block for reproducibility; the current description leaves the reader to infer the precise form from the abstract and later sections.
§4 (Edge-state derivation): the mass-domain-wall ansatz is applied to a semi-infinite geometry; a brief remark on the validity of neglecting higher-order lattice corrections or finite-size effects at the energies where the k³ scaling is claimed would strengthen the analytic-numerical comparison.
Figure 3 (Ribbon spectrum): the momentum range shown for the edge branch should be labeled with the same units and cutoff as the analytic E(k)∝k³ curve to facilitate direct visual comparison; the inset zoom is helpful but the scaling collapse is not shown.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, the clear summary of its contributions, and the recommendation for minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper defines a minimal two-band lattice model, analytically obtains its topological phase diagram, derives the explicit chiral edge-state wavefunction and cubic dispersion E(k)∝k³ via the standard mass-domain-wall construction in semi-infinite geometry, and verifies the result by independent numerical diagonalization of the tight-binding ribbon. All load-bearing steps start from the explicitly constructed Hamiltonian and apply standard analytic/numeric techniques without reducing to self-defined fitted quantities, self-citation chains, or ansatzes smuggled from prior work. The cubic scaling is an output of the derivation, not an input.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on a constructed two-band lattice model whose parameters are chosen to produce semi-Dirac dispersion and nonzero Chern number; the mass-domain-wall method is applied as a standard technique.

free parameters (2)

hopping amplitudes
Tunable parameters in the lattice model required to realize the anisotropic semi-Dirac bulk bands and topological phase.
mass term strength
Parameter controlling the domain wall that defines the edge in the semi-infinite geometry.

axioms (1)

domain assumption The two-band lattice model captures the essential low-energy physics of a semi-Dirac Chern insulator.
Minimal-model assumption standard in condensed-matter theory but not independently verified here.

pith-pipeline@v0.9.0 · 5452 in / 1316 out tokens · 66193 ms · 2026-05-07T15:15:48.926923+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

40 extracted references

[1]

J. K. Asbóth, L. Oroszlány, and A. Pályi,A Short Course on Topological Insulators(Springer, 2016)

2016
[2]

B. A. Bernevig and T. L. Hughes,Topological Insulators and Topological Superconductors(Princeton University Press, Princeton, NJ, 2013)

2013
[3]

parity anomaly

F. D. M. Haldane, Model for a quantum hall effect without landau levels: Condensed-matter realization of the "parity anomaly", Phys. Rev. Lett.61, 2015 (1988)

2015
[4]

Hatsugai, Edge states in the integer quantum hall effect and the riemann surface of the bloch function, Phys

Y. Hatsugai, Edge states in the integer quantum hall effect and the riemann surface of the bloch function, Phys. Rev. B48, 11851 (1993)

1993
[5]

M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys.82, 3045 (2010)

2010
[6]

Qi and S.-C

X.-L. Qi and S.-C. Zhang, Topological insulators and superconductors, Rev. Mod. Phys.83, 1057 (2011)

2011
[7]

Anirban, 15 years of topological insulators, Nature Reviews Physics5, 267 (2023)

A. Anirban, 15 years of topological insulators, Nature Reviews Physics5, 267 (2023)

2023
[8]

X. L. Qi, Y. S. Wu, and S. C. Zhang, Topological quantization of the spin hall effect in two-dimensional paramagnetic semiconductors, Phys. Rev. B74, 085308 (2006)

2006
[9]

B. A. Bernevig, T. L. Hughes, and S. C. Zhang, Quantum spin hall effect and topological phase transition in hgte quantum wells, Science314, 1757 (2006)

2006
[10]

Katayama, A

S. Katayama, A. Kobayashi, and Y. Suzumura, Electric conductivity of the zero-gap semiconducting state inα- (BEDT-TTF)2I3 salt, J. Phys. Soc. Jap.75, 023708 (2006)

2006
[11]

Dietl, F

P. Dietl, F. Piéchon, and G. Montambaux, New magnetic field dependence of landau levels in a graphenelike structure, Phys. Rev. Lett.100, 236405 (2008)

2008
[12]

M. O. Goerbig, J.-N. Fuchs, G. Montambaux, and F. Piéchon, Tilted anisotropic dirac cones in quinoid- typegrapheneandα−(BEDT-TTF) 2i3,Phys.Rev.B78, 045415 (2008)

2008
[13]

Pardo and W

V. Pardo and W. E. Pickett, Half-metallic semi-dirac- point generated by quantum confinement intio 2/vo2 nanostructures, Phys. Rev. Lett.102, 166803 (2009)

2009
[14]

Banerjee, R

S. Banerjee, R. R. P. Singh, V. Pardo, and W. E. Pickett, Tight-binding modeling and low-energy behavior of the semi-dirac point, Phys. Rev. Lett.103, 016402 (2009)

2009
[15]

Zhong, Y

C. Zhong, Y. Chen, Y. Xie, Y.-Y. Sun, and S.Zhang,Semi-diracsemimetalinsiliceneoxide,Physical Chemistry Chemical Physics19, 3820 (2017)

2017
[16]

Y. Shao, S. Moon, A. N. Rudenko, J. Wang, J. Herzog- Arbeitman, M.Ozerov, D.Graf, Z.Sun, R.Queiroz, S.H. Lee, Y. Zhu, Z. Mao, M. I. Katsnelson, B. A. Bernevig, D. Smirnov, A. J. Millis, and D. N. Basov, Semi-dirac fermions in a topological metal, Phys. Rev. X14, 041057 (2024)

2024
[17]

H.Huang, Z.Liu, H.Zhang, W.Duan,andD.Vanderbilt, Emergence of a chern-insulating state from a semi-dirac dispersion, Phys. Rev. B92, 161115 (2015)

2015
[18]

García Olmos, Y

M. García Olmos, Y. Baba, M. Amado, and R. A. Molina, Zero momentum topological insulator in 2d semi- dirac materials, Journal of Physics: Materials7, 045008 (2024)

2024
[19]

García Olmos, Y

M. García Olmos, Y. Baba, A. López, M. Amado, and R. A. Molina, Spin-dependent transport through edge states in 2d semi-dirac materials with rashba spin-orbit coupling and band inversion, 2D Materials12, 045019 (2025)

2025
[20]

Zak, Berry’s phase for energy bands in solids, Phys

J. Zak, Berry’s phase for energy bands in solids, Phys. Rev. Lett.62, 2747 (1989)

1989
[21]

W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in polyacetylene, Phys. Rev. Lett.42, 1698 (1979)

1979
[22]

C.-K. Chiu, J. C. Teo, A. P. Schnyder, and S. Ryu, Classification of topological quantum matter with symmetries, Reviews of Modern Physics88, 035005 (2016). 6

2016
[23]

S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, Topological insulators and superconductors: tenfold way and dimensional hierarchy, New Journal of Physics12, 065010 (2010)

2010
[24]

A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Classification of topological insulators and superconductors in three spatial dimensions, Physical Review B78, 195125 (2008)

2008
[25]

[31, 32], where we provide further details about the numerical calculation of the topological invariants and the analytical edge state dispersion derivation

See Supplemental Material, including Refs. [31, 32], where we provide further details about the numerical calculation of the topological invariants and the analytical edge state dispersion derivation
[26]

M. C. Rechtsman, J. M. Zeuner, Y. Plotnik,et al., Photonic floquet topological insulators, Nature496, 196 (2013)

2013
[27]

L. Lu, J. D. Joannopoulos, and M. Soljačić, Topological photonics, Nature Photonics8, 821 (2014)

2014
[28]

Goldman, J

N. Goldman, J. C. Budich, and P. Zoller, Topological quantum matter with ultracold gases in optical lattices, Nature Physics12, 639 (2016)

2016
[29]

Aidelsburger, M

M. Aidelsburger, M. Lohse, and C. Schweizer, Artificial gauge fields with ultracold atoms in optical lattices, Nature Physics14, 100 (2018)

2018
[30]

Tarruell, D

L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, and T. Esslinger, Creating, moving and merging dirac points with a fermi gas in a tunable honeycomb lattice, Nature 483, 302 (2012)

2012
[31]

Fukui, Y

T. Fukui, Y. Hatsugai, and H. Suzuki, Chern numbers in discretized brillouin zone: Efficient method of computing (spin) hall conductances, Journal of the Physical Society of Japan74, 1674 (2005)

2005
[32]

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

D. Vanderbilt,Berry Phases in Electronic Structure Theory: Electric Polarization, Orbital Magnetization and Topological Insulators(Cambridge University Press, 2018). END MATTER Analytical Berry phase—The Berry phase of the model can be written as Ω−(k) = αBW |E−(k)|3 Bω even(k) +vω odd(k) (19a) withω odd(even)(k)is an odd (even) function ofky: ωeven = (1−...

2018
[33]

the boundary problem forH0 supports a localized bound stateψ 0
[34]

the operatorV 1 has a nonzero matrix element in the bound-state subspace: ⟨ψ0|V1|ψ0⟩ ̸= 0.(25a) Then, the energy of the in-gap state is expanded as E(k∥) =v edgek∥ +· · ·, v edge =⟨ψ 0|V1|ψ0⟩.(25b) By contrast, acubicedge mode arises when the following three conditions hold: 7
[35]

the edge Hamiltonian obeys the symmetry ΓH(k ∥)Γ =−H(−k ∥),(26a) for some unitary operatorΓ, so thatE(k ∥) = −E(−k∥)and all even powers ofk∥ are forbidden
[36]

the linear coefficient vanishes, ⟨ψ0|V1|ψ0⟩= 0,(26b) which happens, for example, whenV1 is derivative- like and integration by parts kills its expectation value on a normalizable bound state
[37]

Under these conditions, E(k∥) =c 3k3 ∥ +O(k 5 ∥).(26c) In the present semi-Dirac model (1), Eq

no additional symmetry forces the cubic coefficient to vanish. Under these conditions, E(k∥) =c 3k3 ∥ +O(k 5 ∥).(26c) In the present semi-Dirac model (1), Eq. (24) reads as H0 =−iv∂ y σx + mso +B W ∂2 y σz,(27a) V1 =−iα∂ y σy,(27b) V2 =B σ x −B W σz ,(27c) forx ∥ =xandx ⊥ =y. Notice that the anisotropic bulk Hamiltonian contains no independent linear tang...
[38]

(26a) is given byΓ =σy and it forces the edge dispersion to be odd inkx

the symmetry in Eq. (26a) is given byΓ =σy and it forces the edge dispersion to be odd inkx
[39]

The zero mode can be chosen as an eigenstate ofσy, σyψ0 =sψ 0, s=±1.(28) This way, the linear coefficient in Eq

letψ 0 be the zero-energy bound state ofH0. The zero mode can be chosen as an eigenstate ofσy, σyψ0 =sψ 0, s=±1.(28) This way, the linear coefficient in Eq. (26b) is ⟨ψ0|V1|ψ0⟩=−iαs⟨ψ 0|∂y|ψ0⟩=−iαs Z dy f(y)f ′(y), wheref(y)isthescalarenvelopeoftheboundstate. Since the state is normalizable, Z dy f(y)f ′(y) = 1 2 Z dy ∂y f(y) 2 = 0,(29) the linear coeffic...
[40]

CUBIC EDGE DISPERSION IN A SEMI-DIRAC CHERN INSULATOR

The exact calculation of the edge dispersion (17) shows that it has the form of Eq. (26c) with leading cubic term c3 =−η αBsgn(B W ) |v| .(30) For comparison, consider a standard Dirac-type Chern insulator near a gap closing, HDirac(k∥, k⊥) =v ∥k∥Σ∥ +v ⊥k⊥Σ⊥ +mΣ m .(31) Upon imposing an edge normal tok ⊥, one replaces k⊥ → −i∂ ⊥ and obtains a Jackiw–Rebbi...