Hamilton operator and the semiclassical limit for scalar particles in an electromagnetic field

We successively apply the generalized Case-Foldy-Feshbach-Villars (CFFV) and the Foldy-Wouthuysen (FW) transformation to derive the Hamiltonian for relativistic scalar particles in an electromagnetic field. In contrast to the original transformation, the generalized CFFV transformation contains an arbitrary parameter and can be performed for massless particles, which allows solving the problem of massless particles in an electromagnetic field. We show that the form of the Hamiltonian in the FW representation is independent of the arbitrarily chosen parameter. Compared with the classical Hamiltonian for point particles, this Hamiltonian contains quantum terms characterizing the quadrupole coupling of moving particles to the electric field and the electric and mixed polarizabilities. We obtain the quantum mechanical and semiclassical equations of motion of massive and massless particles in an electromagnetic field. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 156, No. 3, pp. 398–411, September, 2008.     

Formalism of two potentials for the numerical solution of Maxwell’s equations

A new formulation of Maxwell’s equations based on the introduction of two vector and two scalar potentials is proposed. As a result, the electromagnetic field equations are written as a hyperbolic system that contains, in contrast to the original Maxwell system, only evolution equations and does not involve equations in the form of differential constraints. This makes the new equations especially convenient for the numerical simulation of electromagnetic processes. Specifically, they can be solved by applying powerful modern shock-capturing methods based on the approximation of spatial derivatives by upwind differences. The cases of an electromagnetic field in a vacuum and an inhomogeneous material are considered. Examples are given in which electromagnetic wave propagation is simulated by solving the formulated system of equations with the help of modern high-order accurate schemes.     

Kinetic Theory of Quantum Electrodynamic Plasma in a Strong Electromagnetic Field: I. The Covariant Formalism

We present a covariant approach to the kinetic theory of quantum electrodynamic plasma in a strong electromagnetic field. The method is based on the relativistic von Neumann equation for the nonequilibrium statistical operator defined on spacelike hyperplanes in Minkowski space. We use the canonical quantization of the system on hyperplanes and a covariant generalization of the Coulomb gauge. The condensate mode associated with the mean electromagnetic field is separated from the photon degrees of freedom by a time-dependent unitary transformation of the dynamic variables and the nonequilibrium statistical operator. This allows using expansions of correlation functions and of the statistical operator in powers of the fine structure constant even in the presence of a strong electromagnetic field. We present a general scheme for deriving kinetic equations in the hyperplane formalism.     

New computational method for solving some 2-dimensional nonlinear Volterra integro-differential equations

The aim of this paper is to present an efficient numerical procedure for solving the two-dimensional nonlinear Volterra integro-differential equations (2-DNVIDE) by two-dimensional differential transform method (2-DDTM). The technique that we used is the differential transform method, which is based on Taylor series expansion. Using the differential transform, 2-DNVIDE can be transformed to algebraic equations, and the resulting algebraic equations are called iterative equations. New theorems for the transformation of integrals and partial differential equations are introduced and proved. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments.     

Furtivity and Masking Problems in Time-Dependent Electromagnetic Obstacle Scattering

In this paper, we consider furtivity and masking problems in time-dependent three-dimensional electromagnetic obstacle scattering. That is, we propose a criterion based on a merit function to minimize or to mask the electromagnetic field scattered by a bounded obstacle when hit by an incoming electromagnetic field and, with respect to this criterion, we drive the optimal strategy. These problems are natural generalizations to the context of electromagnetic scattering of the furtivity problem in time-dependent acoustic obstacle scattering presented in Ref. 1. We propose mathematical models of the furtivity and masking time-dependent three-dimensional electromagnetic scattering problems that consist in optimal control problems for systems of partial differential equations derived from the Maxwell equations. These control problems are approached using the Pontryagin maximum principle. We formulate the first-order optimality conditions for the control problems considered as exterior problems defined outside the obstacle for systems of partial differential equations. Moreover, the first-order optimality conditions derived are solved numerically with a highly parallelizable numerical method based on a perturbative series of the type considered in Refs. 2–3. Finally, we assess and validate the mathematical models and the numerical method proposed analyzing the numerical results obtained with a parallel implementation of the numerical method in several experiments on test problems. Impressive speedup factors are obtained executing the algorithms on a parallel machine when the number of processors used in the computation ranges between 1 and 100. Some virtual reality applications and some animations relative to the numerical experiments can be found in the website http://www.econ.unian.it/recchioni/w10/.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

