Superposed Dipole Overview

• Superposed Dipole is a unified configuration formed by the linear addition of dipole sources and fields, used to represent both ideal point dipoles and effective dipoles.
• It employs Green’s function formulations, image-dipole lattices, and distributional methods to satisfy boundary conditions in layered, cylindrical, and bounded media.
• The approach enhances numerical stability and accuracy in electromagnetic modeling, offering practical insights for thin films, superconducting traps, and curved-spacetime scenarios.

A superposed dipole is a dipolar configuration defined by linear addition of dipole sources, dipole fields, or dipole moments. In the cited literature, the term covers several closely related constructions: co-located ideal point dipoles treated as distributions, a physical dipole plus its image lattice in bounded media, superpositions of Hertzian dipoles through dyadic Green’s functions in layered or stratified environments, effective displaced dipoles used to represent interacting material bodies, and surface dipole layers on hypersurfaces in curved spacetime (Parker, 2016, Headley, 29 Jun 2026, Cho et al., 2016, Gürlebeck et al., 2011). This suggests a family of mathematically linear but physically distinct objects, unified by the fact that Maxwell’s equations and the associated source representations remain additive.

1. Linear structure and source representation

The basic analytic structure is the Green-function relation between a dipole source and the resulting field. In layered media, the electric field of a point electric dipole at $\mathbf{r}'$ with dipole moment $\mathbf{p}$ is written as

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$

and for multiple dipoles,

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$

This is the explicit superposition formula used for dipole sources in layered media (Cho et al., 2016).

The same principle appears in thin-film electromagnetics. For an arbitrary dipole orientation at a conducting sheet,

$$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$

the field is expressed as a linear combination of the fields of a horizontal electric dipole and a vertical electric dipole. In that setting, the total current of several dipoles is

$$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$

and the total field is the corresponding sum of Green-tensor contributions (Margetis et al., 2015).

In cylindrically stratified media, the same linearity is used for arbitrarily oriented Hertzian dipoles with current density

$$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$

There the computational problem is not the validity of superposition but its stable evaluation across extreme layer resistivities, frequencies, and source-observer separations; the solution is a pseudoanalytical tensor Green’s function built from range-conditioned cylindrical functions and adaptively chosen integration contours (Moon et al., 2013).

2. Ideal point dipoles and distributional superposition

For the ideal electric dipole at the origin, the full field is not just the familiar $1/r^3$ far field. Under spherical regularization, the complete distributional field is

$$\mathbf{E}(\mathbf{r})= \frac{1}{r^3}\big[3(\mathbf{p}\cdot\hat{\mathbf{r}})\hat{\mathbf{r}}-\mathbf{p}\big] -\frac{4\pi}{3}\,\mathbf{p}\,\delta^{(3)}(\mathbf{r}),$$

and may also be written as

$$E_i(\mathbf{x}) = p_j\,\partial_i\partial_j\!\left(\frac{1}{r}\right), \qquad \partial_i\partial_j\!\left(\frac{1}{r}\right) = \frac{3n_in_j-\delta_{ij}}{r^3} -\frac{4\pi}{3}\delta_{ij}\delta^{(3)}(\mathbf{x}).$$

The source density is

$\mathbf{p}$0

and the divergence of the far-field term alone is

$\mathbf{p}$1

A central result is that the divergence of a singular vector field can contain a derivative of a Dirac delta function even if the field itself does not contain a delta function (Parker, 2016).

This distributional structure resolves the apparent paradox produced by a naive application of Gauss’s law. For the ideal dipole, two thirds of the total $\mathbf{p}$2 comes from the divergence of the far field and one third from the explicit delta term in $\mathbf{p}$3. The far field therefore cannot be discarded at the origin when reconstructing the source.

For superposed ideal dipoles at the same point, linearity is exact at the level of both fields and distributions. If

$\mathbf{p}$4

then

$\mathbf{p}$5

and

$\mathbf{p}$6

For dipoles at distinct points $\mathbf{p}$7,

$\mathbf{p}$8

