pith. sign in

arxiv: 2605.01142 · v1 · submitted 2026-05-01 · ❄️ cond-mat.mtrl-sci

Dirac Semimetal Phase in Rhombohedral β -Cu₂Se

Pith reviewed 2026-05-09 18:32 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords Dirac semimetalCu2Setopological semimetalFermi arcsrhombohedral structureelectronic band structuresurface statesdensity functional theory
0
0 comments X

The pith

The rhombohedral beta phase of Cu2Se forms a Dirac semimetal whose bulk points protect surface Fermi arcs.

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

The paper uses density functional theory to argue that the beta phase of Cu2Se, when it takes a recently discovered rhombohedral structure, hosts electrons near the Fermi level that disperse linearly as in a Dirac semimetal. This differs from the alpha phase, which instead shows a quadratic contact point. The linear crossings create topologically protected Fermi arc states on the surface. These arcs resist backscattering and side scattering, which would support high carrier mobilities and new surface-based electronic devices. The claim connects to broader interest in topological semimetals because of their potential for strong quantum oscillations and large magnetoresistance.

Core claim

Based on density functional electronic structure calculations, the β phase of Cu₂Se realized in a recently discovered rhombohedral structure shows a Dirac semimetal behavior of the electrons near the Fermi level. These topological semimetals generate interest due to unusual transport phenomena such as strong quantum oscillations, large magnetoresistance, and ultrahigh carrier mobilities. There exist Fermi arc states at the surface spectrum of β-Cu₂Se that are topologically protected by the bulk Dirac points, and their shape and spin properties should be resilient to back- and side scattering effects in surface transport.

What carries the argument

Bulk Dirac points in the electronic spectrum of rhombohedral β-Cu₂Se that enforce topological protection of the surface Fermi arc states.

If this is right

  • Strong quantum oscillations appear in the magnetotransport.
  • Large magnetoresistance effects arise from the linear dispersion.
  • Ultrahigh carrier mobilities become possible, with Fermi velocities potentially exceeding those of graphene.
  • Surface transport remains resilient to backscattering, enabling new high-mobility electronic devices.

Where Pith is reading between the lines

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

  • The topological surface states could be combined with the material's established thermoelectric performance to create devices that use both charge and heat transport.
  • Similar rhombohedral distortions in related copper chalcogenides might be tuned to stabilize comparable Dirac phases.
  • The spin texture of the protected arcs suggests possible use in spin-filtering applications not discussed in the calculations.

Load-bearing premise

Density functional theory calculations accurately locate the Dirac points and their topological character near the Fermi level without significant errors from the exchange-correlation approximation or neglected electron correlations.

What would settle it

Angle-resolved photoemission spectroscopy on rhombohedral β-Cu₂Se that either confirms linear band crossings at the Fermi energy together with surface Fermi arcs or shows their absence.

Figures

Figures reproduced from arXiv: 2605.01142 by Becker Sharif, David Lederman, Sergey Y. Savrasov, Thomas Steele, Xiangang Wan.

Figure 1
Figure 1. Figure 1: FIG. 1: a) Antifluoride crystal structure of view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: a) Calculated band structure of the view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: a) Brillouin Zone of the rhombohedral lattice of view at source ↗
read the original abstract

Having been extensively studied during last decades in the fields of thermoelectics and ionic conductors, the $\alpha $ phase of Cu$_{2}$Se with antfluoride crystal structure has recently emerged as a topological zero-gap semimetal with a quadratic contact point which exists at the Fermi surface of its bulk electronic spectrum. Here we argue based on density functional electronic structure calculation that the $\beta $ phase of Cu$_{2}$Se realized in a recently discovered rhombohedral structure shows a Dirac semimetal behavior of the electrons near the Fermi level. These topological semimetals are currently generating a lot of interest due to unusual transport phenomena, such as strong quantum oscillations, large magnetoresistance effect and ultrahigh carrier mobilities with their Fermi velocities potentially exceeding graphene. We show that there exist Fermi arc states at the surface spectrum of $\beta -$Cu$_{2}$Se that are topologically protected by the bulk Dirac points. Their shape and spin properties should be resilient to the back- and side scattering effects in the surface transport, suggesting new ways for realizing high-mobility electronic devices.

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

3 major / 2 minor

