• ### General form of the full electromagnetic Green function in materials physics(1704.07594)

March 11, 2019 physics.class-ph
In this article, we present the general form of the full electromagnetic Green function which is suitable for the application in bulk materials physics. In particular, we show how the seven adjustable parameter functions of the free Green function translate into seven corresponding parameter functions of the full Green function. Furthermore, for both the fundamental response tensor and the electromagnetic Green function, we discuss the reduction of the Dyson equation on the four-dimensional Minkowski space to an equivalent, three-dimensional Cartesian Dyson equation.
• ### An Electron Model with Synchrotron Radiation(1206.0620)

Dec. 4, 2019 physics.class-ph
A classical model of the electron based on Maxwell's equations is presented in which the wave character is described by classical physics. Most properties follow from the description of a classical massless charge circulating with v = c. The magnetic moment of the electron yields the radius of this circulation and the generated synchrotron radiation removes the singularity of the Coulomb field and generates the mass of the electron. Quantum mechanics yields its spin.. The fine structure constant {\alpha} compares this dynamic structure of the electron with the classical point-like static view. The field lines in 3 dimensions within the free electromagnetic fields of the creation process yield possible tracks of the massless charge. In form of a knotted trefoil the electron is described as a stable spin-1/2-particle. The dimensions of an underlying toroid is determined and the size of the charge results in the quantum mechanic environment. When in the radiation the electron is replaced by its whole family the synchrotron radiation becomes strong structures close to the origin. The Coulomb field is then not completely compensated in matter and gravitation results. With its fixed internal revolution frequency this electron appears to the external world as a standing wave with an amplitude propagating like the de Broglie wave. - Keywords: ~ Electron ~ Classical wave model ~ Spherical wave field ~ Elementary charge ~ alpha ~ Mass ~ Wave character - Gravitation
• ### Heat Energy, Entropy and The Second Law: A New Perspective in Thermodynamics(1201.4284)

We present a new state function, denotes the energy of thermal motion, it shows that the energy classification for the internal energy has been completed. In a new theoretical structure, we present the concept of the heat partial pressure, rewrite the equations of the first law, by which, the energy conversion and energy transport can be distinguished by an explicit way, we redefine the concept of the entropy, and confirm a conversion potential related to the dissipative compensation, we discuss in detail on the three sources of irreversibility, and the conversions between the calorimetric entropy and the configurational entropy, an interesting conclusion shows that the second law itself has contained the mechanism of evolution.
• ### Affine and convex spaces: blending the analytic and geometric viewpoints(1104.0032)

This is a short introduction to affine and convex spaces, written especially for physics students. It summarizes different elementary presentations available in the mathematical literature, and blends analytic- and geometric-flavoured presentations. References are also provided, as well as a brief discussion of Grassmann spaces and an example showing the relevance and usefulness of affine spaces in Newtonian physics.
• ### Effective gravity and effective quantum equations in a system inspired by walking droplets experiments(1706.05640)

Jan. 31, 2019 quant-ph, physics.class-ph
In this paper we suggest a macroscopic toy system in which a potential-like energy is generated by a non-uniform pulsation of the medium (i.e. pulsation of transverse standing oscillations that the elastic medium of the system tends to support at each point). This system is inspired by walking droplets experiments with submerged barriers. We first show that a Poincar\'e-Lorentz covariant formalization of the system causes inconsistency and contradiction. The contradiction is solved by using a general covariant formulation and by assuming a relation between the metric associated with the elastic medium and the pulsation of the medium. (Calculations are performed in a Newtonian-like metric, constant in time). We find ($i$) an effective Schr\"odinger equation with external potential, ($ii$) an effective de Broglie-Bohm guidance formula and ($iii$) an energy of the particle' which has a direct counterpart in general relativity as well as in quantum mechanics. We analyze the wave and the particle' in an effective free fall and with a harmonic potential. This potential-like energy is an effective gravitational potential, rooted in the pulsation of the medium at each point. The latter, also conceivable as a natural clock, makes easy to understand why proper time varies from place to place.
• ### Homogenization of Maxwell's equations in layered system beyond static approximation(1708.01661)

