Ultrafast dynamics of a spin-polarized electron plasma with magnetic ions

TL;DR

This paper develops a mean-field model and computational methods to study the ultrafast nonlinear dynamics of spin-polarized electron gases interacting with magnetic ions, revealing stability conditions and instabilities relevant to dense plasma regimes.

Contribution

It introduces a coupled Vlasov-Poisson-Landau-Lifshitz model for spin-electron-ion interactions and develops a high-accuracy Eulerian simulation approach for nonlinear dynamics.

Findings

Identifies conditions for stability and instability of spin modes.

Validates linear response predictions with numerical simulations.

Explores charge and spin instabilities in two-stream equilibria.

Abstract

We construct a mean-field model that describes the nonlinear dynamics of a spin-polarized electron gas interacting with fixed, positively-charged ions possessing a magnetic moment that evolves in time. The mobile electrons are modeled by a four-component distribution function in the two-dimensional phase space $(x,v)$, obeying a Vlasov-Poisson set of equations. The ions are modeled by a Landau-Lifshitz equation for their spin density, which contains ion-ion and electron-ion magnetic exchange terms. We perform a linear response study of the coupled Vlasov-Poisson-Landau-Lifshitz (VPLL) equations for the case of a Maxwell-Boltzmann equilibrium, focussing in particular on the spin dispersion relation. Condition of stability or instability for the spin modes are identified, which depend essentially on the electron spin polarization rate $\eta$ and the electron-ion magnetic coupling constant $K$. We also develop an Eulerian grid-based computational code for the fully nonlinear VPLL equations, based on the geometric Hamiltonian method first developed in [N. Crouseilles et al. Journal of Plasma Physics, 89(2):905890215, 2023]. This technique allows us to achieve great accuracy for the conserved quantities, such as the modulus of the ion spin vector and the total energy. Numerical tests in the linear regime are in accordance with the estimations of the linear response theory. For two-stream equilibria, we study the interplay of instabilities occurring in both the charge and the spin sectors. The set of parameters used in the simulations, with densities close to those of solids ($\approx 10^{29} \rm m^{-3}$) and temperatures of the order of 10 eV, may be relevant to the warm dense matter regime appearing in some inertial fusion experiments.