pith. sign in

arxiv: 2505.12989 · v2 · submitted 2025-05-19 · ❄️ cond-mat.mtrl-sci

Ultrafast Laser Induces Macroscopic Symmetry-Breaking of Diamond Color Centers

Pith reviewed 2026-05-22 14:30 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords nitrogen-vacancy centerultrafast dynamicssymmetry breakingcoherent phononsJahn-Teller distortiondiamonddensity functional theoryelectron-phonon coupling
0
0 comments X

The pith

Ultrafast laser excitation breaks C3v symmetry in diamond NV centers and launches propagating coherent phonons.

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

The paper simulates the response of negatively charged nitrogen-vacancy centers in diamond to an ultrafast laser pulse using real-time time-dependent density functional theory. It establishes that minority-spin electrons are promoted within 100 fs, creating a charge ordering that breaks the C3v symmetry of the defect. Ionic motion on the resulting excited-state potential energy surface then produces both symmetric carbon-nitrogen bond oscillations and dynamic Jahn-Teller distortions that further break C3v symmetry; these distortions drive nonlocal coherent phonons that travel through the diamond lattice at the speed of sound while carrying the broken symmetry. A reader would care because NV centers are leading platforms for quantum sensing and information, so direct atomic-scale insight into how light, vibrations, and spins interact on femtosecond timescales can guide better control of their quantum properties.

Core claim

Laser excitation promotes minority-spin electrons within 100 fs, establishing a C3v-symmetry breaking charge ordering. Subsequently, ionic motion on the potential energy surface of the excited electrons generates both symmetric oscillations of carbon-nitrogen bonds and dynamic Jahn-Teller distortions with a C3v-symmetry breaking. These distortions subsequently induce nonlocal coherent phonons in the diamond lattice, which propagate with the C3v-symmetry breaking at the sound velocity (~2 Å/fs).

What carries the argument

Real-time time-dependent density functional theory simulation of the coupled electron-phonon-spin dynamics that directly tracks symmetry-breaking charge ordering and the subsequent ionic distortions.

If this is right

  • The broken symmetry appears first from electron promotion and then again from ionic motion, on a total timescale of hundreds of femtoseconds.
  • The induced coherent phonons carry the C3v symmetry breaking outward through the lattice at the sound velocity.
  • The simulations supply a direct time-resolved picture of how electrons, phonons, and spins couple inside the defect.
  • These processes occur on timescales short enough to affect spin coherence during quantum operations.

Where Pith is reading between the lines

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

  • The same two-stage symmetry-breaking sequence may appear in other point defects under strong optical excitation.
  • The propagating phonons could serve as a fast channel to couple distant NV centers or to modulate their local strain environment.
  • Ultrafast X-ray or electron diffraction experiments tuned to the predicted ~2 Å/fs wavefront could directly image the moving symmetry breaking.
  • If the phonons preserve coherence over longer distances, they might enable hybrid quantum-acoustic devices based on NV centers.

Load-bearing premise

The chosen density-functional approximation and basis set reproduce the excited-state potential energy surface and the resulting ionic trajectories without large errors that would change the predicted symmetry breaking or phonon propagation.

What would settle it

Time-resolved diffraction or optical measurements that detect no C3v symmetry lowering or no coherent phonons propagating at ~2 Å/fs after femtosecond laser excitation of NV centers would falsify the central sequence of events.

Figures

