Verallgemeinerung der Methode der singulären Normallösungen für isotherme irreversible Vlasov-Systeme

Irreversible Vlasov systems, i.e. systems governed by a Vlasov-type kinetic equation including entropy-producing collision terms, are treated by the techniques of singular normal modes and singular integral equations using a new indirect method which renders possible a straightforward generalization of the Case formalism as developed originally for collision-free Vlasov plasmas. This method is in contrast to a more complex method given by the present authors for the first application of the singular normal mode expansion to irreversible Vlasov systems (1970). The linearized Vlasov operator supplemented by complete Bhatnagar-Gross-Krook collision integrals as the most important model collision terms is analyzed in detail for a nonrelativistic, nondegenerate, stationary electron gas with neutralizing positive ions and neutral particles without a magnetic field at constant temperature, generalizations for more complex irreversible Vlasov systems being possible. The key of the indirect method given is the introduction of a transformed electron distribution function containing as an additive term an integral over the usual distribution function.     

Operatorenmethode zur Lösung der Vlasov-Gleichung mit BGK-Stoßglied und äußerem elektrischem Feld

Irreversible processes in a classical electron plasma are treated on the basis of a linearized Vlasov equation supplemented by Bhatnagar-Gross-Krook terms describing electron-electron and electron-ion collisions correctly. The infinitely extended plasma is under the action of a space- and time-dependent external electric field. A general method of solution with projection operator techniques is given which results in a system of two coupled Volterra integral equations of the convolution type for the internal electric field and the current density. From there follows the electron distribution function, the electric field in the plasma, the electrical conductivity and a very general dispersion relation including Landau and collision damping. The method given can be generalized f. i. for multicomponent plasmas and for strong external electric fields.     

Thermalization of nonequilibrium electrons in quantum wires

We study the problem of energy relaxation in a one-dimensional electron system. The leading thermalization mechanism is due to three-particle collisions. We show that for the case of spinless electrons in a single channel quantum wire the corresponding collision integral can be transformed into an exactly solvable problem. The latter is known as the Schr?dinger equation for a quantum particle moving in a P?schl-Teller potential. The spectrum for the resulting eigenvalue problem allows for bound-state solutions, which can be identified with the zero modes of the collision integral, and a continuum of propagating modes, which are separated by a gap from the bound states. The inverse gap gives the time scale at which counterpropagating electrons thermalize.     

Scattering of long-wavelength acoustic phonons by isotopic impurities

We consider the problem of the relaxation of an arbitrary initial distribution function of a gas of long-wave acoustic phonons scattered by isotopic impurities embedded in a crystalline medium with cubic symmetry. The spectral decomposition of the collision integral of the suitable Boltzmann-Peierls equation is obtained. The spectrum of the collision operator is purely discrete and in addition to the eigenvalue 0 consists of three other eigenvalues. Explicit analytic expressions for these eigenvalues are obtained. Within the Chapman-Enskog approximation we derive the diffusion equation for the density of phonons and obtain the explicit expression for the diffusion coefficient. The dependency of the eigenvalues of the collision operator and the diffusion coefficient on the elastic constants of the medium is studied.     

The boundary value problem in fermion systems

The half-space boundary value problem for fermions near zero temperature in plane geometry is solved for diffuse boundary scattering by numerically constructing the spatial propagator in terms of the eigenfunctions of a generalized eigenvalue problem for the linearized Uehling-Uhlenbeck collision integral. The slip length is calculated for several interparticle scattering laws and compared with a relaxation time ansatz result and the experimental values for normal fluid³He. It is shown that the nonsingular part of the collision operator is relatively compact to the singular part.     

Normal modes of a relativistic quantum plasma; The one-component plasma

The normal modes of a relativistic electron gas are studied on the basis of the Boltzmann-Vlasov kinetic equation via a projection operator formalism. A general framework is constructed in which the fully relativistic Vlasov self-consistent force term appears as a symmetric operator acting in the Hilbert space of one-particle states. The plasma-dynamical equations are obtained by projecting onto the subspace consisting of the charge, energy and momentum densities, plus the nonconserved current density. The eigenmodes of these equations include two transverse and two longitudinal plasma modes, and one damped heat mode. They are explicitly calculated up to second order in the wave vector and to first order in the collision frequency.     

Diffusion in inhomogeneous media