The propagation of electromagnetic waves through disordered layered system is considered in the paradigm of Maxwell's equations homogenization. In spite of the impossibility to describe the system in terms of effective dielectric permittivity and/or magnetic permeability the unified way to describe the propagation and Anderson localization of electromagnetic waves is proposed in terms of the introduced effective wave vector (effective refractive index). It is demonstrated that both real and imaginary parts of the effective wave vector (contrary to effective dielectric permittivity and/or magnetic permeability) are the self-averaging quantities at any frequency. The introduced effective wave vector is analytical function of frequency; corresponding Kramers-Kronig relation generalizes the Jones-Herbert-Thouless formula.
• ### Stochastic Evolution of Augmented Born--Infeld Equations(1705.07645)

This paper compares the results of applying a recently developed method of stochastic uncertainty quantification designed for fluid dynamics to the Born-Infeld model of nonlinear electromagnetism. The similarities in the results are striking. Namely, the introduction of Stratonovich cylindrical noise into each of their Hamiltonian formulations introduces stochastic Lie transport into their dynamics in the same form for both theories. Moreover, the resulting stochastic partial differential equations (SPDE) retain their unperturbed form, except for an additional term representing induced Lie transport by the set of divergence-free vector fields associated with the spatial correlations of the cylindrical noise. The explanation for this remarkable similarity lies in the method of construction of the Hamiltonian for the Stratonovich stochastic contribution to the motion in both cases; which is done via pairing spatial correlation eigenvectors for cylindrical noise with the momentum map for the deterministic motion. This momentum map is responsible for the well-known analogy between hydrodynamics and electromagnetism. The momentum map for the Maxwell and Born-Infeld theories of electromagnetism treated here is the 1-form density known as the Poynting vector. Two Appendices treat the Hamiltonian structures underlying these results.
• ### Faraday's Law and Magnetic Induction: cause and effect(1705.08406)

Jan. 10, 2019 physics.class-ph
Faraday's Law of induction is often stated as "a change in magnetic flux causes an EMF"; or, more cautiously, "a change in magnetic flux is associated with an EMF"; It is as well that the more cautious form exists, because the first "causes" form is incompatible with the usual expression $V = - \partial_t \Phi$. This is not, however, to deny the causality as reasonably inferred from experimental observation - it is the equation for Faraday's Law of induction which does not represent the claimed cause-and-effect relationship. Here I investigate a selection of different approaches, trying to see how an explicitly causal mathematical equation, which attempts to encapsulate the "a change in magnetic flux causes ..." idea, might arise.
• ### Does Bohm's Quantum Force Have a Classical Origin?(1409.8271)

Jan. 4, 2019 quant-ph, physics.class-ph
In the de Broglie - Bohm formulation of quantum mechanics, the electron is stationary in the ground state of the hydrogen atom, because the quantum force exactly cancels the Coulomb attraction of the electron to the proton. In this paper it is shown that classical electrodynamics similarly predicts the Coulomb force can be effectively canceled by part of the magnetic force that occurs between two similar particles each consisting of a point charge moving with circulatory motion at the speed of light. Supposition of such motion is the basis of the {\em Zitterbewegung} interpretation of quantum mechanics. The magnetic force between two luminally-circulating charges for separation large compared to their circulatory motions contains a radial inverse square law part with magnitude equal to the Coulomb force, sinusoidally modulated by the phase difference between the circulatory motions. When the particles have equal mass and their circulatory motions are aligned but out of phase, part of the magnetic force is equal but opposite the Coulomb force. This raises a possibility that the quantum force of Bohmian mechanics may be attributable to the magnetic force of classical electrodynamics. It is further shown that non-relativistic relative motion between the particles leads to modulation of the magnetic force with spatial period equal to the de Broglie wavelength.
• ### Linear stability of periodic three-body orbits with zero angular momentum and topological dependence of Kepler's third law: a numerical test(1705.03728)

We test numerically the recently proposed linear relationship between the scale-invariant period $T_{\rm s.i.} = T |E|^{3/2}$, and the topology of an orbit, on several hundred planar Newtonian periodic three-body orbits. Here $T$ is the period of an orbit, $E$ is its energy, so that $T_{\rm s.i.}$ is the scale-invariant (s.i.) period, or, equivalently, the period at unit energy $|E| = 1$. All of these orbits have vanishing angular momentum and pass through a linear, equidistant configuration at least once. Such orbits are classified in ten algebraically well-defined sequences. Orbits in each sequence follow an approximate linear dependence of $T_{\rm s.i.}$, albeit with slightly different slopes and intercepts. The orbit with the shortest period in its sequence is called the "progenitor": six distinct orbits are the progenitors of these ten sequences. We have studied linear stability of these orbits, with the result that 21 orbits are linearly stable, which includes all of the progenitors. This is consistent with the Birkhoff-Lewis theorem, which implies existence of infinitely many periodic orbits for each stable progenitor, and in this way explains the existence and ensures infinite extension of each sequence.
• ### Field-only integral equation method for time domain scattering of electromagnetic pulses(1708.07934)

