Nonlinear evolution of a relativistically strong electromagnetic wave in self-created electron-positron plasma

The nonlinear interaction between the electron-positron pairs produced by an electromagnetic wave in plasma and the wave leads to damping of the wave, frequency upshift, change of polarization, and particle acceleration. The case of a circularly polarized wave is investigated in the framework of the relativistic Vlasov equation with a source term based on the Schwinger formula for the pair creation rate.     

A Hilbert-Vlasov code for the study of high-frequency plasmabeatwave accelerator

High-frequency beatwave simulations relevant to the University of California at Los Angeles (UCLA) experiment with relativistic eulerian hybrid Vlasov code are presented. These Hilbert-Masov simulations revealed a rich variety of phenomena associated with the fast particle dynamics induced by beatwave experiment for a high ratio of driver frequency to plasma frequency ω_pump/ω_plasma ≈33. The present model allows us to extend detailed modeling to frequency ratios greater than the current practical maximum of 10 or so, for Vlasov or particle-in-cell (PIC) codes, by replacing the Maxwell equations by mode equations for the electromagnetic Vlasov code. Numerical results, including beat frequency chirping (i.e., pump frequency linearly decreasing with time), show that the amplitude limit due to relativistic detuning can be enhanced with accelerated particles up to the ultrarelativistic energies with a high-acceleration gradient of more than 25 GeV/m     

Scattering of light waves by electron electrostatic waves in laser produced plasmas

The propagation of light waves in an underdense plasma is studied using one-dimensional Vlasov-Maxwell numerical simulation.It is found that the light waves can be scattered by electron plasma waves as well as other heavily and weakly damping electron wave modes,corresponding to stimulated Raman and Brilluoin-like scatterings.The stimulated electron acoustic wave scattering is also observed as a high scattering level.High frequency plasma wave scattering is also observed.These electron electrostatic wave modes are due to a non-thermal electron distribution produced by the wave-particle interactions.The collision effects on stimulated electron acoustic wave and the laser intensity effects on the scattering spectra are also investigated.     

Formation and dynamics of relativistic electromagnetic solitons in plasmas containing high-and low-energy electron components

We present theoretical and numerical studies of the nonlinear interactions between intense electromagnetic waves in plasmas containing high-and low-energy electron components. Such plasmas are frequently observed in laser-plasma experiments, where the hot electron component is created by the acceleration of electrons by strong electrostatic waves that are created by the laser-induced Raman forward and backward instabilities. The two-component electron plasma is described by the Vlasov equation for the hot electrons and the hydrodynamic equations for the cold electrons, which are coupled nonlinearly to the electromagnetic wave equation and the Poisson equation for the potential. The present nonlinear system is shown to admit electromagnetic solitary waves correlated with a positive potential and trapped electron islands from the hot electron population. The text was submitted by the authors in English.     

Production of nonisothermal electrons and Langmuir waves because of colliding ion holes and trapping of plasmons in an ion hole

We present new simulation studies exhibiting production of nonisothermal electron distributions and Langmuir waves by colliding ion holes and trapping of plasmons in an ion hole. We find that, during head-on ion hole collisions, streams of accelerated electrons are produced by the electrostatic potentials supporting the ion holes. Subsequently, Langmuir waves are excited by a two-stream instability involving energetic electron beams. The resulting Langmuir waves can be trapped in an ion hole. The present ion-hole-Langmuir wave interactions are unique kinetic phenomena which can be dealt with a Vlasov code, which we developed recently. The results can have relevance to the understanding of particle and field data that are forthcoming from different spacecraft missions in Earth's auroral ionosphere and the magnetosphere.     

激光等离子体二维相对论电磁粒子模拟

研制了激光等离子体二维相对论电磁粒子模拟程序(2DCIC)。追踪几万甚至百万个模拟粒子在外加激光场和自洽场中运动,自洽地计算电荷和电流密度,求解完全的Maxwel方程,电子的相对论运动方程和离子的牛顿运动方程,辅以灵活的诊断研究波-波,波-粒子相互作用的发生、发展和饱和的细节以及时间演化规律。激光可以正入射,也可以斜入射;等离子体可以是均匀密度,也可以具有密度梯度;为了节约机时,还发展了并行运算。物理模型参数可调,既适用于研究激光聚变等离子体相互作用,也适用于超短脉冲超强激光等离子体相互作用和其它等离子体问题。经过多次试算检验,对等离子体平衡态进行了计算研究,对于超短脉冲超强激光的传播也进行了初步模拟计算。     

Investigation of the particle spin properties on the velocity space instability in high‐density plasmas

The nonresonant electromagnetic instabilities of the anisotropic velocity space (Weibel‐like) have always been one of the interesting subjects for researchers. These electromagnetic instabilities play an important role in generating strong magnetic fields in laboratory plasmas for applications such as inertial confinement fusion and space plasmas. In this paper, we investigate the quantum effects of the particle spin on the electromagnetic instabilities. In the case of the presence of a magnetic dipole force and an electron precession frequency like the Vlasov equation, we derive the full quantum equation. This study shows that, in the presence of the spin‐polarized effects, the growth rate of the instabilities is reduced compared to the classical cases and will not arise for low fractions of the temperature anisotropy for different values of the magnetic field. Indeed, it is expected that the probability of electron capture in the background magnetic fields and the effective collision with the particle increase because of the spin effect, so that a high portion of the electron energy is transmitted to the background plasma, and the temperature anisotropy governing the electron distribution is reduced.     