The problem is to establish the correct diffusion equation in a medium that is inhomogeneous and whose temperature also varies in space. As a special model we study particles whose phase space distribution obeys Kramers' equation with a generalized collision operator. In the usual limit of strong collisions a diffusion equation is obtained. This equation contains additional drift terms, which depend on the form of the collision operator. They cannot be expressed as a mobility and a diffusion coefficient, unless the decay law of the velocity happens to be linear. Conclusion: no universal form of the diffusion equation exists, but each system has to be studied individually.Dedicated to Professor Harry Thomas on the occasion of his 60th birthday     

Microscopic theory of sound propagation in dielectric crystals

The time evolution of the atomic displacement field in a dielectric crystal subjected to an external force is studied in the domain of linear response by means of imaginary time Green's functions. For slowly varying disturbances two coupled equations have to be solved: a differential equation for the amplitude of an acoustic wave and a linearized Boltzmann equation. The latter results from the integral equation for the vertex part and includes an additional integral operator. The collision equation is solved for different relative magnitudes of the sound frequency and the frequencies for normal and Umklapp processes using the method developed by Weiss. Some of the expressions showing up in the velocity and damping of the sound wave are estimated numerically for rare gases with two-body forces in the form of the Morse potential.     

A stochastic derivation of the multivariate master equation describing reaction diffusion systems

We derive the multivariate master equation describing reaction diffusion systems from a discrete form master equation in phase space, assuming that the elastic collisions of the chemically active substances with the inert carrier gas have relaxed. In this state of collisional equilibrium the stochastic operator modelling the displacement of the particles between spatial cells reduces to the random wall operator and the reactive collision term yields the usual birth and death operator. Correlation functions are derived and their validity is discussed.     

托卡马克中电子温度梯度不稳定性的数值研究

用本征模积分方程研究了环形等离子体中的电子温度梯度不稳定性，在接近不稳定性阈值的参数区域内对模和输运的特征进行了研究。文中引入了一个新的积分参数，将实平面的积分解析延拓到复平面。这样可以同时对增长模和阻尼模进行研究。对被积函数的奇点进行了处理，给出了临界梯度的拟合公式，与有关的实验进行了比较，计算结果接近实验测量值。     

Remarks on the Calculation of Eigenvalues and Generalized Eigenfunctions Using a Finite Dimensional Approximation for the Inscattering Part of the Boltzmann Collision Operator

From several points of view it is of advantage to know the properties of the collision operators in kinetic equations, e.g. in the well known Boltzmann equation, in particular for the purpose of solving eigenvalue problems. With regard to elastic, exciting and deexciting processes some attemps were recently made to investigate such operators in the Boltzmann equation describing the behaviour of electrons in weakly ionized plasmas. In the following we will prove that in the case of a finite dimensional inscattering operator the eigenvalues and the corresponding eigenfunctions can be represented explicitly. Finite dimensional operators were used successfully in special models to approximate the inscattering operators; they possess the property of compactness and are well suitable for analytical or numerical calculations. The representation has been obtained by solving an adequate linear equation system. The generalized eigenfunctions correspond to the normal solutions used by Case in the neutron transport theory. The regularization of the singular integrals which are necessary to obtain this solution will be given in detail. Further a velocity dependence of the collision frequency which need not be monotonous in the considered case and the dependence on the direction could be included.     

Initialized fractional differential equations with Riemann-Liouville fractional-order derivative

The initial value problem of fractional differential equations and its solving method are studied in this paper. Firstly, for easy understanding, a different version of the initialized operator theory is presented for Riemann-Liouville’s fractional-order derivative, addressing the initial history in a straightforward form. Then, the initial value problem of a single-term fractional differential equation is converted to an equivalent integral equation, a form that is easy for both theoretical and numerical analysis, and two illustrative examples are given for checking the correctness of the integral equation. Finally, the counter-example proposed in a recent paper, which claims that the initialized operator theory results in wrong solution of a fractional differential equation, is checked again carefully. It is found that solving the equivalent integral equation gives the exact solution, and the reason behind the result of the counter-example is that the calculation therein is based on the conventional Laplace transform for fractional-order derivative, not on the initialized operator theory. The counter-example can be served as a physical model of creep phenomena for some viscoelastic materials, and it is found that it fits experimental curves well.     

The numerical method for three-dimensional impact with friction of multi-rigid-body system

