Computer Simulation of Plasmas

Представлен новый метод решения релятивистских уравнений движения заряженных частиц в электромагнитных полях, учитывающий условие постоянства их значений на каждом временном шаге. Проведено сравнение точности и эффективности вычислений при решении тестовых задач в дву- и трехмерной постановках на основе нового метода, метода Бориса и его модификаций. В каждом случае рассматривались варианты аналитически и дискретно заданных значений электрического и магнитного полей. At present, the Boris method is mostly common for the numerical solution of problems of the dynamics of charged particles in electromagnetic fields. In recent years, new modifications of the Boris method have appeared that are capable to simplify, refine, or speed up calculations. However, despite the rather large number of publications which consider algorithms for calculating the trajectories of charged particles, new, more accurate and high-performance algorithms are desired. The article presents a detailed analysis of some explicit methods for solving the problem of the motion of charged particles in electromagnetic fields, which differ in how they set the average values of velocities. A new scheme based on an analytical solution, which has not been studied before, is proposed. For the Boris, Vay, Higuera – Cary schemes and the new VD1 scheme, a comparative analysis of the accuracy, convergence and counting time for the relativistic and nonrelativistic cases is carried out. In the nonrelativistic case, the new scheme allows accurate solving of the equations for motion of charged particles, and in the relativistic case, the accuracy of the new scheme remains higher compared to other schemes. As the relativistic factor increases, the accuracy of all considered schemes for calculating the trajectory of particles in an electromagnetic field decreases. For large values of the relativistic parameter, the new scheme is more accurate than the other schemes and retains the second order. In practical calculations, when particles can have different speeds, the scheme automatically adjusts to the given speed of the particle. The presented new scheme for calculating trajectories can be important for solving a wide range of problems in astrophysics and thermonuclear fusion, where high accuracy in determining particle trajectories in inhomogeneous fields for large moments of time is required.

On progressive blast envelope evolution of charged particles in electromagnetic fields

In this chapter we consider some problems in solid-state physics and in the theory of motion of charged particles in electromagnetic fields, whose solutions lead to incomplete cylindrical functions.

The Wigner formalism is a convenient way of describing quantum mechanical effects through a framework of distribution functions in phase space. Currently, there are stochastic and deterministic approaches in use. In our deterministic method, the critical discretization of the diffusion term is done through the utilization of an integral formulation of the Wigner equation. This deterministic method is studied in the context of superposed quantum states as a precursor to simulations of entangled states.

A nonseparable wave-like integro-differential equation for the time evolution of the Wigner distribution function in phase space is educed from the corresponding separable kinetic equation. By employing the quantum hydrodynamical description, a non-local evolution wave equation is also derived by synthesizing the Hamilton‐Jacobi equation with that of continuity, which predicts the generation of nonlocal and quadrupole quantum phenomena in the propagation of the spatial probability density.

A new scheme for the calculation of the collision integral of the nonlinear Boltzmann kinetic equation is presented, based on a piecewise constant discrete approximation of the distribution function in phase space and corresponding precalculated collision rate coefficients. The computing effort required for evaluation of the fivefold collision integral is significantly reduced by the use of precalculated rate coefficients. The method was combined with the conservative splitting method to solve the Boltzmann equation for a number of test cases to examine the performance of this method and demonstrate its efficiency. 10 refs., 6 figs.

Analytical solution for self-consistent particle equilibrium distribution in RF field with transverse focusing is found. Solution is attained in approximation of high brightness beam. Distribution function in phase space is determined as a stationary function of the energy integral. Equipartitioning for beam distribution between degrees of freedom follows directly from the choice of stationary distribution function. Analytical expressions for r-z equilibrium beam profile and maximum beam current in RF field are obtained.