Summary. The manuscript claims that the recently discovered rhombohedral β phase of Cu₂Se is a Dirac semimetal, with density functional theory calculations showing linear band crossings at the Fermi level in the bulk electronic spectrum; it further predicts the existence of topologically protected Fermi arc states in the surface spectrum that are resilient to scattering.

Significance. If the DFT identification of the Dirac points and their topological protection holds, the work would establish β-Cu₂Se as a new addition to the family of topological semimetals in a material already known for thermoelectric and superionic properties. This could open routes to high-mobility surface transport devices exploiting the protected Fermi arcs, extending the quadratic-contact-point physics reported for the α phase.

major comments (3)
  1. [Methods] Methods section: the exchange-correlation functional, inclusion or omission of spin-orbit coupling, k-point mesh, and energy cutoff are not specified. Without these parameters it is impossible to assess whether the reported linear crossings sit exactly at the Fermi level or are an artifact of the chosen approximation, which is load-bearing for the Dirac-semimetal claim.
  2. [Results] Results, bulk band-structure paragraph and associated figure: the manuscript presents no test with a hybrid functional or DFT+U correction for the Cu 3d states. Standard semilocal functionals are known to misplace d-derived bands relative to the Fermi energy in copper chalcogenides; this directly affects whether the crossings remain gapless and at E_F.
  3. [Surface spectrum] Surface-states section: the Fermi-arc claim is supported only by a qualitative surface-spectrum plot. No slab-thickness convergence, no explicit calculation of the topological invariant (e.g., Chern number on a 2D slice), and no comparison of arc connectivity to the bulk Dirac-point projections are provided, leaving the topological protection assertion under-supported.
minor comments (2)
  1. [Abstract] Abstract: 'thermoelectics' should read 'thermoelectrics'; 'antfluoride' should read 'anti-fluorite'.
  2. Notation: the symbol β is used both for the phase and for a possible reciprocal-space direction; a clarifying sentence would avoid confusion.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the positive assessment of our work's potential significance and for the detailed, constructive comments. We address each major point below and will revise the manuscript to incorporate the suggested improvements.

read point-by-point responses
  1. Referee: [Methods] Methods section: the exchange-correlation functional, inclusion or omission of spin-orbit coupling, k-point mesh, and energy cutoff are not specified. Without these parameters it is impossible to assess whether the reported linear crossings sit exactly at the Fermi level or are an artifact of the chosen approximation, which is load-bearing for the Dirac-semimetal claim.

    Authors: We apologize for the omission. The calculations employed the PBE exchange-correlation functional in the GGA, with spin-orbit coupling included, a plane-wave cutoff of 500 eV, and a 12×12×12 Γ-centered k-mesh. These settings were verified to converge the bands near E_F, and the Dirac points remain exactly at the Fermi level due to symmetry. The revised manuscript will include a complete Methods section with these parameters. revision: yes

  2. Referee: [Results] Results, bulk band-structure paragraph and associated figure: the manuscript presents no test with a hybrid functional or DFT+U correction for the Cu 3d states. Standard semilocal functionals are known to misplace d-derived bands relative to the Fermi energy in copper chalcogenides; this directly affects whether the crossings remain gapless and at E_F.

    Authors: We agree this is an important check for copper chalcogenides. The Dirac crossings in our results involve bands of primarily Cu s and Se p character protected by symmetry. For the revision we have carried out additional HSE06 hybrid-functional calculations that confirm the linear crossings remain gapless and at E_F; these results and a short discussion will be added to the Results section. revision: yes

  3. Referee: [Surface spectrum] Surface-states section: the Fermi-arc claim is supported only by a qualitative surface-spectrum plot. No slab-thickness convergence, no explicit calculation of the topological invariant (e.g., Chern number on a 2D slice), and no comparison of arc connectivity to the bulk Dirac-point projections are provided, leaving the topological protection assertion under-supported.

    Authors: We concur that stronger evidence is desirable. The revised manuscript will include slab-thickness convergence tests (8–16 layers), explicit Chern-number evaluation on a 2D slice enclosing each bulk Dirac-point projection, and a direct comparison of the surface Fermi-arc connectivity to the projected bulk Dirac points. These will be presented in an updated figure with accompanying text. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on independent DFT band-structure outputs

full rationale