Figures reproduced from arXiv: 2505.12989 by Chao-Bo Liu, Chao Lian, Qi Xiao, Qi-Zheng Ji, Yang Gao.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic picture of the [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

The negatively charged nitrogen-vacancy center is a leading quantum platform due to its excellent spin coherence and stable interactions. Understanding its ultrafast dynamics is crucial for quantum applications but presents significant challenges for both experimental characterization and atomic-scale modeling. Here, we employ real-time time-dependent density functional theory to investigate the coupled electron-phonon-spin dynamics in negatively charged nitrogen-vacancy centers. Laser excitation promotes minority-spin electrons within 100~fs, establishing a $C_{3v}$-symmetry breaking charge ordering. Subsequently, ionic motion on the potential energy surface of the excited electrons generates both symmetric oscillations of carbon-nitrogen bonds and dynamic Jahn-Teller distortions with a $C_{3v}$-symmetry breaking. These distortions subsequently induce nonlocal coherent phonons in the diamond lattice, which propagate with the $C_{3v}$-symmetry breaking at the sound velocity ($\sim$2~\AA/fs). Our simulations provide direct time-resolved visualization of these processes, offering novel insights into the microscopic interplay of electrons, phonons, and spins in nitrogen-vacancy centers.

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 employs real-time time-dependent density functional theory (RT-TDDFT) simulations to study laser-driven ultrafast dynamics in negatively charged nitrogen-vacancy (NV-) centers in diamond. It claims that laser excitation promotes minority-spin electrons within 100 fs, producing a C_{3v}-symmetry-breaking charge ordering. Subsequent ionic relaxation on the excited-state potential energy surface generates symmetric carbon-nitrogen bond oscillations together with dynamic Jahn-Teller distortions that also break C_{3v} symmetry; these distortions launch nonlocal coherent phonons that propagate through the diamond lattice while preserving the symmetry breaking at the sound velocity of approximately 2 Å/fs. The work emphasizes direct time-resolved visualization of the coupled electron-phonon-spin processes.

Significance. If the reported symmetry-breaking mechanism and phonon propagation are robust, the study supplies a microscopic, time-resolved picture of how optical excitation couples to lattice and spin degrees of freedom in a leading quantum defect. Such insight is relevant for ultrafast control protocols in NV-based quantum sensing and information processing. The forward-simulation character of the work, which generates the symmetry breaking and phonon velocity as outputs rather than inputs, is a methodological strength.

major comments (2)
  1. Abstract and Methods: The central claim that RT-TDDFT produces a C_{3v}-symmetry-breaking charge ordering via minority-spin promotion within 100 fs rests on the accuracy of the chosen exchange-correlation functional and basis for the localized NV defect states. Semi-local functionals are known to suffer from self-interaction and delocalization errors that can distort spin densities and excited-state potential energy surfaces; the manuscript provides no benchmarks against hybrid functionals, GW, or experimental ultrafast spectra. This validation gap is load-bearing for the reported time scales and symmetry-breaking mechanism.
  2. Abstract: The phonon propagation velocity of ~2 Å/fs with preserved C_{3v} symmetry breaking is presented as a key result. The manuscript does not specify how this velocity was extracted from the ionic trajectories (e.g., via Fourier analysis or real-space propagation) nor does it compare the value to the independently known longitudinal sound velocity in diamond, leaving the quantitative claim without an internal consistency check.
minor comments (2)
  1. Notation for the point-group symmetry should be uniformly rendered as C_{3v} (with subscript) in both text and equations.
  2. The abstract and methods would benefit from explicit statements of the supercell size, k-point sampling, time step, and laser-pulse parameters to allow independent reproduction of the reported dynamics.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback on our manuscript. We address each major comment below with specific responses and indicate where revisions will be made to improve clarity and rigor.

read point-by-point responses
  1. Referee: Abstract and Methods: The central claim that RT-TDDFT produces a C_{3v}-symmetry-breaking charge ordering via minority-spin promotion within 100 fs rests on the accuracy of the chosen exchange-correlation functional and basis for the localized NV defect states. Semi-local functionals are known to suffer from self-interaction and delocalization errors that can distort spin densities and excited-state potential energy surfaces; the manuscript provides no benchmarks against hybrid functionals, GW, or experimental ultrafast spectra. This validation gap is load-bearing for the reported time scales and symmetry-breaking mechanism.

    Authors: We agree that functional choice is important for quantitative accuracy in spin and charge densities. Our simulations use the PBE functional in RT-TDDFT, selected for its balance of computational cost and prior success in describing NV- defect states and phonon dynamics in diamond. We acknowledge the absence of direct benchmarks against hybrids or GW in the current manuscript. Performing hybrid-functional RT-TDDFT for the required supercell sizes and femtosecond timescales remains computationally prohibitive. In revision we will expand the Methods section to discuss known limitations of semi-local functionals for localized defect states, cite literature validating PBE for NV- optical and spin properties, and add a brief note on the robustness of the observed symmetry-breaking mechanism across tested initial conditions. We will also reference available experimental ultrafast spectroscopy data for qualitative comparison of timescales. revision: partial

  2. Referee: Abstract: The phonon propagation velocity of ~2 Å/fs with preserved C_{3v} symmetry breaking is presented as a key result. The manuscript does not specify how this velocity was extracted from the ionic trajectories (e.g., via Fourier analysis or real-space propagation) nor does it compare the value to the independently known longitudinal sound velocity in diamond, leaving the quantitative claim without an internal consistency check.

    Authors: We thank the referee for highlighting this omission. The velocity was obtained by tracking the spatial position of the leading edge of the propagating coherent phonon wavefront in the real-space ionic displacement fields as a function of simulation time. In the revised manuscript we will add an explicit description of this extraction procedure in the Results section, including a supplementary figure showing the wavefront position versus time. We will also include a direct comparison to the known longitudinal sound velocity in diamond (~1.8–2.0 Å/fs), which confirms internal consistency with the simulated value of approximately 2 Å/fs. revision: yes

Circularity Check

0 steps flagged

No circularity: RT-TDDFT forward simulation outputs are independent of reported observables

full rationale

The paper reports results from real-time TDDFT simulations of laser-driven electron-phonon-spin dynamics in NV centers. The central claims—minority-spin promotion within 100 fs, C3v charge ordering, dynamic Jahn-Teller distortions, and coherent phonons propagating at ~2 Å/fs—are direct outputs of the computational model applied to the system. No equations, fitted parameters, or self-citations are shown that reduce these quantities to inputs by construction. The derivation is a standard forward simulation whose results are not presupposed in the setup or definitions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The simulation rests on standard DFT approximations whose accuracy for this excited-state, spin-dependent process is assumed rather than demonstrated in the provided abstract.

axioms (1)
  • domain assumption Real-time time-dependent density functional theory with the chosen exchange-correlation functional and pseudopotentials faithfully captures the ultrafast electron promotion, charge ordering, and subsequent ionic dynamics in the NV- center.
    This is the central modeling premise invoked by the choice of computational method.

pith-pipeline@v0.9.0 · 5729 in / 1457 out tokens · 44975 ms · 2026-05-22T14:30:48.253890+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.

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    Laser excitation promotes minority-spin electrons within 100 fs, establishing a C3v-symmetry breaking charge ordering... dynamic Jahn-Teller distortions with a C3v-symmetry breaking... propagate with the C3v-symmetry breaking at the sound velocity (~2 Å/fs).

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We used the Perdew-Burke-Ernzerhof (PBE) exchange-correlation (XC) functional... plane-wave energy cutoff was set to 80 Ry... electron timestep δt is 4×10^{-4} a.u.

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

42 extracted references · 42 canonical work pages · 1 internal anchor

  1. [1]

    C. L. Degen, F. Reinhard, and P. Cappellaro, Quantum sensing, Reviews of Modern Physics89, 035002 (2017)

  2. [2]

    J. F. Barry, J. M. Schloss, E. Bauch, M. J. Turner, C. A. Hart, L. M. Pham, and R. L. Walsworth, Sensitivity op- timization for NV-diamond magnetometry, Reviews of Modern Physics92, 015004 (2020)

  3. [3]

    Wolfowicz, F

    G. Wolfowicz, F. J. Heremans, C. P. Anderson, S. Kanai, H. Seo, A. Gali, G. Galli, and D. D. Awschalom, Quan- tum guidelines for solid-state spin defects, Nature Re- views Materials6, 906 (2021)

  4. [4]

    Rovny, S

    J. Rovny, S. Gopalakrishnan, A. C. B. Jayich, P. Maletinsky, E. Demler, and N. P. de Leon, Nanoscale diamond quantum sensors for many-body physics, Na- ture Reviews Physics , 1 (2024)

  5. [5]

    J. Du, F. Shi, X. Kong, F. Jelezko, and J. Wrachtrup, Single-molecule scale magnetic resonance spectroscopy using quantum diamond sensors, Reviews of Modern Physics96, 025001 (2024)

  6. [6]

    Gali, Ab initio theory of the nitrogen-vacancy center in diamond, Nanophotonics8, 1907 (2019)

    ´A. Gali, Ab initio theory of the nitrogen-vacancy center in diamond, Nanophotonics8, 1907 (2019)

  7. [7]

    Gali, Recent advances in the ab initio theory of solid- state defect qubits, Nanophotonics12, 359 (2023)

    ´A. Gali, Recent advances in the ab initio theory of solid- state defect qubits, Nanophotonics12, 359 (2023)

  8. [8]

    V. M. Huxter, T. A. A. Oliver, D. Budker, and G. R. Fleming, Vibrational and electronic dynamics of nitrogen–vacancy centres in diamond revealed by two- dimensional ultrafast spectroscopy, Nature Physics9, 744 (2013)

  9. [9]

    Ulbricht, S

    R. Ulbricht, S. Dong, I.-Y. Chang, B. M. K. Mariserla, K. M. Dani, K. Hyeon-Deuk, and Z.-H. Loh, Jahn-Teller- induced femtosecond electronic depolarization dynamics of the nitrogen-vacancy defect in diamond, Nature Com- munications7, 13510 (2016)

  10. [10]

    W. P. Carbery, C. A. Farfan, R. Ulbricht, and D. B. Turner, The phonon-modulated Jahn–Teller distortion of the nitrogen vacancy center in diamond, Nature Commu- nications15, 8646 (2024)

  11. [11]

    Ichikawa, J

    T. Ichikawa, J. Guo, P. Fons, D. Prananto, T. An, and M. Hase, Cooperative dynamic polaronic picture of di- amond color centers, Nature Communications15, 7174 (2024)

  12. [12]

    Gali, Time-dependent density functional study on the excitation spectrum of point defects in semiconductors, physica status solidi (b)248, 1337 (2011)

    A. Gali, Time-dependent density functional study on the excitation spectrum of point defects in semiconductors, physica status solidi (b)248, 1337 (2011)

  13. [13]

    Y. Ma, M. Rohlfing, and A. Gali, Excited states of the negatively charged nitrogen-vacancy color center in dia- mond, Physical Review B81, 041204 (2010)

  14. [14]

    Thiering and A

    G. Thiering and A. Gali,Ab InitioMagneto-Optical Spectrum of Group-IV Vacancy Color Centers in Dia- mond, Physical Review X8, 021063 (2018)

  15. [15]

    De´ ak, B

    P. De´ ak, B. Aradi, T. Frauenheim, E. Janz´ en, and A. Gali, Accurate defect levels obtained from the HSE06 range-separated hybrid functional, Physical Review B 81, 153203 (2010)

  16. [16]

    A. Gali, T. Simon, and J. E. Lowther, Anab initiostudy of local vibration modes of the nitrogen-vacancy center in diamond, New Journal of Physics13, 025016 (2011)

  17. [17]

    Thiering and A

    G. Thiering and A. Gali, Ab initio calculation of spin- orbit coupling for an NV center in diamond exhibiting dy- namic Jahn-Teller effect, Physical Review B96, 081115 (2017)

  18. [18]

    Ulbricht and Z.-H

    R. Ulbricht and Z.-H. Loh, Excited-state lifetime of the$\mathrm{N}{\mathrm{V}}ˆ{\ensuremath{-}}$ infrared transition in diamond, Physical Review B98, 094309 (2018)

  19. [19]

    A. Gali, E. Janz´ en, P. De´ ak, G. Kresse, and E. Kaxi- ras, Theory of Spin-Conserving Excitation of the $N\ensuremath{-}{V}ˆ{\ensuremath{-}}$Center in Di- amond, Physical Review Letters103, 186404 (2009)

  20. [20]

    Iv´ ady, I

    V. Iv´ ady, I. A. Abrikosov, and A. Gali, First princi- ples calculation of spin-related quantities for point de- fect qubit research, npj Computational Materials4, 76 (2018)

  21. [21]

    Wirtitsch, G

    D. Wirtitsch, G. Wachter, S. Reisenbauer, M. Gulka, V. Iv´ ady, F. Jelezko, A. Gali, M. Nesladek, and M. Trupke, Exploiting ionization dynamics in the ni- trogen vacancy center for rapid, high-contrast spin, and charge state initialization, Physical Review Research5, 013014 (2023)

  22. [22]

    Kollarics, B

    S. Kollarics, B. G. M´ arkus, R. Kucsera, G. Thiering, ´A. Gali, G. N´ emeth, K. Kamar´ as, L. Forr´ o, and F. Si- mon, Terahertz emission from diamond nitrogen-vacancy centers, Science Advances10, eadn0616 (2024)

  23. [23]

    C. Lian, M. Guan, S. Hu, J. Zhang, and S. Meng, Pho- toexcitation in Solids: First-Principles Quantum Simula- tions by Real-Time TDDFT, Advanced Theory and Sim- ulations1, 1800055 (2018)

  24. [24]

    Momentum-resolved TDDFT algorithm in atomic basis for real time tracking of electronic excitation

    C. Lian, S.-Q. Hu, M.-X. Guan, and S. Meng, Momentum-resolved TDDFT algorithm in atomic ba- sis for real time tracking of electronic excitation, The Journal of Chemical Physics149, 154104 (2018), arXiv:1702.08163

  25. [25]

    Lian, S.-J

    C. Lian, S.-J. Zhang, S.-Q. Hu, M.-X. Guan, and S. Meng, Ultrafast charge ordering by self-amplified exciton–phonon dynamics in TiSe2, Nature Communi- cations11, 43 (2020)

  26. [26]

    Elliott, K

    P. Elliott, K. Krieger, J. K. Dewhurst, S. Sharma, and E. K. U. Gross, Optimal control of laser-induced spin– orbit mediated ultrafast demagnetization, New Journal of Physics18, 013014 (2016)

  27. [27]

    Shokeen, M

    V. Shokeen, M. Sanchez Piaia, J.-Y. Bigot, T. M¨ uller, P. Elliott, J. K. Dewhurst, S. Sharma, and E. K. U. Gross, 7 Spin Flips versus Spin Transport in Nonthermal Elec- trons Excited by Ultrashort Optical Pulses in Transition Metals, Physical Review Letters119, 107203 (2017)

  28. [28]

    Wang, Y.-F

    S.-F. Wang, Y.-F. Hsu, J.-C. Pu, J. C. Sung, and L. G. Hwa, Determination of acoustic wave velocities and elas- tic properties for diamond and other hard materials, Ma- terials Chemistry and Physics85, 432 (2004)

  29. [29]

    Meng and E

    S. Meng and E. Kaxiras, Real-time, local basis-set im- plementation of time-dependent density functional the- ory for excited state dynamics simulations, The Journal of Chemical Physics129, 054110 (2008)

  30. [30]

    Wang, S.-S

    Z. Wang, S.-S. Li, and L.-W. Wang, Efficient Real-Time Time-Dependent Density Functional Theory Method and its Application to a Collision of an Ion with a 2D Mate- rial, Physical Review Letters114, 063004 (2015)

  31. [31]

    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...

  32. [32]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavaz- zoni, D. Ceresoli, M. c, N. Colonna, I. Carnimeo, 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. Kawa- mura, H.-Y. Ko, A. Kokalj, E. ...

  33. [33]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Physical Review Letters77, 3865 (1996)

  34. [34]

    D. R. Hamann, Optimized norm-conserving Vanderbilt pseudopotentials, Physical Review B88, 085117 (2013)

  35. [35]

    van Setten, M

    M. van Setten, M. Giantomassi, E. Bousquet, M. Ver- straete, D. Hamann, X. Gonze, and G.-M. Rignanese, The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table, Com- puter Physics Communications226, 39 (2018)

  36. [36]

    H. J. Monkhorst and J. D. Pack, Special points for Brillouin-zone integrations, Physical Review B13, 5188 (1976)

  37. [37]

    Popescu and A

    V. Popescu and A. Zunger, Extracting E versus k effec- tive band structure from supercell calculations on alloys and impurities, Physical Review B85, 085201 (2012)

  38. [38]

    Lian and S

    C. Lian and S. Meng, Dirac cone pairs in silicene induced by interface Si-Ag hybridization: A first-principles effec- tive band study, Physical Review B95, 245409 (2017)

  39. [39]

    Pacil` e, C

    D. Pacil` e, C. Cardoso, G. Avvisati, I. Vobornik, C. Mar- iani, D. A. Leon, P. Bonf` a, D. Varsano, A. Ferretti, and M. G. Betti, Narrowing of d bands of FeCo layers in- tercalated under graphene, Applied Physics Letters118, 121602 (2021)

  40. [40]

    Momma and F

    K. Momma and F. Izumi,VESTA 3for three- dimensional visualization of crystal, volumetric and mor- phology data, Journal of Applied Crystallography44, 1272 (2011)

  41. [41]

    T. A. Abtew, Y. Y. Sun, B.-C. Shih, P. Dev, S. B. Zhang, and P. Zhang, Dynamic Jahn-Teller Effect in the${\mathrm{NV}}ˆ{\ensuremath{-}}$Center in Di- amond, Physical Review Letters107, 146403 (2011)

  42. [42]

    Ultrafast Laser Induces Macroscopic Symmetry-Breaking of Diamond Color Centers

    M. S. Dresselhaus, G. Dresselhaus, and A. Jorio,Group Theory: Application to the Physics of Condensed Matter (Springer-Verlag, Berlin, 2008). ACKNOWLEDGMENT We acknowledge the financial supports from National Key Research and Development Program of China, Grant No. 2024YFA1408603 and CAS Project for Young Sci- entists in Basic Research under grant No. YSB...