Maxwell Equations and Magnetic Monopoles

The manuscript introduces an “ab initio” quantum model to deduce the Maxwell equations. After general considerations and laying out the model’s theoretical framework, these equations can be derived alongside a broad variety of other results. Specifically, a corollary of the present model proposes a possible mechanism underlying the formation of magnetic monopoles and allows estimating their formation energy in order of magnitude.

STABILITY AND SUPER CONVERGENCE ANALYSIS OF ADI-FDTD FOR THE 2D MAXWELL EQUATIONS IN A LOSSY MEDIUM

Several new energy identities of the two dimenslonal（2D） Maxwell equations in a lossy medium in the case of the perfectly electric conducting boundary conditions are proposed and proved. These identities show a new kind of energy conservation in the Maxwell system and provide a new energy method to analyze the alternating direction im- plicit finite difference time domain method for the 2D Maxwell equations （2D-ADI-FDTD）. It is proved that 2D-ADI-FDTD is approximately energy conserved, unconditionally sta- ble and second order convergent in the discrete L2 and H1 norms, which implies that 2D-ADI-FDTD is super convergent. By this super convergence, it is simply proved that the error of the divergence of the solution of 2D-ADI-FDTD is second order accurate. It is also proved that the difference scheme of 2D-ADI-FDTD with respect to time t is second order convergent in the discrete H1 norm. Experimental results to confirm the theoretical analysis on stability, convergence and energy conservation are presented.

Noether conserved quantities and Lie point symmetries of difference Lagrange-Maxwell equations and lattices

This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.

SPACE-TIME DISCONTINUOUS GALERKIN METHOD FOR MAXWELL EQUATIONS IN DISPERSIVE MEDIA

In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability and error estimate of order O？τr＋1＋hk＋1/2？ are obtained when polynomials of degree at most r and k are used for the temporal dis-cretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r＋1 in temporal variable t.

A STOCHASTIC GALERKIN METHOD FOR MAXWELL EQUATIONS WITH UNCERTAINTY

In this article,we investigate a stochastic Galerkin method for the Maxwell equations with random inputs.The generalized Polynomial Chaos(gPC)expansion technique is used to obtain a deterministic system of the gPC expansion coefficients.The regularity of the solution with respect to the random is analyzed.On the basis of the regularity results,the optimal convergence rate of the stochastic Galerkin approach for Maxwell equations with random inputs is proved.Numerical examples are presented to support the theoretical analysis.

A wavelet multiscale method for inversion of Maxwell equations

This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is introduced and applied to the inversion of Maxwell equations. The inverse problem is decomposed into multiple scales with wavelet transform, and hence the original problem is reformulated to a set of sub-inverse problems corresponding to different scales, which can be solved successively according to the size of scale from the shortest to the longest. The stable and fast regularized Gauss-Newton method is applied to each scale. Numerical results show that the proposed method is effective, especially in terms of wide convergence, computational efficiency and precision.

ANALYSIS OF FDTD TO UPML FOR MAXWELL EQUATIONS IN POLAR COORDINATES

An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

Complex Maxwell’s Equations

Maxwell’s equations in electromagnetism can be categorized into three dis-tinct groups based on the electromagnetic source when employing quaterni-ons. Each group represents a self-contained system in which Maxwell’s equations are applied and validated concurrently, in contrast to the previous approach that did not account for this. It has been noted that the formulation of these Maxwell equations ultimately results in the formulation of Max-well’s equations utilizing the scalar function.

GLOBAL SUPERCONVERGENCE OF THE MIXED FINITE ELEMENT METHODS FOR 2-D MAXWELL EQUATIONS

Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.

Effective Maxwell Equations from Time-dependent Density Functional Theory

The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.

An Accurate Numerical Scheme for Maxwell Equation with CIP-Method of Characteristics

A new multi-dimensional scheme for the Maxwell equations is established by the CIP method in combination with the method of characteristics(CIP-MOC).In addition,the CIP-MOC can be extended to arbitrary grid system by the Soroban grid without losing the third-order accuracy.With the accuracy fixed,the grid points required for the CIP are 40 times less than the conventional schemes like the FDTD in three dimensions.Numerical solutions obtained by the CIP-MOC are compared with analytical solution and the FDTD in plane-wave scattering by a perfectly-conducting circular cylinder,and the CIP-MOC agrees very well with analytical solutions.The Soroban grid is also applied to the Vlasov equation that describes the kinematics of plasmas that is frequently combined with the Maxwell equation.The adaptively moving points in velocity space are similar to the particle codes but can provide accurate solutions.