The differential equations for planar impacts reduce to an algebraic form, and can be easily solved. For three dimensional impacts with friction, there is no closed-form solution, and numerical integration is required due to the swerve behavior of tangential impulse during collisions. The dynamic governing equations in the impact process are built up in impulse space based on the Lagrangian equation in this paper. The coefficient of restitution defined by Poisson is used as the condition of impact termination. A valid numerical method for solving three-dimensional frictional impact of multi-rigid body system is established. The singular cases of tangential movement in sticking point are especially noticed and analyzed. Several examples are present to reveal the different kinds of tangential movement modes varied with the normal impulse during collision.     

Complex Padé approximant operators for wide-angle beam propagation

The conventional rational Hadley(m, n) approximant of wide-angle beam propagator based on real Padé approximant operators incorrectly propagates the evanescent modes. In order to overcome this problem, two complex Padé approximants of wide-angle beam propagator are presented in this paper. The complex propagators of the first approach are obtained by using the same recurrence formula from the scalar Helmholtz equation of the conventional approximant method with a different initial value while those of the second method derived from Hadley(m, n) approximant of a square-root operator that has been rotated in the complex plane. These resulting approaches allow more accurate approximations to the Helmholtz equation than the well-known real Padé approximant. Furthermore, our proposed complex Padé approximant operators give the evanescent modes the desired damping.     

The Boltzmann Equation for Bose-Einstein Particles: Condensation in Finite Time

The paper considers the problem of the Bose-Einstein condensation in finite time for isotropic distributional solutions of the spatially homogeneous Boltzmann equation for Bose-Einstein particles with the hard sphere model. We prove that if the initial datum of a solution is a function which is singular enough near the origin (the zero-point of particle energy) but still Lebesgue integrable (so that there is no condensation at the initial time), then the condensation continuously starts to occur from the initial time to every later time. The proof is based on a convex positivity of the cubic collision integral and some properties of a certain Lebesgue derivatives of distributional solutions at the origin. As applications we also study a special type of solutions and present a relation between the conservation of mass and the condensation.     

Microlocal Analysis for Waves Propagating in Einstein & de Sitter Spacetime

We consider the waves propagating in the Einstein & de Sitter spacetime, which obey the covariant d’Alembert’s equation. That equation has singular coefficients and belongs to the family of the non-Fuchsian partial differential operators. We introduce the initial value problem for this equation and give the parametrices in the terms of Fourier integral operators. We also discuss the propagation and reflection of the singularities phenomena.     

Collision dynamics of two adjacent oscillators

A linear analysis is made of a single collision between two single-degree-of-freedom systems separated by a gap. The contact is modelled by a spring and a viscous damper. The approach is to describe the motion of the pair as being composed of sum and difference displacements. The equation of motion during contact is found and the solution is obtained from the conditions at initial contact. The main parameters are the ratio of strain energy to kinetic energy at initial contact, and the damping of the contact. The contact time and the energy loss are calculated, which gives an expression for the coefficient of restitution for the collision. This coefficient is found to be dependent on the collision velocity, but becomes constant for strong collisions.     

Collisional gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

Applying gyrokinetic simulations in theoretical turbulence and transport studies for the plasma edge and scrape-off layer (SOL) presents significant challenges. To particularly account for steep density and temperature gradients in the SOL, the “full-f” code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code, is based on an electrostatic model with a linearized field equation, and uses kinetic electrons. In previously published results, we applied PICLS to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) in the SOL without collisions. As an extension to this collision-less case and in preparation for 3D simulations, in this work, a collisional model will be introduced. The implemented Lenard–Bernstein collision operator and its Langevin discretization will be shown. Conservation properties of the collision operator, as well as a comparison of the collisional and non-collisional case, will be discussed.     

Bulk viscosity of low-temperature strongly interacting matter

We study the bulk viscosity of a pion gas in unitarized Chiral Perturbation Theory at low and moderate temperatures, below any phase transition to a quark-gluon plasma phase.We argue that inelastic processes are irrelevant and exponentially suppressed at low temperatures. Since the system falls out of chemical equilibrium upon expansion, a pion chemical potential must be introduced, so we extend the existing theories that include it. We control the zero modes of the collision operator and Landau?s conditions of fit when solving the Boltzmann equation with the elastic collision kernel.The dependence of the bulk viscosity with temperature is reminiscent of the findings of Fernández-Fraile and Gómez Nicola (2009) [1], while the numerical value is closer to that of Davesne (1996) [2]. In the zero-temperature limit we correctly recover the vanishing viscosity associated to a non-relativistic monoatomic gas.