Plasmas have been recently studied as topological materials. However, a comprehensive picture of topological phases and topological phase transitions in cold magnetized plasmas is still missing. Here we systematically map out all the topological phas
es and establish the bulk-edge correspondence in cold magnetized plasmas. We find that for the linear eigenmodes, there are 10 topological phases in the parameter space of density $n$, magnetic field $B$, and parallel wavenumber $k_{z}$, separated by the surfaces of Langmuir wave-L wave resonance, Langmuir wave-cyclotron wave resonance, and zero magnetic field. For fixed $B$ and $k_{z}$, only the phase transition at the Langmuir wave-cyclotron wave resonance corresponds to edge modes. A sufficient and necessary condition for the existence of this type of edge modes is given and verified by numerical solutions. We demonstrate that edge modes exist not only on a plasma-vacuum interface but also on more general plasma-plasma interfaces. This finding broadens the possible applications of these exotic excitations in space and laboratory plasmas.
We develop an Explicitly Solvable Energy-Conserving (ESEC) algorithm for the Stochastic Differential Equation (SDE) describing the pitch-angle scattering process in magnetized plasmas. The Cayley transform is used to calculate both the deterministic
gyromotion and stochastic scattering, affording the algorithm to be explicitly solvable and exactly energy conserving. An unusual property of the SDE for pitch-angle scattering is that its coefficients diverge at the zero velocity and do not satisfy the global Lipschitz condition. Consequently, when standard numerical methods, such as the Euler-Maruyama (EM), are applied, numerical convergence is difficult to establish. For the proposed ESEC algorithm, its energy-preserving property enables us to overcome this obstacle. We rigorously prove that the ESEC algorithm is order 1/2 strongly convergent. This result is confirmed by detailed numerical studies. For the case of pitch-angle scattering in a magnetized plasma with a constant magnetic field, the numerical solution is benchmarked against the analytical solution, and excellent agreements are found.
A method of Parity-Time (PT)-symmetry analysis is introduced to study the high dimensional, complicated parameter space of drift wave instabilities. We show that spontaneous PT-symmetry breaking leads to the Ion Temperature Gradient (ITG) instability
of drift waves, and the collisional instability is the result of explicit PT-symmetry breaking. A new unstable drift wave induced by finite collisionality is identified. It is also found that gradients of ion temperature and density can destabilize the ion cyclotron waves when PT symmetry is explicitly broken by a finite collisionality.