with the corresponding sum of translated dipole fields. The same framework extends to higher multipoles, for which the divergence contains higher derivatives of $\mathbf{p}$9 (Parker, 2016).

3. Image-dipole lattices in bounded geometries

A particularly concrete superposed-dipole construction is the image lattice of a magnetic dipole inside a closed cuboidal superconducting trap. The cavity interior is

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$0

and the Meissner boundary condition is

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$1

Starting from a point dipole $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$2 at $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$3, repeated reflection in the six walls generates an infinite $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$4D lattice of images indexed by $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$5 with

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$6

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$7

The total field is the sum over this lattice, and the paper proves that the resulting superposition satisfies the Meissner boundary condition on all six walls simultaneously (Headley, 29 Jun 2026).

For a centered dipole, the image self-energy is

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$8

and the orientational part reduces to the diagonal quadratic form

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\,\mathbf{G}(\mathbf{r},\mathbf{r}')\cdot\mathbf{p},$$9

The coefficients $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$0 are Epstein-zeta-type lattice sums; off-diagonal terms vanish by symmetry. This makes the equilibrium orientation an axis-selection problem controlled entirely by the image-dipole superposition.

The main qualitative result is that in both infinite and finite rectangular traps the dipole orientation aligns with the short cross-sectional axis over a finite range of aspect ratios. In the tube limit, the sign change occurs at

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$1

and, by symmetry, at

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$2

so that for

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$3

the short cross-sectional axis is preferred. In the fully closed cuboid, the phase diagram in the $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$4 plane contains regions where the easy axis is $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$5, $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$6, or $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$7, as well as four triple points, including the cube and the accidental symmetry orbit

$$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$8

Finite-element verification agrees with the image-lattice predictions to better than $$\mathbf{E}(\mathbf{r}) = i\omega\mu_0\sum_m \mathbf{G}(\mathbf{r},\mathbf{r}_m)\cdot\mathbf{p}_m.$$9 (Headley, 29 Jun 2026).

4. Planar, layered, and cylindrical media

In planar thin-film electromagnetics, superposed dipoles are used as both source decompositions and mode-selective exciters. A vertical electric dipole predominantly excites TM modes, whereas a horizontal electric dipole excites both TM and TE modes. Because the field of an arbitrary dipole orientation is a linear combination of the vertical and horizontal solutions, the TM and TE content of the resulting field can be engineered by source orientation and phase (Margetis et al., 2015).

In layered media, the dyadic Green’s function naturally separates into a singular free-space part and a nonsingular scattering part,

$$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$0

which is especially useful for superposed dipoles because the singular interaction can be treated with standard free-space tools while layered effects appear as an additive correction. The practical advance of the layered-media formulation is that all six independent dyadic components are reduced to a small set of Sommerfeld integrals by Bessel identities, making many-dipole computations significantly more efficient (Cho et al., 2016).

In cylindrically stratified media, the same superposition principle is retained, but numerical stability becomes dominant. Direct evaluation of the canonical cylindrical-wave formulas leads to underflow and overflow because of poor scaling of Bessel and Hankel functions for extreme arguments and high azimuthal order. The stable formulation introduces range-conditioned modified cylindrical functions for small, moderate, and large arguments, together with Sommerfeld Integration Path and Deformed Sommerfeld Integration Path contours in the complex spectral-wavenumber plane. The method remains stable for layer resistivities ranging from $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$1m to $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$2m and for frequencies from $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$3 MHz down to $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$4 Hz (Moon et al., 2013).

These planar and cylindrical formulations are not merely numerical conveniences. They show that a superposed dipole can be represented equivalently as a discrete source ensemble, a modal expansion, or a Green-tensor convolution, depending on the geometry and boundary conditions.

5. Effective and displaced dipoles in material systems

A different use of the term appears when extended bodies are reduced to effective dipoles whose fields are then superposed. For two identical, electrically conducting, permeable spheres in an oscillating homogeneous magnetic field, each sphere is represented by one effective magnetic dipole,

$$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$5

but the dipole origin is displaced from the sphere center and the dipole moment is complex. For a single sphere,

$$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$6

so amplitude and phase are both encoded in the response factor. In the two-sphere problem, the total field is

$$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$7

and the moments and displacements are updated by coupled fixed-point equations (Bönsel et al., 2021).

The displacement is not a minor correction. For small surface-to-surface distance $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$8, the paper reports that $$\mathbf{p}=p_x\hat{\mathbf{x}}+p_y\hat{\mathbf{y}}+p_z\hat{\mathbf{z}},$$9 reaches up to approximately $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$0 of sphere diameter, and that neglecting displacement leaves noticeable phase error. With the full Interacting Displaced Dipole model, a tolerance $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$1 in amplitude, phase, and displacement is reached with $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$2 iterations for typical parameters; in some cases the phase error is reduced to approximately $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$3 and the relative amplitude error to approximately $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$4 (Bönsel et al., 2021).

A depth-resolved version of the same idea appears in electromagnetic induction over horizontally stratified earth. There the primary magnetic dipole excites eddy currents in thin horizontal layers, and each layer contributes a secondary field that behaves as an effective dipole-like source. The total response is the superposition of layer contributions. The damped forward model includes interaction between eddy currents through a conductivity-dependent damping factor and improves the validity range relative to the McNeill low-induction-number model: for a $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$5 error criterion, the induction-number limit increases from $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$6 to $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$7, and for $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$8 error from $$\mathbf{J}(\mathbf{r})=\sum_i \mathbf{p}_i\,\delta(\mathbf{r}-\mathbf{r}_i),$$9 to $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$0 (Delrue et al., 2017). This suggests a superposed-dipole interpretation of conductive layering in which each depth interval acts as a secondary dipolar scatterer.

6. Surface dipole layers, exact source equivalences, and controversy

In curved spacetime, superposed dipoles are formalized as distributions on timelike hypersurfaces. Using the complex self-dual Maxwell tensor

$$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$1

and the complex current

$$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$2

Maxwell’s equations become

$$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$3

An electric dipole layer has a distributional current proportional to $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$4 on the hypersurface $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$5, and its jumps are

$$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$6

The formalism is explicitly linear, so monopole layers, electric dipole layers, magnetic dipole layers, and their combinations are superposed by summing their surface currents and the corresponding jumps of $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$7 and $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$8 (Gürlebeck et al., 2011).

A central exact equivalence relates magnetic dipole layers to electric surface currents: $$\mathbf{J}(\mathbf{r}) = Il\,\hat{\alpha}'\,\delta(\mathbf{r}-\mathbf{r}').$$9 Thus a magnetic dipole layer and a suitable electric monopole current layer generate the same external field. In the Schwarzschild-disk constructions, the literature explicitly combines electric monopole layers with magnetic dipole layers, or electric dipole layers with magnetic monopole layers, by suitable choices of the field above and below the disk (Gürlebeck et al., 2011).

A distinct and explicitly controversial usage appears in “Linear Superposition Effect at Sources and in Waves” (Jiao, 25 Aug 2025). There a superposed dipole is defined as the effective radiation dipole formed by two identical, co-phase, closely spaced radiating dipoles: $1/r^3$0 Under the assumptions

$1/r^3$1

the paper assigns the effective radiation power

$1/r^3$2

The same summary states that this interpretation is at odds with conventional electromagnetic theory and conventional energy conservation, and notes that the cited experiment observed a $1/r^3$3-fold increase rather than a full factor of $1/r^3$4 (Jiao, 25 Aug 2025). Within the broader literature on superposed dipoles, this usage therefore occupies a controversial and non-mainstream position.

Superposed dipole is therefore best understood not as a single invariant object but as a recurrent electromagnetic construction. In distribution theory it is the exact sum of singular ideal dipoles; in bounded media it is a dipole plus its image system; in Green-function methods it is an ensemble of Hertzian sources; in effective-medium models it is a reduced description of interacting bodies or layers; and in curved-spacetime layer theory it is a sum of distributional monopole and dipole currents on hypersurfaces (Parker, 2016, Headley, 29 Jun 2026, Cho et al., 2016, Bönsel et al., 2021, Gürlebeck et al., 2011).