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.

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.

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.

Analysis of Maxwell’s Equations. Cosmic Magnetism

According to Hypersphere World-Universe Model, dark matter particles DIRACs are magnetic dipoles consisting of two Dirac’s monopoles. We conclude that DIRACs are the subject of Maxwell’s equations. So-called “auxiliary” magnetic field intensity H is indeed current density of magnetic dipoles. The developed approach to magnetic field can explain a wealth of discovered phenomena in Cosmic Magnetism: a dark magnetic field, the large-scale structure of the Milky Way’s magnetic field, and other magnetic phenomena which are only partly related to objects visible in other spectral ranges.

A New Formulation of Maxwell’s Equations in Clifford Algebra

A new unification of the Maxwell equations is given in the domain of Clifford algebras with in a fashion similar to those obtained with Pauli and Dirac algebras. It is shown that the new electromagnetic field multivector can be obtained from a potential function that is closely related to the scalar and the vector potentials of classical electromagnetics. Additionally it is shown that the gauge transformations of the new multivector and its potential function and the Lagrangian density of the electromagnetic field are in agreement with the transformation rules of the second-rank antisymmetric electromagnetic field tensor, in contrast to the results obtained by applying other versions of Clifford algebras.

Re-Thinking Biology—I. Maxwell’s Demon and the Spontaneous Origin of Life

A hypothesis is advanced in which life began from a Darwinian selection among a diversification of molecular species containing the phosphate moiety which broke the constraints implicit in the Second Law of Thermodynamics, discussed famously by Schrodinger, by obtaining energy from specific infrared frequencies located in the phosphorus-oxygen vibration at a frequency around 1000 cm-1. We propose the source of this energy was from the internal conversion of solar broadband energy by the phosphate mineral Apatite, present at the bottom of a primitive biogenesis pond. In this scenario, life is re-defined as being hotter than its environment and as using its excess energy, supplied by infra-red conversion, to react with its molecular environment and pump itself up the “entropy slope” thereby;replication is through breakages of increasingly large phosphate containing biopolymers. The idea has implications for modern explanations of living systems.

Solution of Maxwell’s Equations for Cylindrical Symmetry Waveguides

The solution of Maxwell’s equations for a piecewise homogeneous medium of cylindrical symmetry has been obtained. The parameters of the cylindrical waveguide modes have been calculated on its basis. The conclusions are confirmed by numerical calculation of the first four modes of a hollow metal waveguide operating as a mode convector.

Maxwell’s Equations as the Basis for Model of Atoms

A century ago the classical physics couldn’t explain many atomic physical phenomena. Now the situation has changed. It’s because within the framework of classical physics with the help of Maxwell’s equations we can derive Schrödinger’s equation, which is the foundation of quantum physics. The equations for energy, momentum, frequency and wavelength of the electromagnetic wave in the atom are derived using the model of atom by analogy with the transmission line. The action constant A0 = (μ0/ε0)1/2s02e2 is a key term in the above mentioned equations. Besides the other well-known constants, the only unknown constant in the last expression is a structural constant of the atom s0. We have found that the value of this constant is 8.277 56 and that it shows up as a link between macroscopic and atomic world. After calculating this constant we get the theory of atoms based on Maxwell’s and Lorentz equations only. This theory does not require knowledge of Planck’s constant h, which is replaced with theoretically derived action constant A0, while the replacement for the fine structure constant α-1 is theoretically derived expression 2s02 = 137.036. So, the structural constant s0 replaces both constants h and α. This paper also defines the stationary states of atoms and shows that the maximal atomic number is equal to Zmax = 137. The presented model of the atoms covers three of the four fundamental interactions, namely the electromagnetic, weak and strong interactions.

Derivation of Maxwell’s Equations via the Covariance Requirements of the Special Theory of Relativity, Starting with Newton’s Laws

