Citation:

Photonic topological fermi nodal disk in nonHermitian magnetic plasma


  • Light: Science & Applications  9, Article number: (2020)
More Information
  • Corresponding author:
    Yuanjiang Xiang (xiangyuanjiang@126.com)Shuang Zhang (s.zhang@bham.ac.uk)
  • Received: 04 March 2019
    Revised: 13 January 2020
    Accepted: 26 February 2020
    Published online: 11 March 2020

doi: https://doi.org/10.1038/s41377-020-0274-3

  • Topological physics mainly arises as a necessary link between properties of the bulk and the appearance of surface states, and has led to successful discoveries of novel topological surface states in Chern insulators, topological insulators, and topological Fermi arcs in Weyl, Dirac, and Nodal line semimetals owing to their nontrivial bulk topology. In particular, topological phases in non-Hermitian systems have attracted growing interests in recent years. In this work, we predict the emergence of the topologically stable nodal disks where the real part of the eigen frequency is degenerate between two bands in non-ideal magnetohydrodynamics plasma with collision and viscosity dissipations. Each nodal disk possesses continuously distributed topological surface charge density that integrates to unity. It is found that the lossy Fermi arcs at the interface connect to the middle of the projection of the nodal disks. We further show that the emergence, coalescence, and annihilation of the nodal disks can be controlled by plasma parameters and dissipation terms. Our findings contribute to understanding of the linear theory of bulk and surface wave dispersions of non-ideal warm magnetic plasmas from the perspective of topological physics.
  • 加载中
  • [1] Massines, F. et al. On the nature of the luminescence emitted by a polypropylene film after interaction with a cold plasma at low temperature. J. Appl. Phys. 81, 937-943 (1997). doi: 10.1063/1.364186
    [2] Hasegawa, A. Drift-wave instabilities of a compressional mode in a high-β Plasma. Phys. Rev. Lett. 27, 11-14 (1971). doi: 10.1103/PhysRevLett.27.11
    [3] Landau, L. D. & Lifshitz, E. M. Electrodynamics of Continuous Media. (Oxford: Pergamon Press, 1960).
    [4] Ilgisonis, V. I. et al. Geodesic acoustic modes and zonal flows in rotating large-aspect-ratio tokamak plasmas. Plasma Phys. Controlled Fusion 53, 065008 (2011). doi: 10.1088/0741-3335/53/6/065008
    [5] Hirshman, S. P. & Sigmar, D. J. Neoclassical transport of impurities in tokamak plasmas. Nucl. Fusion 21, 1079-1201 (1981). doi: 10.1088/0029-5515/21/9/003
    [6] Lindl, J. Development of the indirect‐drive approach to inertial confinement fusion and the target physics basis for ignition and gain. Plasma Plasmas 2, 3933-4024 (1995).
    [7] Pile, D. Plasma X-ray source. Nat. Photon. 9, 631 (2015).
    [8] Shivarova, A., Stoychev, T. & Russeva, S. Surface wave propagation along a current-carrying warm plasma column. J. Phys. D: Appl. Phys. 8, 383-393 (1975). doi: 10.1088/0022-3727/8/4/010
    [9] Aliev, Y. M. & Brodin, G. Instability of a strongly inhomogeneous plasma. Phys. Rev. A 42, 2374-2378 (1990). doi: 10.1103/PhysRevA.42.2374
    [10] Buti, B. Advances in Space Plasma Physics. (Singapore: World Scientific, 1985).
    [11] Gao, W. L. et al. Photonic Weyl degeneracies in magnetized plasma. Nat. Commun. 7, 12435 (2016). doi: 10.1038/ncomms12435
    [12] Wang, D. Y. et al. Photonic Weyl points due to broken time-reversal symmetry in magnetized semiconductor. Nat. Phys. 15, 1150-1155 (2019). doi: 10.1038/s41567-019-0612-7
    [13] Ida, K. Experimental studies of the physical mechanism determining the radial electric field and its radial structure in a toroidal plasma. Plasma Phys. Controlled Fusion 40, 1429-1488 (1998). doi: 10.1088/0741-3335/40/8/002
    [14] Kriesel, J. M. & Driscoll, C. F. Measurements of viscosity in pure-electron plasmas. Phys. Rev. Lett. 87, 135003 (2001). doi: 10.1103/PhysRevLett.87.135003
    [15] Shen, H. T., Zhen, B. & Fu, L. Topological band theory for non-Hermitian Hamiltonians. Phys. Rev. Lett. 120, 146402 (2018). doi: 10.1103/PhysRevLett.120.146402
    [16] Rudner, M. S. & Levitov, L. S. Topological transition in a non-Hermitian quantum walk. Phys. Rev. Lett. 102, 065703 (2009). doi: 10.1103/PhysRevLett.102.065703
    [17] Lee, T. E. Anomalous edge state in a non-Hermitian lattice. Phys. Rev. Lett. 116, 133903 (2016). doi: 10.1103/PhysRevLett.116.133903
    [18] Leykam, D. et al. Edge modes, degeneracies, and topological numbers in non-Hermitian systems. Phys. Rev. Lett. 118, 040401 (2017).
    [19] Gong, Z. P. et al. Topological phases of non-Hermitian systems. Phys. Rev. X 8, 031079 (2018).
    [20] Yao, S. Y. & Wang, Z. Edge states and topological invariants of non-Hermitian systems. Phys. Rev. Lett. 121, 086803 (2018). doi: 10.1103/PhysRevLett.121.086803
    [21] Chen, W. J. et al. Exceptional points enhance sensing in an optical microcavity. Nature 548, 192-196 (2017). doi: 10.1038/nature23281
    [22] Zhou, H. Y. et al. Observation of bulk Fermi arc and polarization half charge from paired exceptional points. Science 359, 1009-1012 (2018). doi: 10.1126/science.aap9859
    [23] Miri, M. A. & Alù, A. Exceptional points in optics and photonics. Science 363, eaar7709 (2019). doi: 10.1126/science.aar7709
    [24] El-Ganainy, R. et al. Non-Hermitian physics and PT symmetry. Nat. Phys. 14, 11-19 (2018). doi: 10.1038/nphys4323
    [25] Carlström, J. & Bergholtz, E. J. Exceptional links and twisted Fermi ribbons in non-Hermitian systems. Phys. Rev. A 98, 042114 (2018). doi: 10.1103/PhysRevA.98.042114
    [26] Berry, M. V. & Mondragon, R. J. Diabolical points in one-dimensional Hamiltonians quartic in the momentum. J. Phys. A: Math. Gen. 19, 873-885 (1986). doi: 10.1088/0305-4470/19/6/018
    [27] Berry, M. V. & Jeffrey, M. R. Conical diffraction complexified: dichroism and the transition to double refraction. J. Opt. A: Pure Appl. Opt. 8, 1043-1051 (2006). doi: 10.1088/1464-4258/8/12/003
    [28] Moiseyev, N. Non-Hermitian Quantum Mechanics. (Cambridge University Press, Cambridge, 2011).
    [29] Heiss, W. D. The physics of exceptional points. J. Phys. A: Math. Theor. 45, 444016 (2012). doi: 10.1088/1751-8113/45/44/444016
    [30] Cerjan, A. et al. Effects of non-Hermitian perturbations on Weyl Hamiltonians with arbitrary topological charges. Phys. Rev. B 97, 075128 (2018). doi: 10.1103/PhysRevB.97.075128
    [31] Kullig, J., Yi, C. H. & Wiersig, J. Exceptional points by coupling of modes with different angular momenta in deformed microdisks: a perturbative analysis. Phys. Rev. A 98, 023851 (2018). doi: 10.1103/PhysRevA.98.023851
    [32] Zhang, Z. Y. et al. Non-Hermitian optics in atomic systems. J. Phys. B: At. Mol. Opt. Phys. 51, 072001 (2018). doi: 10.1088/1361-6455/aaaf9f
    [33] Kullig, J. et al. Exceptional points of third-order in a layered optical microdisk cavity. N. J. Phys. 20, 083016 (2018). doi: 10.1088/1367-2630/aad594
    [34] Malzard, S. & Schomerus, H. Bulk and edge-state arcs in non-Hermitian coupled-resonator arrays. Phys. Rev. A 98, 033807 (2018). doi: 10.1103/PhysRevA.98.033807
    [35] Xu, Y., Wang, S. T. & Duan, L. M. Weyl exceptional rings in a three-dimensional dissipative cold atomic gas. Phys. Rev. Lett. 118, 045701 (2017). doi: 10.1103/PhysRevLett.118.045701
    [36] Keck, F., Korsch, H. J. & Mossmann, S. Unfolding a diabolic point: a generalized crossing scenario. J. Phys. A: Math. Theor. 36, 2125-2137 (2003).
    [37] Zyuzin, A. A. & Zyuzin, A. Y. Flat band in disorder-driven non-Hermitian Weyl semimetals. Phys. Rev. B 97, 041203 (2018). doi: 10.1103/PhysRevB.97.041203
    [38] González, J. & Molina, R. A. Topological protection from exceptional points in Weyl and nodal-line semimetals. Phys. Rev. B 96, 045437 (2017). doi: 10.1103/PhysRevB.96.045437
    [39] Doveil, F., Vosluisant, Y. & Tsunoda, S. I. Localization of langmuir waves in a fluctuating plasma. Phys. Rev. Lett. 69, 2074-2077 (1992). doi: 10.1103/PhysRevLett.69.2074
    [40] Lee, H. J. & Cho, S. H. Boundary conditions for surface waves propagating along the interface of plasma flow and free space. J. Plasma Phys. 58, 409-419 (1997). doi: 10.1017/S002237789700603X
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(4)

