arxiv: 2605.08236 · v1 · submitted 2026-05-07 · ❄️ cond-mat.mes-hall · quant-ph

Recognition: 2 theorem links

· Lean Theorem

Dynamically Characterizing the Structures of Dirac Points via Wave Packets

Dan-Dan Liang, Xin Shen, Zhi Li

Authors on Pith no claims yet

Pith reviewed 2026-05-12 01:33 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall quant-ph

keywords Dirac pointswave packet dynamicswinding numberZitterbewegungtopological band structuresgraphenehybrid pointsparabolic points

0 comments

The pith

Wave packets' center-of-mass motion and spin texture determine the winding numbers of Dirac and parabolic points.

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

The paper examines how wave packet dynamics can characterize the topological structures around Dirac points in a modified graphene lattice. By including third-nearest-neighbor couplings, the model generates additional Dirac points whose creation and annihilation produce hybrid and parabolic points. The center-of-mass trajectory of a wave packet encodes the winding number at Dirac points, while the spin texture encodes it at parabolic points. For gapped hybrid points the packet undergoes one-dimensional Zitterbewegung. This dynamical approach could provide an experimental route to probe topological features in materials without direct band mapping.

Core claim

In the tight-binding model for graphene augmented with third-nearest-neighbor hoppings, Dirac points appear in pairs that can merge into hybrid or parabolic points. The time evolution of an initial wave packet shows that its center-of-mass displacement directly reflects the winding number of a nearby Dirac point, and its spin polarization texture reflects the winding number of a parabolic point. Gapped hybrid points induce purely one-dimensional oscillatory motion known as Zitterbewegung.

What carries the argument

The center-of-mass motion and spin texture arising from the semiclassical or exact dynamics of wave packets in the extended tight-binding Hamiltonian.

If this is right

The winding number around each Dirac point fixes the direction and magnitude of the long-time center-of-mass drift of a wave packet prepared with appropriate momentum.
Parabolic points produce a characteristic precession or texture in the spin degree of freedom of the wave packet that reveals their topological charge.
A gapped hybrid point confines the wave-packet oscillation to a single spatial direction, producing a distinctive Zitterbewegung signature.
Tracking these dynamical signatures allows identification of the points at which Dirac pairs are created or annihilated during parameter tuning.

Where Pith is reading between the lines

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

If realized in cold-atom or photonic lattices, this method could map topological transitions in real time by watching packet trajectories.
The approach might extend to three-dimensional Weyl points by generalizing the wave-packet observables.
Disorder or interactions would likely blur the clean signatures, implying the need for sufficiently coherent samples.

Load-bearing premise

The ideal tight-binding Hamiltonian with third-nearest-neighbor terms together with semiclassical or exact wave-packet evolution faithfully reproduces the topological winding numbers without disorder, interactions, or finite-size effects changing the observed dynamics.

What would settle it

Launching a wave packet near a Dirac point whose winding number is independently known and finding that its center-of-mass motion does not match the predicted drift would falsify the claimed correspondence.

Figures

Figures reproduced from arXiv: 2605.08236 by Dan-Dan Liang, Xin Shen, Zhi Li.