Modified integral current methods in electrodynamics of nonhomogeneous media

The main difficulty in numerical solution of integral equations of electrodynamics is associated with the need to solve a high-order system of linear equations with a dense matrix. It is therefore relevant to develop numerical methods that lead to linear equation systems of lower order at the cost of more complex evaluation of the coefficients. In this article we propose a method for solving linear equations of electrodynamics which is a modification of the integral current method. The main distinctive feature of the proposed method is double integration of the electric Green’s tensor in the process of algebraization of the original integral equation. The solutions of the system of linear equations are thus integral means of the electric field inside the anomaly constructed by the proposed transformation formula. We prove convergence and derive error bounds for both the solution of the integral equation and the electromagnetic field components evaluated from approximate transformation formulas.     

High order moment closure for Vlasov-Maxwell equations

The extended magnetohydrodynamic models are derived based on the moment closure of the Vlasov-Maxwell (VM) equations. We adopt the Grad type moment expansion which was firstly proposed for the Boltzmann equation. A new regularization method for the Grad’s moment system was recently proposed to achieve the globally hyperbolicity so that the local well-posedness of the moment system is attained. For the VM equations, the moment expansion of the convection term is exactly the same as that in the Boltzmann equation, thus the new developed regularization applies. The moment expansion of the electromagnetic force term in the VM equations turns out to be a linear source term, which can preserve the conservative properties of the distribution function in the VM equations perfectly.     

Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives

In general, we will use the numerical differentiation when dealing with the differential equations. Thus the differential equations can be transformed into algebraic equations and then we can get the numerical solutions. But as we all have known, the numerical differentiation process is very sensitive to even a small level of errors. In contrast it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, based on the Sinc method we provide a new method using Sinc method incorporated with the double exponential transformation based on the interpolation of the highest derivatives (SIHD) for the differential equations. The error in the approximation of the solution is shown to converge at an exponential rate. The numerical results show that compared with the exiting results, our method is of high accuracy, of good convergence with little computational efforts. It is easy to treat nonhomogeneous mixed boundary condition for our method, which is unlike the traditional Sinc method.     

Designable Integrability of the Variable Coefficient Nonlinear Schrödinger Equations

The designable integrability (DI) [ 51 ] of the variable coefficient nonlinear Schrödinger equations (VCNLSEs) is first introduced by construction of an explicit transformation, which maps VCNLSE to the usual nonlinear Schrödinger equations (NLSEs). One novel feature of VCNLSE with DI is that its coefficients can be designed artificially and analytically by using transformation. A special example between nonautonomous NLSEs and NLSEs is given here. Further, the optical super‐lattice potentials (or periodic potentials) and multiwell potentials are designed, which are two kinds of important potential in Bose–Einstein condensation and nonlinear optical systems. There are two interesting features of the soliton of the VCNLSEs indicated by the analytic and exact formula. Specifically, its profile is variable and its trajectory is not a straight line when it evolves with time t.     

The first initial–boundary-value problem for a nonlinear model of the electrodynamics of vibrating elastic media

The first initial–boundary-value problem for nonlinear differential equations describing the interactions of a vibrating electroconductive body and the electromagnetic field is studied. We assume that the motion of the body occurs at velocities that are much smaller than the velocity of propagation of the electromagnetic waves through the elastic medium. The model under study consists of two coupled differential equations; one of them is the hyperbolic equation (an analogue of the Lamé system) and the other is the parabolic equation (an analogue of the diffusion Maxwell system). We prove an existence and uniqueness result. The proof is based on the classical Faedo–Galerkin method.     

Convergence and conditioning of a Nyström method for Stokes flow in exterior three-dimensional domains

Convergence and conditioning results are presented for the lowest-order member of a family of Nyström methods for arbitrary, exterior, three-dimensional Stokes flow. The flow problem is formulated in terms of a recently introduced two-parameter, weakly singular boundary integral equation of the second kind. In contrast to methods based on product integration, coordinate transformation and singularity subtraction, the family of Nyström methods considered here is based on a local polynomial correction determined by an auxiliary system of moment equations. The polynomial correction is designed to remove the weak singularity in the integral equation and provide control over the approximation error. Here we focus attention on the lowest-order method of the family, whose implementation is especially simple. We outline a convergence theorem for this method and illustrate it with various numerical examples. Our examples show that well-conditioned, accurate approximations can be obtained with reasonable meshes for a range of different geometries.     