The scattering of electromagnetic pulses is described using a non-singular boundary integral method to solve directly for the field components in the frequency domain, and Fourier transform is then used to obtain the complete space-time behavior. This approach is stable for wavelengths both small and large relative to characteristic length scales. Amplitudes and phases of field values can be obtained accurately on or near material boundaries. Local field enhancement effects due to multiple scattering of interest to applications in microphotonics are demonstrated.
• ### A Robust Multi-Scale Field-Only Formulation of Electromagnetic Scattering(1611.01978)

We present a boundary integral formulation of electromagnetic scattering by homogeneous bodies that are characterized by linear constitutive equations in the frequency domain. By working with the Cartesian components of the electric, E and magnetic, H fields and with the scalar functions (r*E) and (r*H), the problem is cast as solving a set of scalar Helmholtz equations for the field components that are coupled by the usual electromagnetic boundary conditions at material boundaries. This facilitates a direct solution for E and H rather than working with surface currents as intermediate quantities in existing methods. Consequently, our formulation is free of the well-known numerical instability that occurs in the zero frequency or long wavelength limit in traditional surface integral solutions of Maxwell's equations and our numerical results converge uniformly to the static results in the long wavelength limit. Furthermore, we use a formulation of the scalar Helmholtz equation that is expressed as classically convergent integrals and does not require the evaluation of principal value integrals or any knowledge of the solid angle. Therefore, standard quadrature and higher order surface elements can readily be used to improve numerical precision. In addition, near and far field values can be calculated with equal precision and multiscale problems in which the scatterers possess characteristic length scales that are both large and small relative to the wavelength can be easily accommodated. From this we obtain results for the scattering and transmission of electromagnetic waves at dielectric boundaries that are valid for any ratio of the local surface curvature to the wave number. This is a generalization of the familiar Fresnel formula and Snell's law, valid at planar dielectric boundaries, for the scattering and transmission of electromagnetic waves at surfaces of arbitrary curvature.
• ### Non-singular field-only surface integral equations for electromagnetic scattering(1608.02058)

A boundary integral formulation of electromagnetics that involves only the components of $\boldsymbol{E}$ and $\boldsymbol{H}$ is derived without the use of surface currents that appear in the classical PMCHWT formulation. The kernels of the boundary integral equations for $\boldsymbol{E}$ and $\boldsymbol{H}$ are non-singular so that all field quantities at the surface can be determined to high precision and also geometries with closely spaced surfaces present no numerical difficulties. Quadratic elements can readily be used to represent the surfaces so that the surface integrals can be calculated to higher numerical precision than using planar elements for the same numbers of degrees of freedom.
• ### Predictions of homogeneous nucleation rates for $n$-alkanes accounting for the diffuse phase interface and capillary waves(1705.07588)

Homogeneous droplet nucleation has been studied for almost a century but has not yet been fully understood. In this work, we used the density gradient theory (DGT) and considered the influence of capillary waves (CW) on the predicted size-dependent surface tensions and nucleation rates for selected $n$-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) and compared to the classical nucleation theory and the Peng--Robinson EoS. It was found that the critical clusters are practically free of CW because they are so small that even the smallest CW wavelengths do not fit into their finite dimensions. The CW contribute to the entropy of the system and thus decrease the surface tension. A correction for the effect of CW on the surface tension is presented. The effect of the different EoSs is relatively small because by a fortuitous coincidence their predictions are similar in the relevant range of critical cluster sizes. The difference of the DGT predictions to the classical nucleation theory computations is important but not decisive. Of the effects investigated, the most pronounced is the suppression of the CW which causes a sizable decrease of the predicted nucleation rates. The major difference between experimental nucleation rate data and theoretical predictions remains in the temperature dependence. For normal alkanes, this discrepancy is much stronger than observed, e.g., for water. The theoretical corrections developed here have a minor influence on the temperature dependency. We provide empirical equations correcting the predicted nucleation rates to values comparable with experiments.
• ### Application of axiomatic formal theory to the Abraham--Minkowski controversy(1502.05997)