**Figure 2.** Figure 2: Band structures (left panel) of the Hamiltonian Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** Figure 3: (Color online) Analytically (lines) and numerically [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

**Figure 4.** Figure 4: (Color online) Analytically (lines) and numerically [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

read the original abstract

Topological non-trivial band structures are the core problem in the field of topological materials. In this paper, we investigate the topological band structure in a system with controllable Dirac points from the perspective of wave packet dynamics. By adding a third-nearest-neighboring coupling to the graphene model, additional pairs of Dirac points emerge. The emergence and annihilation of Dirac points result in hybrid and parabolic points, and we show that these band structures can be revealed by the dynamical behaviors of wave packets. Particularly, for the gapped hybrid point, the motion of the wave packet shows a one-dimensional \emph{Zitterbewegung} motion. Furthermore, we also show that the winding number associated with the Dirac point and parabolic point can be determined via the center-of-mass and spin texture of wave packets, respectively. The results of this work could motivate new experimental methods to characterize the system's topological signatures through wave packet dynamics, which may also find application in systems of other exotic topological materials.

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

The paper uses wave-packet dynamics to characterize Dirac, hybrid, and parabolic points in a modified graphene lattice, but the link from observables to winding numbers remains more suggestive than rigorously derived.

read the letter

Your colleague should know that this work takes a standard graphene tight-binding model, adds third-nearest-neighbor terms to split and create hybrid and parabolic points, and then runs wave-packet simulations to look at center-of-mass trajectories and spin textures. The new part is applying these dynamical probes specifically to the hybrid and parabolic cases that appear when the extra hopping is tuned. They report that gapped hybrids show one-dimensional Zitterbewegung, and they claim the winding numbers at the other points can be read off from the center-of-mass motion for Dirac points and from the spin texture for parabolic points. This is a reasonable extension of existing wave-packet techniques in topological materials. The numerics are the usual kind for these models, and the qualitative distinction in dynamics is plausible given the different dispersions. The main limitation is that the abstract and the stress-test both indicate no explicit formula is given that turns, say, an oscillation amplitude or a texture winding into the integer winding number. Without that step, or without checks that the mapping survives changes in initial wave-packet width or small disorder, the determination claim is hard to evaluate. The paper also assumes the ideal lattice without interactions or finite-size effects, which is common but limits how directly it applies to real samples. This is the kind of paper that could interest experimental groups working on graphene or similar 2D systems who want dynamical ways to map band features. It is not a major theoretical advance, but the concrete demonstration is worth a look for specialists. I would recommend sending it for peer review. The calculations are accessible and the topic is timely, so referees can check the details and ask for the missing mapping if it is not there.

Referee Report

2 major / 2 minor

Summary. The manuscript studies wave-packet dynamics in a tight-binding graphene model augmented by third-nearest-neighbor hopping. Adding this term produces additional Dirac points whose pairwise emergence and annihilation create hybrid and parabolic points. The authors report that these structures are visible in the time evolution of wave packets, with gapped hybrid points exhibiting one-dimensional Zitterbewegung, and claim that the winding numbers of Dirac points and parabolic points are quantitatively recoverable from the center-of-mass trajectory and the spin texture, respectively.

Significance. If an explicit, invertible mapping from the cited observables to the integer winding numbers can be established and shown to be robust, the work would supply a dynamical route to extracting topological invariants that complements static probes such as ARPES or transport. The approach could be relevant to both real materials and engineered lattices where wave-packet control is feasible.

major comments (2)

[Abstract] Abstract: the assertion that 'the winding number associated with the Dirac point and parabolic point can be determined via the center-of-mass and spin texture of wave packets' is not accompanied by any explicit formula, derivation, or numerical protocol that converts the observable (e.g., COM oscillation amplitude, frequency, or spin-texture winding) into the integer winding number. Without this mapping the claim reduces to a qualitative observation rather than a quantitative determination, which is load-bearing for the central result.
[Model and Results] Model definition and results sections: the third-nearest-neighbor hopping strength appears as a free parameter that controls the creation of hybrid/parabolic points, yet no systematic scan or error analysis is provided showing that the reported COM or spin-texture signatures remain one-to-one with the winding number when this parameter, initial-state choice, or weak disorder is varied. This leaves the robustness of the extraction unverified.

minor comments (2)

[Introduction] Notation for the winding number and for the spin texture is introduced without a clear reference to the standard definition used in the tight-binding literature; a brief equation or citation would improve clarity.
[Figures] Figure captions and axis labels for the wave-packet trajectories should explicitly state the initial conditions and the value of the third-nearest-neighbor parameter used in each panel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important points regarding the quantitative nature of our claims and the need to demonstrate robustness. We address each major comment below and have prepared revisions to strengthen the manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: the assertion that 'the winding number associated with the Dirac point and parabolic point can be determined via the center-of-mass and spin texture of wave packets' is not accompanied by any explicit formula, derivation, or numerical protocol that converts the observable (e.g., COM oscillation amplitude, frequency, or spin-texture winding) into the integer winding number. Without this mapping the claim reduces to a qualitative observation rather than a quantitative determination, which is load-bearing for the central result.

Authors: We agree that the abstract claim would benefit from greater precision. The manuscript presents numerical results in which the COM oscillation amplitude for Dirac points and the winding of the spin texture for parabolic points match the known integer winding numbers in the chosen parameter regimes. However, we did not supply an explicit analytic mapping or a reproducible numerical protocol that converts the raw observables into the integer. In the revised version we will add a dedicated subsection deriving the relation from the semiclassical wave-packet equations of motion (using the effective low-energy Hamiltonian around each point) and will include a clear step-by-step protocol for extracting the winding number from simulated or measured COM trajectories and spin textures. revision: yes
Referee: [Model and Results] Model definition and results sections: the third-nearest-neighbor hopping strength appears as a free parameter that controls the creation of hybrid/parabolic points, yet no systematic scan or error analysis is provided showing that the reported COM or spin-texture signatures remain one-to-one with the winding number when this parameter, initial-state choice, or weak disorder is varied. This leaves the robustness of the extraction unverified.

Authors: We acknowledge that the present manuscript illustrates the signatures only for representative values of the third-nearest-neighbor hopping and for a limited set of initial wave-packet parameters. To verify robustness we will add supplementary material containing (i) a systematic scan of the hopping strength across the emergence and annihilation points, (ii) results for several initial-state widths and momenta, and (iii) simulations with weak on-site disorder. These additions will demonstrate that the extracted winding numbers remain stable within the reported precision. revision: yes

Circularity Check

0 steps flagged

Wave-packet dynamics derivation is self-contained; no reduction to inputs by construction

full rationale

The paper starts from an explicit tight-binding Hamiltonian (graphene plus third-nearest-neighbor terms), evolves wave packets either semiclassically or exactly, and extracts center-of-mass trajectories and spin textures as observables. These observables are shown to encode the winding numbers of Dirac and parabolic points. No step equates a fitted parameter to a prediction, renames a known result, or relies on a self-citation chain whose uniqueness theorem is unverified outside the present work. The mapping from dynamics to winding number is derived from the equations of motion rather than presupposed in the Hamiltonian definition, satisfying the requirement that the central claim retain independent content.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The work rests on the standard nearest-neighbor graphene tight-binding Hamiltonian plus an added third-nearest-neighbor term whose strength is treated as a tunable parameter; no new particles or forces are postulated.

free parameters (1)

third-nearest-neighbor hopping strength
Tuned to control the emergence and annihilation of Dirac-point pairs; its specific value determines the locations of hybrid and parabolic points.

axioms (2)

domain assumption Electrons in the lattice are described by a non-interacting tight-binding Hamiltonian on a honeycomb lattice.
Invoked throughout the model construction and wave-packet evolution.
domain assumption Wave-packet dynamics can be obtained from the time-dependent Schrödinger equation on the lattice without decoherence or many-body effects.
Underlying all dynamical signatures reported.

pith-pipeline@v0.9.0 · 5466 in / 1328 out tokens · 41550 ms · 2026-05-12T01:33:42.182761+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
the winding number associated with the Dirac point and parabolic point can be determined via the center-of-mass and spin texture of wave packets, respectively
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
By adding a third-nearest-neighboring coupling to the graphene model, additional pairs of Dirac points emerge

Reference graph

Works this paper leans on

50 extracted references · 50 canonical work pages

[1]

N. W. Ashcroft and N. D. Mermin,Solid State Physics (Saunders College Publishing, 1976)

work page 1976
[2]

Wen, Rev

X.-G. Wen, Rev. Mod. Phys.89, 041004 (2017)

work page 2017
[3]

Bansil, H

A. Bansil, H. Lin, and T. Das, Rev. Mod. Phys.88, 021004 (2016)

work page 2016
[4]

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

work page 2010
[5]

Qi and S.-C

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

work page 2011
[6]

S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. Mele, and A. M. Rappe, Phys. Rev. Lett.108, 140405 (2012)

work page 2012
[7]

S. M. Young and C. L. Kane, Phys. Rev. Lett.115, 126803 (2015)

work page 2015
[8]

X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B83, 205101 (2011)

work page 2011
[9]

Lu and S.-Q

H.-Z. Lu and S.-Q. Shen, Frontiers of Physics12, 127201 (2017)

work page 2017
[10]

Z. Yang, F. Gao, X. Shi, X. Lin, Z. Gao, Y. Chong, and B. Zhang, Phys. Rev. Lett.114, 114301 (2015)

work page 2015
[11]

Peri, Z.-D

V. Peri, Z.-D. Song, M. Serra-Garcia, P. Engeler, R. Queiroz, X. Huang, W. Deng, Z. Liu, B. A. Bernevig, and S. D. Huber, Science367, 797 (2020)

work page 2020
[12]

Ozawa, H

T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zil- berberg, and I. Carusotto, Rev. Mod. Phys.91, 015006 (2019)

work page 2019
[13]

Zhang, Y.-Q

D.-W. Zhang, Y.-Q. Zhu, Y. X. Zhao, H. Yan, and S.-L. Zhu, Adv. Phys.67, 253 (2018)

work page 2018
[14]

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, Rev. Mod. Phys.93, 025001 (2021)

work page 2021
[15]

Nayak, S

C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys.80, 1083 (2008)

work page 2008
[16]

Shen and Z

X. Shen and Z. Li, Phys. Rev. A97, 013608 (2018). 5

work page 2018
[17]

Shen, Y.-Q

X. Shen, Y.-Q. Zhu, and Z. Li, Phys. Rev. B106, L180301 (2022)

work page 2022
[18]

C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Rev. Mod. Phys.88, 035005 (2016)

work page 2016
[19]

Tarruell, D

L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, and T. Esslinger, Nature483, 302 (2012)

work page 2012
[20]

L. Duca, T. Li, M. Reitter, I. Bloch, M. Schleier-Smith, and U. Schneider, Science347, 288 (2015)

work page 2015
[21]

Nielsen and M

H. Nielsen and M. Ninomiya, Phys. Lett. B130, 389 (1983)

work page 1983
[22]

Bena and L

C. Bena and L. Simon, Phys. Rev. B83, 115404 (2011)

work page 2011
[23]

Montambaux, Eur

G. Montambaux, Eur. Phy. J. B85, 375 (2012)

work page 2012
[24]

Montambaux, L.-K

G. Montambaux, L.-K. Lim, J.-N. Fuchs, and F. Pié- chon, Phys. Rev. Lett.121, 256402 (2018)

work page 2018
[25]

Schrödinger, Sitzungsber

E. Schrödinger, Sitzungsber. Preuss. Akad. Wiss. Phys. Math. Kl.24, 418 (1930)

work page 1930
[26]

N. R. Cooper, J. Dalibard, and I. B. Spielman, Rev. Mod. Phys.91, 015005 (2019)

work page 2019
[27]

Zhang, J.-B

Q. Zhang, J.-B. Gong, and C.-H. Oh, Chin. Phys. Lett. 30, 080301 (2013)

work page 2013
[28]

Deng, S.-T

D.-L. Deng, S.-T. Wang, K. Sun, and L.-M. Duan, Chin. Phys. Lett.35, 013701 (2017)

work page 2017
[29]

Guo, X.-X

G.-F. Guo, X.-X. Bao, L. Tan, and H.-Q. Gu, Chin. Phys. Lett.38, 040302 (2021)

work page 2021
[30]

Z. Li, H. Cao, and L.-B. Fu, Phys. Rev. A91, 023623 (2015)

work page 2015
[31]

Li, H.-Q

Z. Li, H.-Q. Wang, D.-W. Zhang, S.-L. Zhu, and D.-Y. Xing, Phys. Rev. A94, 043617 (2016)

work page 2016
[32]

Shen, D.-W

X. Shen, D.-W. Zhang, H. Yan, Z. Li, and S.-L. Zhu, Phys. Rev. Res.2, 013037 (2020)

work page 2020
[33]

Merkl, F

M. Merkl, F. E. Zimmer, G. Juzeli¯ unas, and P. öhberg, Europhys. Lett.83, 54002 (2008)

work page 2008
[34]

Zawadzki and T

W. Zawadzki and T. M. Rusin, J. Phys.: Condens. Mat- ter23, 143201 (2011)

work page 2011
[35]

Z. Meng, L. Wang, W. Han, F. Liu, K. Wen, C. Gao, P.Wang, C.Chin, andJ.Zhang,Nature615,231(2023)

work page 2023
[36]

Fallani, L

L. Fallani, L. D. Sarlo, J. E. Lye, M. Modugno, R. Saers, C. Fort, and M. Inguscio, Phys. Rev. Lett.93, 140406 (2004)

work page 2004
[37]

Kling, T

S. Kling, T. Salger, C. Grossert, and M. Weitz, Phys. Rev. Lett.105, 215301 (2010)

work page 2010
[38]

L. J. LeBlanc, M. C. Beeler, K. Jiménez-Garcĺa, A. R. Perry, S. Sugawa, R. A. Williams, and I. B. Spielman, New J. Phys.15, 073011 (2013)

work page 2013
[39]

E. Alba, X. Fernandez-Gonzalvo, J. Mur-Petit, J. K. Pa- chos, and J. J. Garcia-Ripoll, Phys. Rev. Lett.107, 235301 (2011)

work page 2011
[40]

Z. Wu, L. Zhang, W. Sun, X.-T. Xu, B.-Z. Wang, S.-C. Ji, Y. Deng, S. Chen, X.-J. Liu, and J.-W. Pan, Science 354, 83 (2016)

work page 2016
[41]

Hasan, C

M. Hasan, C. S. Madasu, K. D. Rathod, C. C. Kwong, C. Miniatura, F. Chevy, and D. Wilkowski, Phys. Rev. Lett.129, 130402 (2022)

work page 2022
[42]

X. Shen, F. Wang, Z. Li, and Z. Wu, Phys. Rev. A100, 062514 (2019)

work page 2019
[43]

Fläschner, B

N. Fläschner, B. S. Rem, M. Tarnowski, D. Vogel, D.- S. Lühmann, K. Sengstock, and C. Weitenberg, Science 352, 1091 (2016)

work page 2016
[44]

T. Li, L. Duca, M. Reitter, F. Grusdt, E. Demler, M. En- dres, M. Schleier-Smith, I. Bloch, and U. Schneider, Sci- ence352, 1094 (2016)

work page 2016
[45]

Xie, H.-X

B. Xie, H.-X. Wang, X. Zhang, P. Zhan, J.-H. Jiang, M. Lu, and Y. Chen, Nat. Rev. Phys.3, 520 (2021)

work page 2021
[46]

H. C. Po, H. Watanabe, and A. Vishwanath, Phys. Rev. Lett.121, 126402 (2018)

work page 2018
[47]

Okuma and M

N. Okuma and M. Sato, Annu. Rev. Condens. Matter Phys.14, 83 (2023)

work page 2023
[48]

M. A. J. Herrera and D. Bercioux, Commun. Phys.6, 42 (2023)

work page 2023
[49]

Sun, C.-R

W. Sun, C.-R. Yi, B.-Z. Wang, W.-W. Zhang, B. C. Sanders, X.-T. Xu, Z.-Y. Wang, J. Schmiedmayer, Y. Deng, X.-J. Liu, S. Chen, and J.-W. Pan, Phys. Rev. Lett.121, 250403 (2018)

work page 2018
[50]

C. Wang, P. Zhang, X. Chen, J. Yu, and H. Zhai, Phys. Rev. Lett.118, 185701 (2017). 6

work page 2017