Article Metrics

Article views(320) PDF downloads(50) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

Photonic topological fermi nodal disk in nonHermitian magnetic plasma

  • 1. International Collaborative Laboratory of 2D Materials for Optoelectronic Science & Technology of Ministry of Education, Institute of Microscale Optoelectronics (IMO), Shenzhen University, 518060 Shenzhen, China
  • 2. School of Physics & Astronomy, University of Birmingham, Birmingham B15 2TT, UK
  • 3. National Key Laboratory for Laser Fusion, Research Centre of Laser Fusion, CAEP, 621900 Mianyang, China
  • 4. School of Physics and Electronics, Hunan University, 410082 Changsha, China
  • Corresponding author:

    Yuanjiang Xiang, xiangyuanjiang@126.com

    Shuang Zhang, s.zhang@bham.ac.uk

doi: https://doi.org/10.1038/s41377-020-0274-3

Abstract: Topological physics mainly arises as a necessary link between properties of the bulk and the appearance of surface states, and has led to successful discoveries of novel topological surface states in Chern insulators, topological insulators, and topological Fermi arcs in Weyl, Dirac, and Nodal line semimetals owing to their nontrivial bulk topology. In particular, topological phases in non-Hermitian systems have attracted growing interests in recent years. In this work, we predict the emergence of the topologically stable nodal disks where the real part of the eigen frequency is degenerate between two bands in non-ideal magnetohydrodynamics plasma with collision and viscosity dissipations. Each nodal disk possesses continuously distributed topological surface charge density that integrates to unity. It is found that the lossy Fermi arcs at the interface connect to the middle of the projection of the nodal disks. We further show that the emergence, coalescence, and annihilation of the nodal disks can be controlled by plasma parameters and dissipation terms. Our findings contribute to understanding of the linear theory of bulk and surface wave dispersions of non-ideal warm magnetic plasmas from the perspective of topological physics.