We treat continuum electrodynamics as an axiomatic formal theory based on the macroscopic Maxwell--Minkowski equations applied to a thermodynamically closed system consisting of an antireflection-coated block of a simple linear dielectric material situated in free-space that is illuminated by a quasimonochromatic field. We prove that valid theorems of the formal theory of Maxwellian continuum electrodynamics are inconsistent with conservation laws for the inviscid incoherent flow of non-interacting particles (photons) in the continuum limit (light field) in the absence of external forces, pressures, or constraints. We also show that valid theorems of Maxwellian continuum electrodynamics are contradicted by the refractive index-independent Lorentz factor of von Laue's application of Einstein's special relativity to a dielectric medium. Obviously, the fundamental physical principles in the vacuum are not affected. However, the extant theoretical treatments of electrodynamics, special relativity, and energy--momentum conservation must be regarded as being mutually inconsistent in a simple linear dielectric in which the effective speed of light is $c/n$. Having proven that the established applications of fundamental physical principles to dielectric materials lead to mutually inconsistent theories, we derive, from first principles, a mutually consistent, alternative theoretical treatment of electrodynamics, special relativity, and energy--momentum conservation in an isotropic, homogeneous, linear dielectric-filled, flat, non-Minkowski, continuous material spacetime.
• ### Accurate and Efficient Evaluation of Characteristic Modes(1709.09976)

A new method to improve the accuracy and efficiency of characteristic mode (CM) decomposition for perfectly conducting bodies is presented. The method uses the expansion of the Green dyadic in spherical vector waves. This expansion is utilized in the method of moments (MoM) solution of the electric field integral equation to factorize the real part of the impedance matrix. The factorization is then employed in the computation of CMs, which improves the accuracy as well as the computational speed. An additional benefit is a rapid computation of far fields. The method can easily be integrated into existing MoM solvers. Several structures are investigated illustrating the improved accuracy and performance of the new method.
• ### Dynamics of non-holonomic systems with stochastic transport(1707.04741)

This paper formulates a variational approach for treating observational uncertainty and/or computational model errors as stochastic transport in dynamical systems governed by action principles under nonholonomic constraints. For this purpose, we derive, analyze and numerically study the example of an unbalanced spherical ball rolling under gravity along a stochastic path. Our approach uses the Hamilton-Pontryagin variational principle, constrained by a stochastic rolling condition, which we show is equivalent to the corresponding stochastic Lagrange-d'Alembert principle. In the example of the rolling ball, the stochasticity represents uncertainty in the observation and/or error in the computational simulation of the angular velocity of rolling. The influence of the stochasticity on the deterministically conserved quantities is investigated both analytically and numerically. Our approach applies to a wide variety of stochastic, nonholonomically constrained systems, because it preserves the mathematical properties inherited from the variational principle. Keywords: Nonholonomic constraints, Stochastic dynamics, Transport noise.
• ### On Noether's theorem for the Euler-Poincar\'e equation on the diffeomorphism group with advected quantities(1206.2976)

We show how Noether conservation laws can be obtained from the particle relabelling symmetries in the Euler-Poincar\'e theory of ideal fluids with advected quantities. All calculations can be performed without Lagrangian variables, by using the Eulerian vector fields that generate the symmetries, and we identify the time-evolution equation that these vector fields satisfy. When advected quantities (such as advected scalars or densities) are present, there is an additional constraint that the vector fields must leave the advected quantities invariant. We show that if this constraint is satisfied initially then it will be satisfied for all times. We then show how to solve these constraint equations in various examples to obtain evolution equations from the conservation laws. We also discuss some fluid conservation laws in the Euler-Poincar\'e theory that do not arise from Noether symmetries, and explain the relationship between the conservation laws obtained here, and the Kelvin-Noether theorem given in Section 4 of Holm, Marsden and Ratiu, {\it Adv. in Math.}, 1998.
• ### Amorphous topological insulators constructed from random point sets(1612.09267)