基于神经网络的一类非仿射非线性系统自适应控制

针对一类带有不确定性的非仿射非线性系统,利用Backstepping设计方法,设计了一种神经网络自适应控制器．该控制器可以实现跟踪特性．基于Lyapunov函数,得出稳定的权学习算法．并利用Lyapunov稳定性理论证明了闭环系统是一致最终有界的．仿真结果表明,这种控制器具有良好的鲁棒性和跟踪特性．     

A New Jacobi Elliptic Function Expansion Method for Solvinga Nonlinear PDE Describing Pulse Narrowing Nonlinear Transmission Lines

In this article, we apply the first elliptic function equation to find a new kind of solutions of nonlinear partial differential equations (PDEs) based on the homogeneous balance method, the Jacobi elliptic expansion method and the auxiliary equation method. New exact solutions to the Jacobi elliptic functions of a nonlinear PDE describing pulse narrowing nonlinear transmission lines are given with the aid of computer program, e.g. Maple or Mathematica. Based on Kirchhoff's current law and Kirchhoff's voltage law, the given nonlinear PDE has been derived and can be reduced to a nonlinear ordinary differential equation (ODE) using a simple transformation. The given method in this article is straightforward and concise, and can be applied to other nonlinear PDEs in mathematical physics. Further results may be obtained.     

Effect of interconvertibility of electromagnetic and gravitational waves in strong external electromagnetic fields and the propagation of waves in the field of a charged “black hole” : PMM vol. 38, n 6, 1974, pp. 1122–1129

It is shown that the behavior of an arbitrary wave propagating in the field of a nonrotating charged black hole is defined (with the use of quadratures) by four functions. Each of these functions obeys its second order equation of the wave kind. Short electromagnetic waves falling onto a black hole are reflected by its field in the form of gravitational and electromagnetic waves whose amplitude was explicitly determined. In the case of the wave carrying rays winding around the limit cycle the reflection and transmission coefficients were obtained in the form of analytic expressions.Various physical processes taking place inside, as well as outside a collapsing star, may induce perturbations of the gravitational, electromagnetic and other fields, and lead to the appearance in the surrounding space of waves of various kinds which propagate over a distorted background and are dissipated along its inhomogeneities.In the absence of rotation and charge in a star, the analysis of small perturbations of the gravitational fields is based on the system of Einstein equations linearized around the Schwarzschild solution. In [1, 2] this system of equations, after expansion of perturbations in spherical harmonics and Fourier transformation with respect to time, was reduced to two independent linear ordinary differential equations of second order of the form of the stationary Schrödinger equation for a particle in a potential force field. Each of these equations defines one of two possible independent perturbation kinds: “even” and “odd” (the different behavior of spherical tensor harmonics at coordinate inversion is the deciding factor in the determination of the kind of perturbation [1, 2]). Although these equations were derived with the superposition on the perturbations of the metric of specific coordinate conditions, they define, as shown in [4], the behavior of invariants of the perturbed gravitational field, which imparts to the potential barriers appearing in these equations an invariant meaning.The system of Maxwell equations on the background of Schwarzschild solution also reduces to similar equations, which differ from the above only by the form of potential barriers appearing in these [5].In the presence in the unperturbed solution of a strong electromagnetic field the gravitational and electromagnetic waves interact with each other, and transmutation takes place. The train of short periodic electromagnetic waves generates the accompanying train of gravitational waves. This phenomenon was first analyzed in [6] on and arbitrary background. It was shown in [7, 8] that dense stars surrounded by hot plasma may acquire a charge owing to splitting of charges by radiation pressure and the “sweeping out” of positrons nascent in vapors in strong electrostatic fields. The interaction of waves becomes particularly clearly evident in the neighborhood of black holes which may serve as “valves” by maintaining equilibrium between the relict electromagnetic and gravitational radiation in the Universe. Rotation of black holes intensifies this effect [6].If a nonrotating star possesses an electrostatic charge, the definition of perturbations of the electromagnetic and gravitational fields must be based on the complete system of Einstein-Maxwell equations linearized around the Nordström-Reissner solution. (Small perturbations of electromagnetic field outside a charged black hole were considered in [9, 10] on the basis of the system of Maxwell equations on a “rigid” background of the Nordström-Reissner solution, without taking into account the interconvertibility of gravitational and electromagnetic waves, which materially affects their behavior in the neighborhood of a charged black hole). Here this system of equations which define the interacting gravitational and electromagnetic perturbations are reduced to four independent second order differential equations, two for each kind of perturbations (an importsnt part is played here by the coordinate conditions imposed on the perturbations of the metric, proposed by the authors in [4]). Perturbation components of the metric and of the electromagnetic field are determined in quadratures by the solutions of these equations. If the charge of a star tends to vanish, two of the derived equations convert to equations for gravitational waves on the background of the Schwarzschild solution [1, 2], while the twoothers become equations which are equivalent to Maxwell solutions on the same background. The short-wave asymptotics of derived equations is determined throughout including the neighborhood of the limit cycle for the wave carrying rays. These solutions far away from the point of turn coincide with those obtained in [6] for any arbitrary background. Approximation of geometric optics does not provide correct asymptotics for impact parameters of rays which are close to critical for which the Isotropie and geodesic parameters wind around the limit cycle. This case is investigated below.A similar situation in the Schwarzschild field was analyzed in [11], where analytic expressions for the wave reflection and transmission coefficients were determined, and the integral radiation stream trapped by a black hole produced by another radiation component of the dual system was calculated.     