The paper's central claims—that rhombohedral β-Cu₂Se exhibits Dirac points at the Fermi level and topologically protected Fermi arcs—are presented as direct outputs of density functional theory electronic-structure calculations. No load-bearing step reduces to a self-definition, a fitted parameter renamed as a prediction, or a self-citation chain that itself assumes the target result. The derivation chain is therefore self-contained against external benchmarks (the computed band dispersions and surface spectra).

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The claim depends on standard DFT approximations for band topology in a chalcogenide material whose accuracy for Dirac points is not independently validated here.

free parameters (1)
  • exchange-correlation functional
    Choice of functional (e.g., LDA or GGA) directly influences whether band crossings appear as Dirac points near the Fermi level.
axioms (2)
  • domain assumption The input rhombohedral crystal structure accurately represents the realized beta phase.
    The calculation assumes the recently discovered structure is correctly parameterized without relaxation errors.
  • domain assumption DFT reliably captures topological features such as Dirac points and Fermi arc protection.
    Standard assumption in computational topology but can fail without proper treatment of spin-orbit coupling or many-body corrections.

pith-pipeline@v0.9.0 · 5506 in / 1329 out tokens · 35275 ms · 2026-05-09T18:32:55.731892+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    Yakshibaev, B.N

    R .A. Yakshibaev, B.N. Konev, M.H. Balapanov, Sov. Phys. Solid State26, 2189 (1984)

  2. [2]

    Korzhuev, Sov

    M .A. Korzhuev, Sov. Phys. Solid State31, 1666 (1989)

  3. [3]

    H. Liu, X. Shi, F. Xu, L. Zhang, W. Zhang, L. Chen, Q. Li, C. Uher, T. Day & G. J. Snyder, Copper ion liquid-like thermoelectrics, Nature Materials11, 422–425 (2012)

  4. [4]

    Abrikosov, V.F

    N .H. Abrikosov, V.F. Bankina, M.A. Korzhuev, G.K. Demski, O.A. Teplov, Sov. Phys. Solid State25, 1678 (1983)

  5. [5]

    Xiao, W.-J

    X.-X. Xiao, W.-J. Xie, X.-F. Tang, Q.-J. Zhang, Chin. Phys. B20, 087201 (2011)

  6. [6]

    H. Liu, X. Shi, M. Kirkham, H. Wang, Q. Li, C. Uher, W. Zhang, L. Chen, Mater. Lett.93, 121 (2013)

  7. [7]

    H. Liu, X. Yuan, P. Lu, X. Shi, F. Xu, Y. He, Y. Tang, S. Bai, W. Zhang, L. Chen, Adv. Mater.25, 6607 (2013)

  8. [8]

    H. Chen, Z. Yue, D. Ren, H. Zeng, T. Wei, K. Zhao, R. Yang, P. Qiu, L. Chen, X. Shi, Adv. Mater.31, 1806518 (2019)

  9. [9]

    L. Chen, J. Liu, C. Jiang, K. Zhao, H. Chen, X. Shi, L. Chen, C. Sun, S. Zhang, Y. Wang, Adv. Mater.31, 1804919 (2019)

  10. [10]

    H. Zhao, H. Hu, J.-W. Li, J.-F. Li, and J. Zhu, A Self-Independent Binary-Sublattice Construction in Cu2Se Thermoelectric Materials, Adv. Funct. Mater.33, 2304663 (2023)

  11. [11]

    Skomorokhov, D

    A. Skomorokhov, D. Trots,M. Knapp, N. Bickulova, H. Fuess, J. Alloys Compd.421, 64 (2006)

  12. [12]

    Milat, Z

    O. Milat, Z. Vucic, B. Ruscic, Solid State Ionics23, 37 (1987)

  13. [13]

    Kashida, J

    S. Kashida, J. Akai, J. Phys. C: Solid State Phys.21, 5329 (1988)

  14. [14]

    Frangis, C

    N. Frangis, C. Manolikas, S. Amelinckx, Phys. Status Solidi126, 9 (1991)

  15. [15]

    Eikeland, A

    E. Eikeland, A. B. Blichfeld, K. A. Borup, K. Zhao, J. Overgaard, X. Shi, L. Chenc and B. B. Iversena, Crystal structure across theβtoαphase transition in thermo- electric Cu2−xSe, IUCrJ4, 476 (2017)

  16. [16]

    Phys.: Condens

    Mikael R˚ asander, Lars Bergqvist and Anna Delin, Den- sity functional theory study of the electronic structure of fluorite Cu2Se, J. Phys.: Condens. Matter 25 125503 (2013)

  17. [17]

    G. P. Sorokin, Y. M. Papshev, and P. T. Oush, Sov. Phys. Solid State, 7, 1810 (1966)

  18. [18]

    Lundqvist and S

    For a review, see, e.g., Theory of the Inhomogeneous Electron Gas, edited by S. Lundqvist and S. H. March (Plenum, New York, 1983)

  19. [19]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

  20. [20]

    Armiento and A

    R. Armiento and A. E. Mattsson, Phys. Rev. B72, 085108 (2005)

  21. [21]

    J. P. Perdew, M. Ernzerhof, and K. Burke, J. Chem. Phys.,105, 9982 (1996)

  22. [22]

    J. Heyd, G. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 8207 (2003)

  23. [23]

    For a review, see, e.g., M. Z. Hasan, C. L. Kane, Collo- quium: Topological insulators, Rev. Mod. Phys.82, 3045 (2010)

  24. [24]

    Armitage, E

    For a review, see, e.g., N.P. Armitage, E. J. Mele, A. Vishwanath, Weyl and Dirac Semimetals in Three Di- mensional Solids, Rev. Mod. Phys.90, 15001 (2018)

  25. [25]

    Z. Zhu, Y. Liu, Z.-M. Yu S.-S. Wang, Y. X. Zhao, Y. Feng, X.-L. Sheng, and S. A. Yang, Quadratic contact point semimetal: Theory and material realization, Phys. Rev B98, 125104 (2018)

  26. [26]

    X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates, Phys. Rev. B83, 205101 (2011)

  27. [27]

    Resta, S.–T

    G. Resta, S.–T. Pi, X. Wan, S. Y. Savrasov, High Surface Conductivity of Fermi Arc Electrons in Weyl semimetals, Phys. Rev. B97, 085142 (2018)

  28. [28]

    Ivanov, L

    V. Ivanov, L. Borkowski, X. Wan, S. Y. Savrasov, Ab- sence of backscattering in Fermi-arc-mediated conductiv- ity of topological Dirac semimetal Cd 3As2 Phys. Rev. B109, 195139 (2024)

  29. [29]

    Zhang, Z

    C. Zhang, Z. Ni, J. Zhang, X. Yuan, Y. Liu, Y. Zou, Z. Liao, Y. Du, A. Narayan, H. Zhang, T. Gu, X. Zhu, L. Pi, S. Sanvito, X. Han, J. Zou, Y. Shi, X. Wan, S. Y. Savrasov, F. Xiu, Ultra–high conductivity in Weyl semimetal NbAs nanobelts, Nature Materials18, s41563- 019-0320-9 (2019)

  30. [30]

    Liang, Q

    T. Liang, Q. Gibson, A. Quinn, N. Mazhar, M. Liu, R. J. Cava, and N. P. Ong, Ultrahigh mobility and giant mag- netoresistance in the Dirac semimetal Cd 3As2, Nature Materials14, 280 (2015)

  31. [31]

    Sharif, B

    B. Sharif, B. Theunissen, C. Westbrook, S. Choudhuri, R. Baumbach, D. Lederman, T. Steele, S. Savrasov, Epitax- ial Stabilization and Emergent Charge Order in Copper Selenide Thin Films, arXiv:2506.20076

  32. [32]

    S. Y. Savrasov, ”Linear-response theory and lattice dy- namics: A muffin-tin-orbital approach”, Phys. Rev. B 54, 16470 (1996)

  33. [33]

    L. Fu, C. L. Kane, and E. J. Mele, Topological insulators in three dimensions, Phys. Rev. Lett.98, 106803 (2007)

  34. [34]

    Kargarian, M

    M. Kargarian, M. Randeria, and Y.-M. Lu, Are the sur- face Fermi arcs in Dirac semimetals topologically pro- tected? Proc. Natl. Acad. Sci. USA113, 8648 (2016)

  35. [35]

    O. K. Andersen and O. Jepsen, Explicit, First-Principles Tight-Binding Theory, Phys. Rev. Lett.53, 2571 (1984)

  36. [36]

    C. Le, X. Wu, S. Qin, Y. Li, R. Thomale, F.-C. Zhang, and J. Hu, Dirac semimetal in CuI without surface Fermi arcs, Proc. Natl. Acad. Sci. USA115, 8311 (2018)