The discovery that the band structure of electronic insulators may be topologically non-trivial has unveiled distinct phases of electronic matter with novel properties. Recently, mechanical lattices have been found to have similarly rich structure in their phononic excitations, giving rise to protected uni-directional edge modes whose existence was demonstrated in lattices of interacting gyroscopes and coupled pendula. In all these cases, however, as well as in other topological metamaterials, the underlying structure was finely tuned, be it through periodicity, quasi-periodicity or isostaticity. Here we show that amorphous mechanical Chern insulators consisting of interacting gyroscopes can be readily constructed from arbitrary underlying structures, including hyperuniform, jammed, quasi-crystalline, and uniformly random point sets. While our findings apply to mechanical and electronic systems alike, we focus on networks of interacting gyroscopes as a model system. Local decoration control the topology of the vibrational spectrum, endowing amorphous structures with protecting edge modes -- with a chirality of choice. Using a real-space generalization of the Chern number, we investigate the topology of our structures numerically, analytically and experimentally. The robustness of our approach enables the topological design and self-assembly of non-crystalline topological metamaterials on the micro and macro scale.
• ### On the nature of discrete space-time: The atomic theory of space-time and its effects on Pythagoras's theorem, time versus duration, inertial anomalies of astronomical bodies, and special relativity at the Planck scale(1608.08506)

In this work, resolutions will be given for commonly stated problems associated with a model that assumes that space and time are discretized (i.e., atomized). This model is in contrast to the continuous space-time model that is used in all common physical theories and equations -- a model that assumes that spatial coordinates and time are continuous variables. The resolutions to the problems are arrived at, not by proposing any new theories or postulates, but by strictly adhering to: Ernst Mach's principle of non-absolute space, the tenets of logical positivism, quantum mechanics and general relativity. The problems associated with discrete space-time addressed in this paper include: Lorentz contraction (time dilation) of the ostensibly smallest spatial (temporal) interval, maintaining isotropy, violations of causality, and conservation of energy and momentum. Importantly, this work yields modifications to the standard formulae for time dilation and length contraction, with these modifications preserving the quantums of space and time and allowing for temporary travel at the speed of light. Also given are: a resolution to Weyl's tile argument, a modification of the 2500 year old Pythagoras's theorem, a reassessment of Henri Bergson's theory advocating a distinction between the time durations measured by scientists and an immutable "Time", and a discussion of whether Einstein's light-clocks are as ideal as most scientist believe. Also included is a demonstration of how discrete space imposes order upon John Wheeler's quantum foam such that the foam becomes a gravity crystal permeating all space and producing measurable inertial anomalies of astronomical bodies.
• ### High order moment model for polydisperse evaporating sprays towards interfacial geometry(1608.07148)

In this paper we propose a new Eulerian modeling and related accurate and robust numerical methods, describing polydisperse evaporating sprays, based on high order moment methods in size. The main novelty of this model is its capacity to describe some geometrical variables of the droplet-gas interface, by analogy with the liquid-gas interface in interfacial flows. For this purpose, we use fractional size-moments, where the size variable is taken as the droplet surface. In order to evaluate the evaporation of the polydisperse spray, we use a smooth reconstruction which maximizes the Shannon entropy. However, the use of fractional moments introduces some theoretical and numerical difficulties, which need to be tackled. First, relying on a study of the moment space, we extend the Maximum Entropy (ME) reconstruction of the size distribution to the case of fractional moments. Then, we propose a new accurate and realizable algorithm to solve the moment evolution due to evaporation, which preserves the structure of the moment space. This algorithm is based on a mathematical analysis of the kinetic evolution due to evaporation, where it shown that the evolution of some negative order fractional moments have to be properly predicted, a peculiarity related to the use of fractional moments. The present model and numerical schemes yield an accurate and stable evaluationof the moment dynamics with minimal number of variables, as well as a minimal computational cost as with the EMSM model, but with the very interesting additional capacity of coupling with diffuse interface model and transport equation of averaged geometrical interface variables, which are essential in oder to describe atomization.
• ### Quantum mechanics as a solution to the classical self-force problem(0902.4606)

It is argued that, contrary to conventional wisdom, no trustworthy universal self-force/radiative corrections to the Lorentz force equation, can be derived from the basic tenets of classical electrodynamics. This concords with the apparent randomness observed in quantum mechanical scattering experiments and with the absence of any experimental support for such universality. In a recent paper [11], the statistical effect of radiative corrections to the motion of charged bodies has been derived from the basic tenets and does take a universal form, described by quantum mechanical wave equations---again conforming with experiment. As that derivation assumes nothing about the size, mass or composition of the body, it is conjectured that quantum mechanics is the appropriate framework for dealing also with radiative corrections to the motion of macroscopic bodies.
• ### Novel Practicable GLHUA Invisible Cloak Absorbs the Incident Wave and creates Outgoing Wave Without Exceeding Light Speed Propagation(1706.10147)

