Sufficient uniqueness conditions for the solution of the time harmonic Maxwell’s equations associated with surface impedance boundary conditions

We consider the time-harmonic Maxwell’s equations for the scattering or radiating problem from a 3-D object that is either a metallic surface coated with material layers (MCS) or a dichroic structure (DS) made up of multiple frequency selective surfaces (FSS) embedded in materials. Low or high order impedance boundary conditions (IBC) are employed to reduce the numerical complexity of the solution of this problem via an integral equation or a finite element formulation. An IBC links the tangential components of the electric field to those of the magnetic field on the outer surface of the MCS, or on the FSSs, and avoids the solution of Maxwell’s equations inside the inhomogeneous domain for a MCS or, for a DS, the meshing of the numerous unit cells in a FSS. Sufficient uniqueness conditions (SUC) are established for the solutions of Maxwell’s equations associated with these IBCs, the performances of which, when constrained by the corresponding SUCs, are numerically evaluated for an infinite or finite planar structure.     

A linear nonconforming finite element method for Maxwell’s equations in two dimensions. Part I: Frequency domain

We suggest a linear nonconforming triangular element for Maxwell’s equations and test it in the context of the vector Helmholtz equation. The element uses discontinuous normal fields and tangential fields with continuity at the midpoint of the element sides, an approximation related to the Crouzeix–Raviart element for Stokes. The element is stabilized using the jump of the tangential fields, giving us a free parameter to decide. We give dispersion relations for different stability parameters and give some numerical examples, where the results converge quadratically with the mesh size for problems with smooth boundaries. The proposed element is free from spurious solutions and, for cavity eigenvalue problems, the eigenfrequencies that correspond to well-resolved eigenmodes are reproduced with the correct multiplicity.     

Complex symmetries of electrodynamics

We prove that a set of nonsingular free solutions of Maxwell's equations forms a representation of the group obtained by analytic continuation of the Poincaré group to complex values of the group parameters, and that a set of singular solutions forms a representation of the group obtained by analytic continuation of the conformal group to complex values of the group parameters. These results are obtained by constructing a theory governing 2 × 2 complex matrix fields defined for complex values of position and time; the equations of this theory are invarient with respect to complex Poincaré transformations and complex conformal transformations, but the set of nonsingular solutions is in one-to-one correspondence with a set of nonsingular solutions of Maxwell's equations, and a similar correspondence exists for the singular solutions. Certain collections of solutions of Maxwell's equations for the field of a current form representations of these complex groups if both magnetic and electric currents are permitted, in which case complex transformations provide a natural connection between electric and magnetic charge. A class of complex transformations also yield natural relations between sources moving slower than light and sources moving faster than light.     

Superconvergence of mixed finite element approximations to 3-D Maxwell’s equations in metamaterials

Numerical simulation of metamaterials has attracted more and more attention since 2000, after the first metamaterial with negative refraction index was successfully constructed. In this paper we construct a fully-discrete leap-frog type finite element scheme to solve the three-dimensional time-dependent Maxwell’s equations when metamaterials are involved. First, we obtain some superclose results between the interpolations of the analytical solutions and finite element solutions obtained using arbitrary orders of Raviart–Thomas–Nédélec mixed spaces on regular cubic meshes. Then we prove the superconvergence result in the discrete l₂ norm achieved for the lowest-order Raviart–Thomas–Nédélec space. To our best knowledge, such superconvergence results have never been obtained elsewhere. Finally, we implement the leap-frog scheme and present numerical results justifying our theoretical analysis.     

Generalization of the FDTD algorithm for simulations of hydrodynamic nonlinear Drude model

In this paper we present a numerical method for solving a three-dimensional cold-plasma system that describes electron gas dynamics driven by an external electromagnetic wave excitation. The nonlinear Drude dispersion model is derived from the cold-plasma fluid equations and is coupled to the Maxwell’s field equations. The Finite-Difference Time-Domain (FDTD) method is applied for solving the Maxwell’s equations in conjunction with the time-split semi-implicit numerical method for the nonlinear dispersion and a physics based treatment of the discontinuity of the electric field component normal to the dielectric-metal interface. The application of the proposed algorithm is illustrated by modeling light pulse propagation and second-harmonic generation (SHG) in metallic metamaterials (MMs), showing good agreement between computed and published experimental results.     