The Rouse model for dilute polymer solutions undergoing homogeneous flows has been generalized to include the inertia of the beads in the equations of motion. To obtain the correct ‘‘diffusion equation’’ for the probability density distribution function in phase space, we generalize the diffusion equation derived by Murphy and Aguirre [J. Chem. Phys. 57, 2098 (1972)] from Hamilton’s equations of motion for an arbitrary number of interacting Brownian particles at equilibrium. Material functions are found, and the noninertial case is seen to be obtained as the zero mass limit in all steps of the development. In particular, the steady-state shear results are unaffected by the inclusion of inertia. It is also shown how two assumptions, ‘‘equilibration in momentum space,’’ and the neglect of acceleration, made independently by Curtiss, Bird, and Hassager in their phase-space kinetic theory, are actually the result of assuming zero mass.

Self-accelerating beam is a kind of light beam capable of self-bending in free space without any external potential, of which a typical one is the well-known Airy beam. Such a beam has gained great attention for its extraordinary properties, including nondiffracting, self-accelerating and self-healing, which may have versatile applications in the delivery and guiding of energy, information and objects using light, such as particle manipulation, micro-machining, optical routing, super-resolution imaging, etc. However, since Airy beam can only propagate along parabolic trajectory, which reduces the flexibility in practical applications, thus how to design accelerating beams propagating along arbitrary trajectory is still a crucial problem in this area. One scheme is to keep on finding other analytical solutions of the wave equation besides Airy beam, such as semi-Bessel accelerating beams, Mathius beams, and Weber beams, moving along circular, elliptical, or parabolic trajectories, but it becomes increasingly difficult to find out any more solutions. A more effective solution to this problem is based on the caustic method, which associates the predesigned trajectory with an optical caustics and then obtains the necessary initial field distribution by performing a light-ray analysis of the caustics. This method has been implemented in real space and Fourier space based on Fresnel diffraction integral and angular-spectrum integral, respectively. It has been found recently that they can be unified by constructing Wigner distribution function in phase space. Based on the caustic method, accelerating beams were constructed to propagate along arbitrary convex trajectories in two-dimensional space at first. With continuous development of this method, the types of accelerating beams available have been extending from convex trajectories to nonconvex trajectories, from two-dimensional trajectories to three-dimensional trajectories, and from one main lobe to multiple main lobes, which opens up more possibilities for emerging applications based on accelerating beams. In future, previous researches and applications based on Airy beams will certainly be generalized to all these new types of accelerating beams, and owing to the great flexibility in designing accelerating beams, more application scenarios may emerge in this process with huge development potential. Thus in this paper, we review the principle and progress of the caustic method in designing accelerating beams.

A semiclassical independent electron model for ion-atom collisions

Landau theory of Fermi liquids and nuclear fluid dynamics

The entropy of an electronic system is defined in terms of the Husimi function, a nonnegative distribution function in phase space. The Husimi function is calculated by maximizing the entropy subject to the constraints that the Husimi function give a Gaussian convolution of the desity when integrating over the momentum coordinates and that its second moment with respect to momentum give a sum of Gaussian convolutions of the density and the kinetic energy density. The result is compared with the Wigner function. Equations are given for calculating the density matrix from the Husimi function. The resulting equation for the exchange energy requires a difficult numerical integration. An alternate method is used to obtain the density matrix from an approximate partially collapsed Husimi matrix that gives the maximum entropy Husimi function as its diagonal. The results are exact for the harmonic oscillator ground state. Exchange energies calculated for H and the He isoelectronic series through C+4 show slight improvements over those calculated using a maximum entropy Wigner function.

This note is concerned with an optimal control problem governed by the relativistic Maxwell--Newton--Lorentz equations, which describe the motion of charged particles in electro-magnetic fields and consists of a hyperbolic PDE system coupled with a nonlinear ODE. An external magnetic field acts as control variable. Additional control constraints are incorporated by introducing a scalar magnetic potential which leads to an additional state equation in the form of a very weak elliptic PDE. Existence and uniqueness for the state equation is shown and the existence of a global optimal control is established. Moreover, first-order necessary optimality conditions in the form of Karush--Kuhn--Tucker conditions are derived. A numerical test illustrates the theoretical findings.

Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model

Charged Particles in Electromagnetic Fields