EIGENVALUE INTEGRAL RELATION AND STABILITY CONDITION FOR A MODE

In this paper,we point out that for a static system the eigenvalue integral relation(iσ)~2I +iσΦ+Σ=0 derived from an eigenequation (?)Ψ=[(iσ)~2(?)_2+(iσ)(?)_1+(?)_0]Ψ(r)=0 only gives the suffi- ciency of condition Σ>0 for the stability of a certain mode,but it can not provide the necessity to the condition in general.Three theorems presented in this paper show that under some conditions, Σ>0 makes the mode stable,while Σ<0 makes the mode unstable.

Determination of the Range of Magnetic Interactions from the Relations between Magnon Eigenvalues at High-Symmetry k Points

Magnetic exchange interactions(MEIs) define networks of coupled magnetic moments and lead to a surprisingly rich variety of their magnetic properties. Typically MEIs can be estimated by fitting experimental results.Unfortunately, how many MEIs need to be included in the fitting process for a material is unclear a priori,which limits the results obtained by these conventional methods. Based on linear spin-wave theory but without performing matrix diagonalization, we show that for a general quadratic spin Hamiltonian, there is a simple relation between the Fourier transform of MEIs and the sum of square of magnon energies(SSME). We further show that according to the real-space distance range within which MEIs are considered relevant, one can obtain the corresponding relationships between SSME in momentum space. By directly utilizing these characteristics and the experimental magnon energies at only a few high-symmetry k points in the Brillouin zone, one can obtain strong constraints about the range of exchange path beyond which MEIs can be safely neglected. Our methodology is also generally applicable for other Hamiltonian with quadratic Fermi or Boson operators.

PDRK:A General Kinetic Dispersion Relation Solver for Magnetized Plasma

A general,fast,and effective approach is developed for numerical calculation of kinetic plasma linear dispersion relations.The plasma dispersion function is approximated by J-pole expansion.Subsequently,the dispersion relation is transformed to a standard matrix eigenvalue problem of an equivalent linear system.Numerical solutions for the least damped or fastest growing modes using an 8-pole expansion are generally accurate;more strongly damped modes are less accurate,but are less likely to be of physical interest.In contrast to conventional approaches,such as Newton＇s iterative method,this approach can give either all the solutions in the system or a few solutions around the initial guess.It is also free from convergence problems.The approach is demonstrated for electrostatic dispersion equations with one-dimensional and twodimensional wavevectors,and for electromagnetic kinetic magnetized plasma dispersion relation for bi-Maxwellian distribution with relative parallel velocity flows between species.