Introduction
  • Plasma, widely considered as the fourth states of matter, is a completely ionized gas consisting of freely moving ions, nuclei, and electrons. Plasma is believed to account for over 99% of the matters in the universe. Understanding the bulk and the surface wave dispersions of plasma's magnetic hydrodynamics is of great importance for both fundamental interests in topological matter states and broad applications in fusion energy harvesting, plasma diagnostics etc. In applications with low particle density and negligible temperature effect such as neon tubes1, known as "cold plasma", the charged particles are assumed to oscillate without significant inter-particle interactions due to the diminishing Lamor radius2. Landau3 considered magnetic fluids as continuous media and used the cold plasma theory to describe the excitations in plasma as frequency dependent, local (meaning no spatial dispersion) dielectric functions. Such electrodynamics of continuous media theory has been experimentally testified to give satisfying descriptions of zonal flows, Geodesic acoustic mode, and Alfven waves4.

    However, in higher energy regimes, such as thermonuclear weapons, tokamak configuration with magnetic and inertial confinement, and intense X-ray sources, the particle density is immensely higher, and the temperature effect cannot be neglected5-7. As a consequence, the dynamics of the plasma are thermal-dependant and described by the Navier-Stokes equations, and are known as "warm plasma", whose dielectric functions are not only frequency dependant but also spatially dispersive (nonlocal). Study of the bulk and surface wave dispersions of the "warm plasma" can be used for plasma diagnostics and for understanding astrophysical problems in the magnetosphere, the solar corona and surface waves enhanced electron heating, ion acceleration, synchrotron radiation, and static magnetic generation. Furthermore, surface waves are also involved in laser fusion, which is valuable for energy productions through controlled thermonuclear reactions (CTR)8-10.

    It has been shown previously that Weyl points may exist in cold plasma, leading to observation of interesting topological effects such as presence of Fermi arcs11, 12. Since large dissipation terms could arise in warm plasma13, 14, they provide an interesting platform for studying non-Hermitian topological physics, which has become a vibrant research field in recent years15-25. Topologically protected degeneracies in Hermitian systems such as Weyl points are stable, due to the participation of all three Pauli matrices in constructing the low-energy Hamiltonian. However, with additional degree of freedom from the loss/gain, the probability of finding a degeneracy is lower than is dictated by the Von Neumann-Wigner theorem26. This means the two band Hamiltonian may experience level repulsion lifting the degeneracies. For instance, the conical refraction phenomenon could be destroyed by introducing another "loss optical axis"27. However, interestingly, rather than simply lifting the degeneracy, the non-Hermicity may turn a Weyl point into a two-dimensional nodal disk where the real part of the complex eigen frequencies becomes degenerate, which is surrounded by an exceptional loop (EL), where both the complex eigen frequency and eigen-states coalesce. Till now, besides being theoretically investigated in a minimum 2-by-2 Hamiltonian model15, 28-31, there have been only a few studies on topological ELs in optical lattices and woodpile photonic crystals32-34.

    In this work, we predict the emergence of topologically stable nodal disks and ELs in non-ideal magnetohydrodynamics (MHD) plasma by taking into account the dissipations, including collisions and viscosity forces. While previous works have focused on ELs, we find the nodal disks to be more topologically important than the ELs mainly in two aspects: (1). The Berry curvatures are primarily emitted from the nodal disks. (2). The Fermi arcs connect to the centers of the nodal disks projected onto the interface. These nodal disks are non-Hermitian generalizations of the Weyl points whose topological charge can be confirmed by the distribution of non-Hermitian Berry curvature under bi-orthogonal basis15, 35-38. Momentum-space positions of nodal disks are sensitive to the plasma temperature, magnitude of the magnetic field, and the viscosity force. They could coalesce and annihilate through variation of these parameters, and this phenomenon has not been readily found in other systems. In the lossless limit, these nodal disks reduce into type-II Weyl points and could further transform into transitional Weyl points, which lies between type-I and type-II with the dispersion of one participating mode being perfectly flat, by reducing the plasma temperature to zero. Dissipative surface states (Fermi arcs) are also found to emanate from the topologically charged nodal disks. Interestingly, viscosity force is found to introduce extra type-I nodal disks at larger momentum. By tuning the magnitude of the viscosity loss, the extra type-I nodal disks could coalesce and annihilate with type-II nodal disks. Our study introduces the modern topological band theory into the century-long research field of warm plasma.