共振吸收的自洽场结构和密度轮廓分布

 用二维多时标全电磁相对论粒子模拟程序，对斜入射P极化强激光与具有线性密度分布的非均匀等离子体平板相互作用中的共振吸收进行了模拟计算，给出了冕区自洽场的结构和等离子体密度轮廓分布，粒子模拟结果观察到了临界面附近等离子体波的激发、等离子体密度轮廓的变陡和空穴的产生。     

共振吸收的自洽场结构和密度轮廓分布

用二维多时标全电磁相对论粒子模拟程序,对斜入射P极化强激光与具有线性密度分布的非均匀等离子体平板相互作用中的共振吸收进行了模拟计算,给出了冕区自洽场的结构和等离子体密度轮廓分布,粒子模拟结果观察到了临界面附近等离子体波的激发、等离子体密度轮廓的变陡和空穴的产生。     

A conservative high order semi-Lagrangian WENO method for the Vlasov equation

In this paper, we propose a novel Vlasov solver based on a semi-Lagrangian method which combines Strang splitting in time with high order WENO (weighted essentially non-oscillatory) reconstruction in space. A key insight in this work is that the spatial interpolation matrices, used in the reconstruction process of a semi-Lagrangian approach to linear hyperbolic equations, can be factored into right and left flux matrices. It is the factoring of the interpolation matrices which makes it possible to apply the WENO methodology in the reconstruction used in the semi-Lagrangian update. The spatial WENO reconstruction developed for this method is conservative and updates point values of the solution. While the third, fifth, seventh and ninth order reconstructions are presented in this paper, the scheme can be extended to arbitrarily high order. WENO reconstruction is able to achieve high order accuracy in smooth parts of the solution while being able to capture sharp interfaces without introducing oscillations. Moreover, the CFL time step restriction of a regular finite difference or finite volume WENO scheme is removed in a semi-Lagrangian framework, allowing for a cheaper and more flexible numerical realization. The quality of the proposed method is demonstrated by applying the approach to basic test problems, such as linear advection and rigid body rotation, and to classical plasma problems, such as Landau damping and the two-stream instability. Even though the method is only second order accurate in time, our numerical results suggest the use of high order reconstruction is advantageous when considering the Vlasov–Poisson system.     

Numerical approximation of the Vlasov–Maxwell–Fokker–Planck system in two dimensions

A numerical method is developed for solving the Vlasov–Maxwell–Fokker–Planck system in two spatial dimensions. This system of equations is a model for a collisional plasma in the presence of a self consistent electromagnetic field. The numerical procedure is a type of deterministic particle method and is an extension to include the full electromagnetic field of the approximation method of Wollman and Ozizmir [S. Wollman, E. Ozizmir, Numerical approximation of the Vlasov–Poisson–Fokker–Planck system in two dimensions, J. Comput. Phys. 228 (2009) 6629–6669]. In addition, the long time asymptotic behavior of solutions is studied. It is determined that the solution to the Vlasov–Maxwell–Fokker–Planck system converges to the same steady state solution as that for the Vlasov–Poisson–Fokker–Planck system.     

Localization of intense electromagnetic waves in a relativistically hot plasma

We consider nonlinear interactions between intense short electromagnetic waves (EMWs) and a relativistically hot electron plasma that supports relativistic electron holes (REHs). It is shown that such EMW-REH interactions are governed by a coupled nonlinear system of equations composed of a nonlinear Schro dinger equation describing the dynamics of the EMWs and the Poisson-relativistic Vlasov system describing the dynamics of driven REHs. The present nonlinear system of equations admits both a linearly trapped discrete number of eigenmodes of the EMWs in a quasistationary REH and a modification of the REH by large-amplitude trapped EMWs. Computer simulations of the relativistic Vlasov and Maxwell-Poisson system of equations show complex interactions between REHs loaded with localized EMWs.     

Connecting shock velocities to electron-injection mechanisms

Electrons can be accelerated by their interaction with nonlinearly saturated electrostatic waves up to speeds with which they can undergo diffusive acceleration across supernova remnant shocks. Here, we model this wave-electron interaction by particle-in-cell and Vlasov simulations. We find that the lifetime of the saturated wave is considerably longer in the Vlasov simulation, due to differences in how these simulation methods approximate the plasma. Electron surfing acceleration which requires a stable saturated wave may thus be more important for electron acceleration at shocks than previously thought. For beam speeds above a critical value, which we estimate here, both simulation codes exclude surfing acceleration due to a rapid wave collapse.     

Novel features of non-linear Raman instability in a laser plasma