Computing for second harmonic generation in one-dimensional nonlinear photonic crystals

A numerical study of second harmonic generation (SHG) in one-dimensional nonlinear photonic crystals based on full nonlinear system of equations, implemented by a combination of the method of finite elements and fixed-point iterations, is reported. This model is derived from a nonlinear system of Maxwell’s equations, which partly overcomes the known shortcoming of some existing models relied on the undepleted-pump approximation. We derive a general solution of SHG in one-dimensional nonlinear photonic crystals structures. The convergence of our method is fast. Numerical simulations also show the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental wave are located at the photonic band edges or are assigned to the designed defect states.     

Photon spin and Bose-Einstein statistics

It is shown that Maxwell’s equations for the electric and magnetic fields free of sources can be inferred from Dirac’s pair of first-order equations for a zero-mass, zero-charge particle. This result is interpreted as a Lorentz invariant form of the transverse nature of photonic propagation in which only two components of the spin-1 field exist in nature.Canonical quantization of Dirac’s equations leads to a time average of the electromagnetic energy in agreement with the standard result of quantum electrodynamics. It is shown that the spin-statistics theorem is not violated for canonical quantization of the Dirac field provided the mass of the particle is zero.     

Simulation of laser–plasma interactions and fast-electron transport in inhomogeneous plasma

A new framework is introduced for kinetic simulation of laser–plasma interactions in an inhomogeneous plasma motivated by the goal of performing integrated kinetic simulations of fast-ignition laser fusion. The algorithm addresses the propagation and absorption of an intense electromagnetic wave in an ionized plasma leading to the generation and transport of an energetic electron component. The energetic electrons propagate farther into the plasma to much higher densities where Coulomb collisions become important. The high-density plasma supports an energetic electron current, return currents, self-consistent electric fields associated with maintaining quasi-neutrality, and self-consistent magnetic fields due to the currents. Collisions of the electrons and ions are calculated accurately to track the energetic electrons and model their interactions with the background plasma. Up to a density well above critical density, where the laser electromagnetic field is evanescent, Maxwell’s equations are solved with a conventional particle-based, finite-difference scheme. In the higher-density plasma, Maxwell’s equations are solved using an Ohm’s law neglecting the inertia of the background electrons with the option of omitting the displacement current in Ampere’s law. Particle equations of motion with binary collisions are solved for all electrons and ions throughout the system using weighted particles to resolve the density gradient efficiently. The algorithm is analyzed and demonstrated in simulation examples. The simulation scheme introduced here achieves significantly improved efficiencies.     

The research on polarization states with local pumping in two-dimensional active random medium

The polarization-dependent competition on inversed population with local pumping in two-dimensional (2D) active random media is investigated by simultaneously solving Maxwell’s equations and rate equations of electronic population for both transverse magnetic (TM) and transverse electric (TE) modes. The threshold property of the two polarization states with local pumping is analyzed. Results show that the lasing threshold of TM modes is lower than that of TE ones when the radius of the local pumping region is suitably selected, and the case would be contrary when the selected radius of the local pumping region is excessive.     

Scattering by random rough surfaces: Study of direct and inverse problem

In order to study the problems of scattering by rough metallic surfaces, we have used Maxwell’s equations in covariant form within the framework of a non-orthogonal coordinates system adapted to the geometry of the problem. Electromagnetic fields are written in Fourier’s integral form. The solution is found by using a perturbation method applied to the smooth surface problem; this is fully justified when the defects are of small magnitude.For the direct problem, the mean value of diffraction intensity is obtained for random rough surfaces of finite conductivity by computer simulation.In the case of the inverse problem, the reconstruction of the profile of the metal surface from values of the diffraction intensity, obtained by simulation, is found using an iterative algorithm.     

Tensor Lagrangians,Lagrangians Equivalent to the Hamilton-Jacobi Equation and Relativistic Dynamics

We deal with Lagrangians which are not the standard scalar ones. We present a short review of tensor Lagrangians, which generate massless free fields and the Dirac field, as well as vector and pseudovector Lagrangians for the electric and magnetic fields of Maxwell’s equations with sources. We introduce and analyse Lagrangians which are equivalent to the Hamilton-Jacobi equation and recast them to relativistic equations.     