Results
  • In a warm plasma, the charged particles' macroscopic motion can be described by the fluid dynamics (Navier-Stokes equation) with parameters, such as pressure and density. In contrast to conventional fluid where only longitudinal wave exists due to the absence of tangential resilience, warm plasma's dispersion is more complicated due to the presence of electromagnetic force. In this work, MHD in combination with Maxwell equations are used to describe the fluid motions of the plasma under electromagnetic interactions. Taking into account the dissipation terms, which include the collision and viscosity, the plasma fluid dynamics equation is given as,

    $$ \frac{{\partial {\boldsymbol{u}}_{\boldsymbol{e}}}}{{\partial t}} + {\boldsymbol{u}}_{\boldsymbol{e}} \cdot \nabla {\boldsymbol{u}}_{\boldsymbol{e}} = q_en_{\mathrm{e}}\left( {{\boldsymbol{E}} + \frac{1}{c}{\boldsymbol{u}}_{\boldsymbol{e}} \times {\boldsymbol{B}}} \right) - \frac{{\nabla P_{\mathrm{e}}}}{{{\it{n}}_{\it{e}}{\it{m}}_{\it{{\mathrm{e}}}}}} + \frac{{{\boldsymbol{u}}_{\boldsymbol{e}}}}{\tau } + \eta \nabla ^2{\boldsymbol{u}}_{\boldsymbol{e}} $$ (1)

    Where ne is the plasma density, me is the electron mass, qe is the elementary charge, Pe is the electrons thermal pressure, η is the viscosity coefficient, which arises from the frictions between counter-propagating electrons, τ is the mean free time between electrons collision. The corresponding Hamiltonian formalism is given as (Supplementary Note 1 for detail).

    $$ \omega _p\left[ {\begin{array}{*{20}{c}} {} & { - {\boldsymbol{K}} \times /k_P} & {} & { - i} \\ {{\boldsymbol{K}} \times /k_P} & {} & {} & {} \\ {} & {} & {} & {\kappa {\boldsymbol{K}}/k_P} \\ i & {} & {\kappa {\boldsymbol{K}}/k_P} & {\omega _c\Delta - iI/\tau - iI\eta K^2} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\boldsymbol{E}} \\ {\boldsymbol{H}} \\ n \\ {\boldsymbol{j}} \end{array}} \right] = \omega \left[ {\begin{array}{*{20}{c}} {\boldsymbol{E}} \\ {\boldsymbol{H}} \\ n \\ {\boldsymbol{j}} \end{array}} \right] $$ (2)

    where $\omega _p = \left({n_0e^2/m_e} \right)^{1/2}$ is the equilibrium plasma frequency, $\omega _c = \frac{{q_eB}}{{m_ec}}$ is the electron cyclotron frequency, $k_p = \omega _p/c$, Δ=$\left[{\sigma _y\, 0;0\, 0} \right]$, I is the identity matrix, and ${\rm{\kappa }}^2 = \frac{{P_0\gamma }}{{n_0mc^2}}$ with γ being the ratio of specific heats, n is the fluctuation of plasma density, j the polarized current. In the rest of the paper, the frequencies are normalized by ωp, η has unit of ${\mathrm{kg}}/\left({{\mathrm{m}} \!\cdot\! {\mathrm{s}}^2} \right)$ corresponding to actual viscosity coefficient to η/ωp whose unit is ${\mathrm{kg}}/\left({{\mathrm{m}} \!\cdot\! {\mathrm{s}}} \right)$. The traditional mathmatically complicated calculation of the nonlocal dielectric functions is thus transformed into an eigen-value problem.

    We start by considering only the collision dissipation while neglecting the effect of viscosity. The dispersion of MHD plasma can be solved by diagonalizing the Hamiltonian in Eq. (1). The absolute value of the complex frequency of the non-ideal magnetized warm plasma is shown in Fig. 1a with $\omega _p = 3.5 \times 10^{12}\, {\mathrm{rad}}/{\mathrm{s}}$, $\omega _c = 2\omega _p$, κ=0.4, τ=2, corresponding to a magnetized plasma system with an electron density of $3.9 \times 10^{19}\, 1/{\mathrm{m}}^3$, an external magnetic field B of 8.01T, an electron thermal temperature of 27.3keV and an electronic collision frequency of $1.75 \times 10^{12}\, {\mathrm{rad}}/{\mathrm{s}}$. The band structure shows four bands, which, ordered from low to high frequencies, are the first right-handed circularly polarized band (R-wave), the longitudinal Langmuir (LM) wave (see detail in Supplementary Fig. 1), the left-handed circularly polarized band (L-wave) and the second band of R-wave. In the lossless limit R-wave and L-wave's dispersions are expressed as $k_{{\mathrm{R/L}}} = \frac{\omega }{c}\left[{1- \frac{{\omega _{pe}^2}}{{\omega ^2\left({1 \mp \omega _{ce}/\omega } \right)}}} \right]^{1/2}$. The LM wave's dispersion is given by $\omega ^2 = \omega _{pe}^2 + \kappa ^2c^2k^2$, which can reduce to the cold plasma's4 oscillation frequency by taking κ=0. (Supplementary Note 2 for detail and Supplementary Fig. 2). The positions of the Weyl points are caculated in Supplementary Fig. 3. To better understand these linear band crossing, we apply the p theory (effective Hamiltonian theory) under collisionless conditions to obtain the approximate Hamiltonian near the degeneracies. Expanding to first order in the vicinity of the outer degeneracy, we find the effective Hamiltonian as:

    $$ H_1 = \frac{{A + N}}{2}k_z\sigma _z + Mk_x\sigma _x - Mk_y\sigma _y + \frac{{A - N}}{2}Ik_z $$ (3)

    Fig. 1  Existence of Nodal disks in magnetohydrodynamic plasma.

    a Dispersion of the absolute value of eigen frequency ω with $\omega _c = 2, \, \kappa = 0.4, \, \tau = 2$. Along the kz axis, the bands are differentiated as right/left handedness circular polarization wave (R/L-wave) and Langmuir wave whose polarization is only in the z-direction. The bands along transverse directions (x-direction) have degeneracies denoted by the red/blue dots. Owing to the rotational symmetry in the plasma system, the dots in a form the loops in b. Akin to the Weyl points, there are two species of Nodal disks in the plasma system with opposite topological charge (+/-1 charge denoted by the red/blue color) that are connect by the spatial inversion operation. c-f Real and imaginary part of the eigen frequency in the kx-kz and kx-ky plane. Note that in the kx-kz plane, real and imaginary parts are simultaneously degenerated at isolated points away from the kz axis. While in the kx-ky plane, the degeneracies form a ring

    Fig. 2  Phase map, the emergence, and annihilations of the exceptional loops at τ=2.

    a The number of nodal disks is due to change by the strength of magnetic field ωc and the electron thermal parameter κ. Along the path crossing the border line, for example the path shown as the yellow arrow in a, the annihilation or emergence of nodal disks could be observed. Three-dimensional (3-D) band structures before and at the critical transition point are given in b and c. Difference of the real part of the eigen frequencies around the nodal disks before, at and after the annihilation of the nodal disks are given in d-f. The dashed lines represent the locations where the difference goes to zero. The positions of the nodal disks are located at the edge of the dashed lines. d Before the coalesce of the two nodal disks. e At the border between the two phases, the two nodal disks would coalesce and dissipate their topological charges. f After the coalesce, the zeros in the e would shift from zero

    Fig. 3  Berry curvature distribution and dissipative fermi arcs of the nodal disks with $\omega _c = 2, \, \kappa = 0.4, \, \tau = 2$.

    a Berry curvature for the upper band of collision warm magnetic plasma in the kx-kz plane with. The red rings are the surrounding of the Nodal disks in kx-ky plane and the blue solid line is the projection of the bulk nodal disk in kx-kz plane. Scales of the arrow represent the intensity of Berry curvatures and orientation indicates the direction. The arrows show that the bulk nodal disks are pure source of Berry curvature. b fermi arcs connect at the larger (1.96), equals (1.86), and smaller (1.76) bulk nodal disk frequency. For the Nodal disks frequency, Fermi arc connects at center of the real dependency line while Fermi arc emits from the bulk states when the frequency shifts. The color represents the value of the corresponding imaginary frequency. c dissipating fermi arc between plasma and vacuum, which connects the two bulk nodal disk with opposite chirality, the color represents the value of the corresponding imaginary frequency

    where $A = \frac{{\sqrt {\varepsilon _{12}} }}{{2\left({\varepsilon _{12}^2\, -\, \varepsilon _{12} \, +\, 2} \right)}}$, $N = \frac{{\omega _p^2\sqrt {\omega _c} }}{{2\left({\omega _p \, +\, \omega _c} \right)\sqrt {\omega _c \, -\, \omega _p} }}\kappa ^2$, $M = \frac{{\sqrt {\varepsilon _{12}} }}{{2\sqrt {\varepsilon _{12}^2 \, -\, \varepsilon _{12} \, +\, 2} }}$, and ε12 is sum of the diagonal and off-diagonal elements of the dielectric matrix evaluated at the Weyl's frequency, which can be expresses as $\omega _{{\mathrm{weyl}}} = \omega _p + \frac{{\omega _p^2\omega _c}}{{2\left({\omega _c^2 \, -\, \omega _p^2} \right)}}\kappa ^2$ (Supplementary Note 3 for detail). These linear degeneracies are Weyl points that function as the sources and drains of Berry curvature flux lines. Owing to the electrons collision, the new degeneracy points where both the real and the imaginary part of the complex frequency spectrum are degenerate shift away from the kz axis, as denoted by the colored dots in Fig. 1a. Owing to the rotational symmetry of the system around kz axis, the degeneracy points form exceptional loops in the kx-ky plane in Fig. 1b. Figure 1c, d show the real and imaginary parts of the eigen frequency in kx-kz direction, where the degeneracy line for the real part of the eigen frequency in the kx-kz plane forms a bulk nodal disk due to the rotational symmetry. The dispersion of the real and imaginary parts of the eigen frequency on the kx-ky plane are presented in Fig. 1e, f, respectively, which show the presence of an EL. As shown in Fig. 1, the modified dispersion of the LM wave is curved due to the thermal effect, providing the possibility of forming three pairs of nodal disks between the LM Wave and the L/R-wave. In contrast the cold plasma can at most possess two pairs of Weyl points due to its dispersionless LM mode39. As shown in Fig. 1a, two types of band-crossing exist: the band crossings between the LM wave and the L-wave at small k, and the other two between the LM wave and the first R-wave at larger k's. These band crossings are guaranteed by the polarization orthogonality between the longitudinal LM waves and the transverse L and R waves.

  • By tuning the thermal parameter κ and the cyclotron frequency ωc under a fixed collision coefficient τ=2, the nodal disks' location in the momentum space can be shifted along the kz direction, or even merged or annihilated. In Fig. 2a, the transition of the number of nodal disks is illustrated in a phase map in the parameter space (ωc, κ), which shows that the system acquires either three pairs or a single pair of nodal disks. Right at the phase boundary, the two nodal disks at larger kz join with each other. The difference in the real part of the frequency between the LM wave and R-wave constituting the nodal disks for different κ values is given in Fig. 2d-f. As shown in Fig. 2d, when κ=0.52, ωc=2 (below the critical transition value), the nodal disks are separate (white dash lines). While at the critical transition value κ=0.5228, ωc=2 the exceptional loops merge into each other and the nodal disks form a closed nodal surface, as shown in Fig. 2e. The corresponding 3-D band structures to Fig. 2d, e are given in Fig. 2b, c. It should be noted that even though the Chern number of the closed nodal surface is zero since the topological charges of the two nodal disks cancel each other, the Berry curvatures do not vanish because of the nonzero topological charge density on the surface. Further increase of κ beyond the critical value leads to simply shrinkage of the closed nodal surface and lifting of the degeneracy as is shown in Fig. 2f.

  • Similar to Weyl points in Hermitian systems, nodal disks are sources and drains of Berry curvatures in non-Hermitian systems. In such a dissipative system, Berry curvatures can be derived from any entry of the Berry connection matrix: $A_\mu ^{I, J} = i\langle I{\mathrm{|}}\partial _\mu {\mathrm{|}}J\rangle, \, I, J \in \left\{ {L, R} \right\}$. In Supplementary Note 5, it is shown that this matrix is Hermitian, whose diagonal entries are real, and off-diagonal complex. The matrix entries satisfy the relation: $A_\mu ^{L, L} = A_\mu ^{R, R}\, \ne\, A_\mu ^{L, R} = A_\mu ^{L, R \ast }$, where the sign * stands for the complex conjugate. It could be further shown that the Berry curvature calculated from different matrix entries are locally different, however, their integral over the nodal disk, the Chern numbers, are the same (for the complex term $A_\mu ^{L, R}$ and $A_\mu ^{L, R}$, only real parts are integrated)15. In dissipative MHD plasma, Berry curvature flux distribution from the nodal disk between the LM and the R waves is shown in Fig. 3a, where the red rings represent the ELs in kx-ky plane and the blue solid line is the projection of the nodal disk in kx-kz plane. Figure 3a clearly shows that the nodal disk is the source of Berry curvature. In proximity to the EL, the derived Berry curvatures diverage is at the rate of $1/\!\sqrt r$, where r is the distance to the exceptional loop (see detail in Supplementary Note 5). This means the Berry curvature emitted by the exceptional loop has neligible value. Hence it is confirmed that the topological charge is continuously distributed across the bulk nodal disks, rather than concentrated at the ELs.

    Owing to the nontrivial topological charge of the nodal disk, presence of Fermi arcs is expected. To solve for the Fermi arc, we consider an abrupt interface between the MHD plasma and vacuum, which can be found in magnetic confinement plasmas such as Tokamak configurations. Because of the nonlocality of the system, an additional boundary condition that the current density diminishes in the surface normal direction should be imposed, together with the conventional electromagnetic wave continuity conditions. Considering the continuity of bulk displacement current normal to the interface, this additional boundary condition is equivalent to continuity of Ex across the interface40. The dissipating Fermi arc near the outer the nodal disk at the real frequency parts of 1.76ωp, 1.86ωp, and 1.96ωp is given in Fig. 3b (more details are shown in Supplementary Fig. S7). Interestingly, at the EL frequency, the dissipative Fermi arc connects between the projected nodal disks right in the middle, as shown in Fig. 3c. However, away from the exception loop frequency, the Fermi arc could exist inside the bulk state continuum, as shown in Fig. 3b. Detailed derivations are provided in the Supplementary Note 4, and the surface modes at different frequencies are presented in Supplementary Fig. S4.

  • In the above, we have studied the warm magnetic plasma with collision loss. However, in tokamaks and helical devices, there also exists strong anomalous viscosity14. As indicated by Eq. (2), viscosity leads to an extra quadratic term in the dispersion of the LM wave, which becomes significant at large momentum. The significantly modified dispersion of the LM wave could introduce an extra pair of nodal disks by intersecting the R-wave again. The real part of the eigen frequency spectrum around the nodal disks is given in Fig. 4a. Comparing with that in a plasma with only collisions, the nodal disk in a plasma with viscosity dissipation is highly curved due to the presence of the higher order term in the momentum. The curved nodal disks can be attributed to the $k_{x, y}^2$ terms present in the effective Hamiltonian model (Supplementary Note 7). Increasing the strength of the viscosity force could induce coalescence and annihilation of the nodal disks as shown in Fig. 4b-d. Similar to the effects shown in previous sections, at a certain viscosity force, tuning the magnetic field strength and electron thermal parameter, the 2nd and 3rd nodal disks could annihilate, leaving only the 1st and the 4th exceptional loops.

    Fig. 4  Viscosity force induced extra Nodal disks and their phase transitions.

    a According to different viscosity parameter range at $\omega _c = 2, \, \kappa = 0.4$, the magnetohydrodynamic plasma could have either four pairs of Nodal disks or two pairs by the annihilations of the 2nd and the 3rd, or the 3rd and the 4th. White dash lines represent the real parts' degeneracy line in the momentum plane of x-y. Difference of the real part of the eigen frequencies b before, c at, and d after near to the exceptional loops' annihilation by increasing the viscosity force. The 3-D band structures are give in e-g