The electron phase space evolution in a non-relativistic and homogeneous laser plasma generated by a nanosecond laser in a near infrared region in the presence of stimulated Raman scattering is investigated by a numerical simulation. The mechanism of electron acceleration in the potential wells of the plasma wave accompanying the Raman back-scattering is analyzed in a 1D Vlasov-Maxwell model. The dominant wave modes are both the backward and the forward propagating Raman waves, each accompanied by a daughter electrostatic wave. In addition to a strong interaction of plasma electrons with the primary electrostatic wave in the case of back-scattering, a cascading is observed consisting in a secondary scattering of the primary Raman back-scattered wave. This phenomenon reduces the Raman reflectivity and causes an acceleration of electrons against the direction of the heating laser beam. Moreover, the strong trapping in the primary electrostatic wave generated by the Raman back-scattering leads due to the trapped particle instability to a significant spectral broadening of the original plasma wave and a subsequent intermittent behaviour of the scattering process. The high phase velocity electrostatic daughter wave of the forward Raman scattering cannot trap the electrons directly, but there is an indication of non-resonant quasi-modes combined of this wave and of the simultaneously existing electrostatic daughter wave accompanying the Raman back-scattering. The transform method is used for a solution of the set of partial differential equations, which consists of the Vlasov equation and of the full set of Maxwell equations in a 1D approximation. A simplified Fokker-Planck collision term is added to overcome the numerical instabilities during the simulation. The model has relevance to a long scale plasma geometry, such as occurring in the indirect drive experiments near the light entrance holes of target hohlraum.     

A Vlasov code for the numerical simulation of stimulated raman scattering

Numerical simulations of the stimulated Raman scattering are presented using an Eulerian relativistic Vlasov code. Such a code allows a finer resolution in phase space than a particle code and provides a better understanding of the acceleration process for the particles at relativistically high energy. Forward Raman scattering as well as backward Raman scattering are considered to illustrate the possibilities of the Eulerian Vlasov code.     

3cm相对论返波荡器的粒子模拟研究

 运用全电磁的相对论的25维粒子模拟程序MAGIC模拟分析了3cm相对论返波管中注波互作用的非线性过程，得到了器件微波输出功率，效率，工作频率等参量，模拟结果与实验数据较好地一致。对电子注特性参数，慢波结构尺寸，引导磁场进行了优化，获得了器件最佳化运行的工作参数。     

基于等离子体粒子模拟的喷气Z箍缩过程物理研究

给出了喷气Z箍缩动力学过程在二维柱坐标系下的等离子体粒子模拟物理模型,编写了相应的程序.对低电流驱动下的稀薄喷气Z箍缩动力学过程进行了验证性的等离子体粒子模拟,得到了许多微观的Z箍缩物理信息,如负载中的电流(密度)、电磁场、粒子位置和密度的时空演化,以及总的z箍缩拖尾质量和拖尾电流等信息.发现在Z箍缩过程中,模拟得到的等离子体电流随时间的变化反映出了等离子体箍缩到心和反弹的过程特征,磁场随径向的变化与长直导线电流给出的磁场很接近;电子所受到的电场力和磁场力(洛伦兹力)是相当的,而离子所受到的力主要是电场力;电子首先在z方向加速,然后在自身运动产生电流的磁场的作用下向轴心箍缩,而离子是在电子和离子电荷分离所产生的电场力的作用下向轴心运动;在压缩到轴心附近时,电子首先因静电排斥而飞散,而离子则在惯性的作用下继续向轴心箍缩,而后滞止飞散.Z箍缩等离子体的拖尾质量在20%左右,拖尾电流最大时在7%左右.     

量子电动力学—粒子模拟极端等离子体动力学研究

新一代拍瓦激光装置有望将激光强度提升至10²³~10²⁴ W·cm^-2,在此极端强场条件下非线性量子电动力学效应对等离子体动力学过程产生重要影响.相对论电子在强电磁场作用下会同步辐射大量伽马光子,当后者穿过超强电磁场时会级联产生正负电子对.与此同时,这些量子电动力学效应也会反作用于激光等离子体相互作用过程,如辐射阻尼严重影响电子运动过程.为了研究这样极端的等离子体动力学,我们介绍最近几年发展的量子电动力学数值模拟模块,并将其耦合到传统的粒子模拟程序中,即量子电动力学-粒子模拟程序.由于大量新辐射的光子和产生的正负电子对会造成模拟粒子数目的不断增加,我们发展了粒子融合技术来减小模拟规模.利用此量子电动力学-粒子模拟程序,我们对极端强场激光物质相互作用以及极端天体物理现象开展了数值模拟研究.     

电子回旋共振放电电离特性的PIC/MCC模拟

A theoretical and computational model is presented to study the ionization of the argon electron cyclotron resonance(ECR)microwave discharge using a quasi-three-dimensional electromagnetic particle-in- cell plus Monte Carlo collision method.The interaction between the charged particles and microwave fields are described by the electromagnetic mode of particle-in-cell method.The collision processes are treated with Monte Carlo method.The simulation code is the original work.The results of the particle simulation for the ECR discharge of argon gas which include the microscopic features of charged particles and the electromagnetic characteristics of the ECR discharge plasma,and also the transient phenomena have been presented.