Weakly Collisional Landau Damping and BGK Modes New

Weakly Collisional Landau Damping and BGK Modes: New Results on Old Problems A. Bhattacharjee, C. S. Ng, and F. Skiff Space Science Center University of New Hampshire Davidson Symposium, Princeton, June 11 -12, 2007

What have I learned from Ron? 1. 2. 3. 4. Be broad in your perspectives---plasma physics is diverse, and yet deeply interconnected. Learn plasma theory by doing it----calculate, and do so as rigorously as possible. Don’t worry if the calculation is long---if you are careful, the terms will all cancel out in the end leaving a nice result. If you want to serve your community and do your calculations, be disciplined and organized. Advice from one Editor to another: choose your reviewers well, and let the review process of a paper be a dialogue between the authors and reviewers.

Introduction High-temperature plasmas are nearly collisionless. Two classic results in a collisionless plasma: Landau damping of linear electrostatic waves Landau (1946) Exact nonlinear 1 D electrostatic solutions (BGK modes) ---- (undamped) Bernstein, Greene & Kruskal (1957) This talk: two new results on these old and classic problems.

Vlasov-Poisson equations

Landau Damping of Plasma Oscillations • Vlasov (1938): Collisionless kinetic theory, asssuming normal modes of the form • Landau (1946): “…most of his (A. A. Vlasov’s) results turn out to be incorrect. Vlasov looked for solutions of the form and determined the dependency of the frequency on the wave vector k. Actually there exists no dependence of on k at all, and for given value of k, arbitrary values of are possible. ” • � Landau solved the initial-value problem (using Laplace transforms) and obtained collisionless damping for monotonic distribution functions. Landau-damped solutions are not eigenmodes, but represent the linear response of the system in the asymptotic limit

Vexing question For a linear perturbation of the form if the distribution function is non-monotonic, there are unstable eignemodes with. In this case Landau’s analysis coincides with Vlasov’s. But if the distribution function is monotonic (such as a Maxwellian), the solutions with are not eigenmodes. What is the physical reason for this strange asymmetry?

Case-Van Kampen modes [Van Kampen (1955), Case (1959)] Vlasov-Poisson equation for a single Fourier mode: Exact solution (singular and undamped): These are the Case-Van Kampen eigenmodes. Landau-damped modes are not eigenmodes, but long-time remnants of an arbitrary smooth initial condition.

Case-Van Kampen modes form a complete set A general solution : Phase-mixing of undamped Case-Van Kampen modes produces Landau damping of the plasma wave. Thus, the linear Vlasov-Poisson is completely solved. Question: what is the effect of collision, even if it is weak? A model equation : Lenard & Bernstein (1958)

Collision as a singular perturbation --- normalized collision frequency Collisions have to be included if there are sharp gradients in velocity space. Each Case-Kampen mode has an infinitely sharp gradient in velocity space: Eigenmodes of the system change greatly even if collisions are extremely weak.

Discrete eigenmodes in experiments Skiff et al. (1998) measured distribution functions with great accuracy in a weakly collisional plasma using laser-induced florescence Eigenmodes decay exponentially, instead of as expected from classic theory [Su & Oberman (1968)] or

A complete set of discrete eigenmodes For an initial value problem: Ng et al. , PRL, (1999, 2004) gn found in closed form involving incomplete gamma function. cn can be found by integration involving the initial data.

Properties of the complete set of eigenmodes The new set replaces the Case-Van Kampen modes. It is discrete, unlike the Case-Van Kampen continuous spectrum Dispersion relation tends to Landau’s roots in the limit of weak collisions: except near the collisional modes All eigenfunctions are non-singular, unlike Case-Van Kampen modes.

Calculating the eigenvalues

Shape of the eigenfunctions

What is a BGK mode? An exact undamped nonlinear solution of the steady -state Vlasov-Poisson system of equations. Normalized; uniform ion background (for simplicity) 1 D solution: Bernstein, Greene & Kruskal (1957) Traveling wave solution by change of reference frame: v v x x

Construction method of 1 D BGK mode BGK (1957) 1 st equation can be solved by: Integro-differential equation: Can be solved by given f(w) or (x).

Physical picture of 1 D BGK mode Electron velocity increases in the center, so electron density decreases

3 D features in BGK mode observations Ergun et al. (1998) Traveling direction bipolar B unipolar

3 D BGK mode for finite B? For infinitely strong B: Chen & Parks (2002) basic assumption: electrons moving along B only back to 1 -D problem Conjecture: No solution for B = 0? [Chen 2002] Proof: conjecture is true if f = f(w), i. e. , f depends only on energy [Ng & Bhattacharjee 2005]

No 2 D/3 D solution if f depends only on w - 1 D: - 2 D: + + + + - 3 D: - - In 2 D/3 D, -

3 D solution depending on energy as well as angular momentum - For a spherical potential with solution of the Vlasov equation. Possible to support the electric field self-consistently --- satisfying the Vlasov-Poisson equation. - - is also a + + + + - - -

Conclusion In Landau damping, even weak collisions can have profound implications, and change completely the nature of the spectrum of eigenmodes. A new complete spectrum of discrete eigenmodes is found that replaces the Case-Van Kampen continuous spectrum, where Landau solutions now become the true eigenmodes. 2 D/3 D BGK modes cannot exist if the distribution function depends only on energy. 3 D BGK modes for B=0, and 2 D BGK modes for finite B are constructed when f also depends on angular momentum. See following papers for more details: C. S. Ng, A. Bhattacharjee, and F. Skiff, Phys. Rev. Lett. 83, 1974 (1999). C. S. Ng, A. Bhattacharjee, and F. Skiff, Phys. Rev. Lett. 92, 065002 (2004). C. S. Ng and A. Bhattacharjee, Phys. Rev. Lett. 95, 245004 (2005).