A high accuracy algorithm for 3D periodic electromagnetic scattering

We develop a highly accurate numerical method for scattering of 3D electromagnetic waves by doubly periodic structures. We approximate scattered fields using the Müller boundary integral formulation of Maxwell’s equations. The accuracy is achieved as singularities are isolated through the use of partitions of unity, leaving smooth, periodic integrands that can be evaluated with high accuracy using trapezoid sums. The removed singularities are resolved through a transformation to polar coordinates. The method relies on the ideas used in the free space scattering algorithm of Bruno and Kunyansky.     

Numerical method for hydrodynamic modulation equations describing Bloch oscillations in semiconductor superlattices

We present a finite difference method to solve a new type of nonlocal hydrodynamic equations that arise in the theory of spatially inhomogeneous Bloch oscillations in semiconductor superlattices. The hydrodynamic equations describe the evolution of the electron density, electric field and the complex amplitude of the Bloch oscillations for the electron current density and the mean energy density. These equations contain averages over the Bloch phase which are integrals of the unknown electric field and are derived by singular perturbation methods. Among the solutions of the hydrodynamic equations, at a 70 K lattice temperature, there are spatially inhomogeneous Bloch oscillations coexisting with moving electric field domains and Gunn-type oscillations of the current. At higher temperature (300 K) only Bloch oscillations remain. These novel solutions are found for restitution coefficients in a narrow interval below their critical values and disappear for larger values. We use an efficient numerical method based on an implicit second-order finite difference scheme for both the electric field equation (of drift-diffusion type) and the parabolic equation for the complex amplitude. Double integrals appearing in the nonlocal hydrodynamic equations are calculated by means of expansions in modified Bessel functions. We use numerical simulations to ascertain the convergence of the method. If the complex amplitude equation is solved using a first order scheme for restitution coefficients near their critical values, a spurious convection arises that annihilates the complex amplitude in the part of the superlattice that is closer to the cathode. This numerical artifact disappears if the space step is appropriately reduced or we use the second-order numerical scheme.     

Generalized Split-Octonion Electrodynamics

Starting with the usual definitions of octonions and split octonions in terms of Zorn vector matrix realization, we have made an attempt to write the consistent form of generalized Maxwell’s equations in presence of electric and magnetic charges (dyons). We have thus written the generalized potential, generalized field, and generalized current of dyons in terms of split octonions and accordingly the split octonion forms of generalized Dirac Maxwell’s equations are obtained in compact and consistent manner. This theory reproduces the dynamic of electric (magnetic) in the absence of magnetic (electric) charges.     

电磁波导的辛分析与对偶棱边元

将电磁波导的控制方程导向了Hamilton体系、辛几何的形式.以电磁场的横向分量组成对偶向量并采用分离变量法，可以得到Hamilton算子矩阵的辛本征值问题.共轭辛正交归一关系、辛本征解展开定理等均可在此应用.对于复杂横截面和填充非均匀材料的电磁波导，提出对偶棱边元，对截面半解析离散后即可进行数值求解.对偶棱边元克服了结点基有限元求解电磁场问题的困难，与常规棱边元相比在某些方面具有一定的优势. 关键词： 电磁波导 Hamilton体系 对偶变量 棱边元     

Standstill Electric Charge Generates Magnetostatic Field under Born-Infeld Electrodynamics

The Abelian Born-Infeld classical non-linear Electrodynamic has been used to investigate the electric and magnetostatic fields generated by a point-like electric charge at rest in an inertial frame. The results show a rich internal structure for the charge. Analytical solutions have also been found. Such fields configurations have been interpreted in terms of vacuum polarization and magnetic-like charges produced by the very high strengths of the electric field considered. Apparently non-linearity is responsible for the emergence of an anomalous magnetostatic field suggesting a possible connection to that created by a magnetic dipole composed of two magnetic charges with opposite signs. Consistently in situations where the Born-Infeld field strength parameter is free to become infinite, Maxwell’s regime takes over, the magnetic sector vanishes and the electric field assumes a Coulomb behavior with no trace of a magnetic component. The connection to other monopole solutions, like Dirac’s and ’tHooft’s Poliakov’s types are also discussed. Finally, some speculative remarks are presented in an attempt to explain such fields.