The purpose of this paper is to establish a connection between Maxwell’s equations, Newton’s laws, and the special theory of relativity. This is done with a derivation that begins with Newton’s verbal enunciation of his first two laws. Derived equations are required to be covariant, and a simplicity criterion requires that the four-vector force on a charged particle be linearly related to the four-vector velocity. The connecting tensor has derivable symmetry properties and contains the electric and magnetic field vectors. The Lorentz force law emerges, and Maxwell’s equations for free space emerge with the assumption that the tensor and its dual must both satisfy first-order partial differential equations. The inhomogeneous extension yields a charge density and a current density as being the source of the field, and yields the law of conservation of charge. Newton’s third law is reinterpreted as a reciprocity statement, which requires that the charge in the source term can be taken as the same physical entity as that of the test particle and that both can be assigned the same units. Requiring covariance under either spatial inversions or time reversals precludes magnetic charge being a source of electromagnetic fields that exert forces on electric charges.

Determination of Photon and Elementary Particles Rest Masses Using Maxwell’s Equations and Generalized Potential Dependent Special Relativity

The nature and origin of the photon and elementary rest masses are some of the challeng-ing problems that physics face. The approaches used to solve these problems are complex and time-consuming. Specifically, the photon rest mass pays attention to theoretical physi-cists. Many experimental works show that the photon rest mass is non zero. This problem can be solved using generalized potential dependent special relativity, which has been de-rived using simple arguments, and Maxwell’s equations, besides the conventional Einstein energy-momentum relation. The results obtained show that the rest mass of photons and elementary particles are strongly dependent on the vacuum energy and a universal con-stant. This result conforms with the models that predict time decaying vacuum energy as-sociated with production of smaller rest mass particles followed by larger masses. The two potential dependent mass expressions conform with the cosmological models that suggest the photon is generated first by assuming the universe consisting of total constant vacuum with decaying cosmological part and mass generating part. Using Maxwell’s equations, beside plank and De Broglie hypothesis together with special relativity energy-momentum relation the photon rest mass is estimated. It was shown that the photon rest mass is ex-tremely small compared to the electron mass.

The Dynamic Gravitation of Photons from the Perspective of Maxwell’s Wave Equations

Although the gravitational constant (G) does not explicitly occur in the Maxwell Wave Equations, this paper will show that G is indeed implicitly contained in them. The logical consequence hereby is that electromagnetic radiation is associated with dynamic gravitation and not—as assumed in Einstein’s Special Theory of Relativity—with “static” gravitation, dynamic gravitation being at the time unknown. According to the Maxwell Wave Equations, gravitation experiences the same dynamic (speed of light c) as electromagnetic radiation and must therefore also be of a quantum nature. There must exist an equal number of gravitational quanta as there are photons. Since photons do not possess a baryonic rest mass but only a relativistic mass, this mass must be nonbaryonic in nature—precisely as their dynamic gravitation.

A Fresh View for Maxwell’s Equations and Electromagnetic Wave Propagation

Equations related with wave propagation are reexamined as in certain circumstances law of conservation of energy is not fulfilled even though it is cautiously clarified with the help of Heisenberg’s uncertainty principle. Recently, attempt has also been made to understand certain discrepancies in optical phenomena like diffraction or interference. The purpose of the present investigation, therefore, is to overcome some discrepancies by introducing constants of integration in Maxwell’s Equation. It turns out that the presence of vibrating strings (or store energy) in the medium becomes essential to understand several details of the wave propagation.

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.

Recent Progress on the Maxwell’s Equations for Describing a Mechano-Driven Medium System with Multiple Moving Objects/Media

Maxwell’s equations for a mechano-driven media system(MEs-f-MDMS)have been used to characterize the electromagnetism of multislow-moving media that may be accelerated with complex trajectories.Such an approach starts from the integral forms of the four physics laws and is different from the classical approach of using the Lorentz transformation for correlating the electromagnetic phenomena observed in two inertial reference frames with relative motion.The governing equations inside the moving object/medium are the MEs-f-MDMS,and those in vacuum are the classical Maxwell’s equations;the full solutions of both reconcile at the medium surface/interface and satisfy the boundary conditions.This paper reviews the background,physical principle,and mathematical derivations for formulating the MEs-f-MDMS.Strategies are also presented for mathematically solving the MEs-f-MDMS.The unique advances made by the MEs-f-MDMS have been systematically summarized,as are their potential applications in engineering.We found that the Lorentz transformation is perfect for treating the electromagnetic phenomena of moving point charges in vacuum;however,for moving objects,the covariance of Maxwell’s equations may not hold,and use of the MEs-f-MDMS may be required if the velocity is low.Finally,recent advances for treating the boundary conditions at the nanoscale without assuming an abrupt boundary are also reviewed.