Solution of time-domain Maxwell equation with PML by using modified Laguerre polynomials

In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of the modified Laguerre polynomials, and then the original equations of the initial value and boundary value can be transformed into a series of problems independent of the time variable. In this case the method of finite difference （FD）, the finite element method （FEM）, the method of moment （MoM）, etc. or the combination of these methods can be used to solve the problems. Finally, a numerical model is provided for the scattering problem with perfect matched layer （PML） by using FD. The comparison between the results of the proposed method and FDTD is presented to verify the proposed new method.

Dirac and Maxwell Systems in Split Octonions

The known equivalence of 8-dimensional chiral spinors and vectors, also referred to as triality, is discussed for (4 + 4)-space. Split octonionic representation of SO(4, 4) and Spin(4, 4) groups and the trilinear invariant form are explicitly written and compared with Clifford algebraic matrix representation. It is noted that the complete algebra of split octonionic basis units can be recovered from the Moufang and Malcev relations for the three vector-like elements. Lagrangians on split octonionic fields that generalize Dirac and Maxwell systems are constructed using group invariant forms. It is shown that corresponding equations are related to split octonionic analyticity conditions.

Adaptive Maxwell’s Equations Derived Optimization and Its Application in Antenna Array Synthesis

In this paper,a self-adaptive method for the Maxwell’s Equations Derived Optimization(MEDO)is proposed.It is implemented by applying the Sequential Model-Based Optimization(SMBO)algorithm to the iterations of the MEDO,and achieves the automatic adjustment of the parameters.The proposed method is named as adaptive Maxwell’s equations derived optimization(AMEDO).In order to evaluate the performance of AMEDO,eight benchmarks are used and the results are compared with the original MEDO method.The results show that AMEDO can greatly reduce the workload of manual adjustment of parameters,and at the same time can keep the accuracy and stability.Moreover,the convergence of the optimization can be accelerated due to the dynamical adjustment of the parameters.In the end,the proposed AMEDO is applied to the side lobe level suppression and array failure correction of a linear antenna array,and shows great potential in antenna array synthesis.

Nonlinear waves of the Hirota and the Maxwell Bloch equations in nonlinear optics

In this paper, considering the Hirota and the Maxwell–Bloch （H-MB） equations which are governed by femtosecond pulse propagation through a two-level doped fiber system, we construct the Darboux transformation of this system through a linear eigenvalue problem. Using this Daurboux transformation, we generate multi-soliton, positon, and breather solutions （both bright and dark breathers） of the H-MB equations. Finally, we also construct the rogue wave solutions of the above system.

HIGH-ORDER EXPLICIT TIME-INTEGRATORS FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE MAXWELL EQUATIONS

We investigate the practical implementation of a high-order explicit time-stepping method based on polynomial approximations,for possible application to large-scale problems in electromagnetics.After the spatial discretization by a high-order discontinuous Galerkin method,we obtain a linear system of differential equations of the form,∂tY(t)=HY(t)+S(t),where H is a matrix containing the spatial derivatives and t is the time variable.The formal solution can be written in terms of the matrix exponential,exp(tH),acting on some vectors.We introduce a general family of time-integrators based on the approximation of exp(tH)by Jacobi polynomial expansions.We discuss the efficient implementation of this technique,and based on some test problems,we compare the virtues and shortcomings of the algorithm.We also demonstrate how these schemes provide an efficient alternative to standard explicit integrators for computing solutions over long time intervals.

Splitting Finite Difference Methods on Staggered Grids for the Three-Dimensional Time-Dependent Maxwell Equations

In this paper,we study splitting numerical methods for the three-dimensional Maxwell equations in the time domain.We propose a new kind of splitting finitedifference time-domain schemes on a staggered grid,which consists of only two stages for each time step.It is proved by the energy method that the splitting scheme is unconditionally stable and convergent for problems with perfectly conducting boundary conditions.Both numerical dispersion analysis and numerical experiments are also presented to illustrate the efficiency of the proposed schemes.

High-Order Low Dissipation Conforming Finite-Element Discretization of the Maxwell Equations

In this paper,we study high order discretization methods for solving the Maxwell equations on hybrid triangle-quad meshes.We have developed high order finite edge element methods coupled with different high order time schemes and we compare results and efficiency for several schemes.We introduce in particular a class of simple high order low dissipation time schemes based on a modified Taylor expansion.