Discussion
  • In this paper, we have studied the nodal disks and the topological band theory in MHD plasma by considering the collision and viscosity effects. Electron thermal pressures drives the plasma oscillation (longitude mode) propagation and modifies the dispersion slope. Topological charged nodal disks are observed due to the dissipation effects arising from collision and viscosity. Topological protected dissipating fermi arcs are observed, which connects to the bulk nodal disk. Importantly, the emerge, coalesce, and annihilation of these nodal disks can be tuned by thermal parameter and the external magnetic field. Our discoveries could facilitate studies of both the linear theory of MHD plasma and the non-Hermitian topological band theories.

Methods
  • Plasma encompass a very large range of scales in length, density and temperature. Our theoretically predicted topological phase transitions could be observed in a broad range of the above parameters. For simplicity we normalize the cyclotron frequency ωc, the electron collision frequency 1/τ and the viscosity coefficient η by the plasma oscillation frequency ωp. In our work, the electron density is $3.9 \times 10^{19}\, 1/{\mathrm{m}}^3$, which corresponds to the plasma frequency $\omega _{\mathrm{p}} = 3.5 \times 10^{12}\, {\mathrm{rad}}/{\mathrm{s}}$, the electron thermal temperature of 27.3keV with an external magnetic field of 8.01T.

    The dissipative surface states in our work are obtained by two methods: (1) solving the plane wave solution of Eq. (2) followed by matching the boundary condition and (2) numerical simulation with COMSOL, a commercially available software. In the first method, extra boundary condition demanding the continuity of the surface normal direction's electric field is added. This extra condition ensures the plasma will not spill out into the surrounding medium. In COMSOL, the surface states are solved using the generalized PDE module, and the results can be found in Supplementary Fig. S4 and Supplementary Note 6 in the Supplementary Information.

Acknowledgements
  • This work was partially supported by the National Natural Science Foundation of China (Grant Nos. 11874269 and 61875133). SZ acknowledges financial support from H2020 European Research Council Project Nos. 734578 (D-SPA), 777714 (NOCTURNO) and 648783 (TOPOLOGICAL), the Royal Society, the Wolfson Foundation.

Author contributions
  • W.W., W.G., Y.X., and S.Z. conceived the initial idea; W.W. and W.G. performed the analytical and the numerical analysis; Y.X., S.Z. directed the project. All authors discussed the results and commented in the manuscript.

Conflict of interest
  • The authors declare that they have no conflict of interest.

Supplementary information
Reference (40)

Catalog

    /

    DownLoad:  Full-Size Img PowerPoint
    Return
    Return