ON DISCRETE ENERGY DISSIPATION OF MAXWELL’S EQUATIONS IN A COLE-COLE DISPERSIVE MEDIUM

A simple criterion is studied for the first time for identifying the discrete energy dissipation of the Crank-Nicolson scheme for Maxwell’s equations in a Cole-Cole dispersive medium.Several numerical formulas that approximate the time fractional derivatives are investigated based on this criterion,including the L1 formula,the fractional BDF-2,and the shifted fractional trapezoidal rule(SFTR).Detailed error analysis is provided within the framework of time domain mixed finite element methods for smooth solutions.The convergence results and discrete energy dissipation law are confirmed by numerical tests.For nonsmooth solutions,the method SFTR can still maintain the optimal convergence order at a positive time on uniform meshes.Authors believe this is the first appearance that a second-order time-stepping method can restore the optimal convergence rate for Maxwell’s equations in a Cole-Cole dispersive medium regardless of the initial singularity of the solution.

Implicit DG Method for Time Domain Maxwell’s Equations Involving Metamaterials

An implicit discontinuous Galerkin method is introduced to solve the timedomain Maxwell’s equations in metamaterials.The Maxwell’s equations in metamaterials are represented by integral-differential equations.Our scheme is based on discontinuous Galerkin method in spatial domain and Crank-Nicolson method in temporal domain.The fully discrete numerical scheme is proved to be unconditionally stable.When polynomial of degree at most p is used for spatial approximation,our scheme is verified to converge at a rate of O(τ^(2)+h^(p)+1/2).Numerical results in both 2D and 3D are provided to validate our theoretical prediction.

Adaptive hp-FEM with Arbitrary-Level Hanging Nodes for Maxwell’s Equations

Adaptive higher-order finite element methods(hp-FEM)are well known for their potential of exceptionally fast(exponential)convergence.However,most hp-FEM codes remain in an academic setting due to an extreme algorithmic complexity of hp-adaptivity algorithms.This paper aims at simplifying hpadaptivity for H(curl)-conforming approximations by presenting a novel technique of arbitrary-level hanging nodes.The technique is described and it is demonstrated numerically that it makes adaptive hp-FEM more efficient compared to hp-FEM on regular meshes and meshes with one-level hanging nodes.

A High-OrderDiscontinuous GalerkinMethod for the Two-Dimensional Time-Domain Maxwell’s Equations on Curved Mesh

In this paper,a DG(Discontinuous Galerkin)method which has been widely employed in CFD(Computational Fluid Dynamics)is used to solve the twodimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh.The element interfaces on solid boundary are treated in both curved way and straight way.Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids,where the high-order convergence in accuracy can be observed.Both the curved and the straight solid boundary implementation can give accurate RCS(Radar Cross-Section)results with sufficiently small mesh size,but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size.More importantly,this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.

Analysis and Application of Stochastic Collocation Methods for Maxwell’s Equations with Random Inputs

In this paper we develop and analyze the stochastic collocation method for solving the time-dependent Maxwell’s equations with random coefficients and subject to random initial conditions.We provide a rigorous regularity analysis of the solution with respect to the random variables.To our best knowledge,this is the first theoretical results derived for the standard Maxwell’s equations with random inputs.The rate of convergence is proved depending on the regularity of the solution.Numerical results are presented to confirm the theoretical analysis.

Realization of Quantum Maxwell's Demon with Solid-State Spins

We report experimental realization of a quantum version of Maxwell＇s demon using solid state spins where the information acquiring and feedback operations by the demon are achieved through conditional quantum gates.A unique feature of this implementation is that the demon can start in a quantum superposition state or in an entangled state with an ancilla observer. Through quantum state tomography, we measure the entropy in the system, demon, and the ancilla, showing the influence of coherence and entanglement on the result. A quantum implementation of Maxwell＇s demon adds more controllability to this paradoxical thermal machine and may find applications in quantum thermodynamics involving microscopic systems.