The method of normal splines for linear DAEs on the number semi-axis

The method of normal spline-collocation (NSC), applicable to a wide class of ordinary linear singular differential and integral equations, is specified for the boundary value problems for differential-algebraic equations of second order on the number semi-axis. The method consists in minimization of a norm of the collocation systems' solutions in an appropriate Hilbert–Sobolev space. The NSC method does not use the notion of differentiation index and it is applicable to DAEs of any index as well as to equations not reducible to the normal form. The problems on the infinite interval can be solved in two ways. The first way is based on the use of the original space of functions defined on the semi-axis, and the second way is based on a singular transformation of the semi-axis into the unit segment. A new reproducing kernel, that provides the first way, is presented. An algorithm to create a non-uniform collocation grid is described.     

Solitons in coupled nonlinear Schrödinger equations with variable coefficients

We study the coupled nonlinear Schrodinger equation with variable coefficients (VCNLS), which can be used to describe the interaction among the modes in nonlinear optics and Bose–Einstein condensation. By constructing an explicit transformation, which maps VCNLS to the classical coupled nonlinear Schrödinger equations (CNLS), we obtain Bright–Dark and Bright–Bright solitons for VCNLS. Furthermore, the optical super-lattice potentials (or periodic potentials) and hyperbolic cosine potentials with parameters are designed, which are two kinds of important potentials in physics. This method can be used to design a large variety of external potentials in VCNLS, which could be meaningful for manipulating solitons experimentally.     

On the Comparisons of Unit Dual Quaternion and Homogeneous Transformation Matrix

This paper reveals the differences and similarities between two popular unified representations, i.e. the UDQ (unit dual quaternion) and the HTM (homogeneous transformation matrix), for transformation in the solution to the kinematic problem, in order to provide a clear, concise and self-contained introduction into dual quaternions and to further present a cohesive view for the UDQ and HTM representations as used in robotics. Specifically, after investigating some fundamental algebraic properties of the UDQ, it is revealed that the kinematical equations represented by the UDQ and the HTM are accordant, and afterwards the direct relationship of UDQ-based error kinematical models in spatialframe and in body-frame are further discussed, with conclusion that either error kinematic model can be chosen for designing kinematical control laws. Finally, the comparative study on the proportional control algorithms based on the logarithmical mapping of the HTM and the UDQ shows that the UDQ-based control law is indeed higher in computational efficiency.     

Stability analysis of FDTD to UPML for time dependent Maxwell equations

We study an finite-difference time-domain (FDTD) system of uniaxial perfectly matched layer (UPML) method for electromagnetic scattering problems. Particularly we analyze the discrete initial-boundary value problems of the transverse magnetic mode (TM) to Maxwell’s equations with Yee’s algorithm. An exterior domain in two spacial dimension is truncated by a square with a perfectly matched layer filled by a certain artificial medium. Besides, an artificial boundary condition is imposed on the outer boundary of the UPML. Using energy method, we obtain the stability of this FDTD system on the truncated domain. Numerical experiments are designed to approve the theoretical analysis.