In this paper, we propose novel practicable GLHUA invisible cloak, our cloak absorbs the incident wave and creates outgoing wave that make it to be invisible cloak and analytical electromagnetic wave without exceeding light speed propagation. Ten more discoveries and creations are reported. 1. New GLHUA annular layer electromagnetic (EM) invisible cloak with relative parameter not less than 1. 2. Analytical EM wave solution of EM full wave equation with the GLHUA 1-3 cloak materials is found, In annular layer cloak,$R_1 \le r \le R_2$, $R_1>0$,$R_2>R_1$, $\varepsilon_r = \mu_r = \frac{{R_2 ^2 }}{{r^2 }}$, $\varepsilon_\theta = \varepsilon_\phi = \mu_\theta = \mu_\phi = \frac{{R_2 - R_1 }}{{r - R_1 }}$. 3. A GLHUA expansion method to find an exact analytical EM wave propagation in the GLHUA double cloak is propsed. 4.Novel negative space is proposed.5. GLHUANP transformation is proposed, which maps negative infinity,$- \infty$ to $R_1$ and $-R_2$ to $R_2$. 6. The GLHUANP transformation creates GLHUA 2 and GLHUA3 cloaks. the 7.We find analytical EM wave without exceeding light speed propagation in the cloak. 8. GLHUA 1 and GLHUA 2 and GLHUA3 cloaks absorb the incident wave and create outgoing wave that make them invisible cloak and analytical electromagnetic wave through them. 9. GLHUAF transformation and .GLHUAF invisible cloak with double negative materials is created. 11. The relative parameters of the GLHUA1 are large than 1, radial parameter in GLHUA 3 cloak is large than zero, refractive index of the GLHUA 1-3 cloak is large than 1. 12. The wave front and ray tracing in GLHUA cloak are discontinuous, the ray redirection can not be used for GLHUA invisible cloak. In my arXiv:1612.02857, I refer my this paper,that is not overlap. The arXiv admin note is made by its software. The patent and copyright and all rights are reserved by authors in GL Geophysical Laboratory in USA.
• ### The Theoretical Proof for GLHUA EM Invisible Double Layer Cloak By Using GL No Scattering Modeling and Inversion(1701.00534)

In this paper, we proposed the GLHUA double layer cloak; proved the properties of GLHUA double layer cloak; Using GL no scattering inversion and the pre cloak condition 6.1 to 6.4 in paper [1], we create GLHUA outer layer cloak radial relative parameter and angular relative parameter theoretically. We proved theorem 4.1 to theorem 4.4 that the phase velocity of the electromagnetic wave in GLHUA outer layer cloak is less than light speed and tends to zero at the boundary $r=R_1$, the EM wave tends to zero at $r=R1$; When source $r_s > R_2$ in the outside of the cloak and observer $r_o < R_o$ in the concealment, then $\vec E(\vec r_o)=0$ and $\vec H(\vec r_o)=0$ that is proved in theorem 5.1 and 5.3. We prove that the EM wave excited in outside of the cloak can not propagation penetrate into the concealment. In theorem 5.4 to 5.6, we rigorously proved the EM wave excited in outside of the cloak can not be disturbed by the cloak. In theorem 6.1 to theorem 6.6, we prove that the EM wave excited in the concealment can not propagate to outside of GLHUA inner cloak. we prove that the EM wave excited in the concealment can not be disturbed by the cloak. We theoretically prove that the GLHUA double cloak is invisible cloak with concealment and with relative parameter not less than 1; the GLHUA double layer cloak is practicable. In December of 2016, we submitted 3 papers to arXiv for GLHUA double layer cloak with relative parameters not less than 1. This is the third paper. The task and content of the three paper are different from each other. The second paper [1] and this Paper are theoretical base and proof of the subjects in paper arXiv:1612.02857. We find an exact analytical EM wave field of Maxwell EM equation in GLHUA cloak and mirage bridge wave. Patent of the GLHUA EM cloaks,GLHUA sphere and GL modeling and inversion methods are reserved by authors in GL Geophysical Laboratory.
• ### Electric field and electric forces in a spontaneously polarized nonpolar isotropic dielectric(1702.04290)

Based on the microscopic Maxwell equations, we develop a method of description of the electric field in a spontaneously polarized isotropic nonpolar dielectric. We find the solution for the electric field $\textbf{E}(\textbf{r})$ for several typical examples. Moreover, we generalize Helmholtz's formula for the electric force acting on a volume element of a dielectric with regard for the contribution of the spontaneous polarization.