Pith. sign in

REVIEW 4 major objections 4 minor 46 references

Full Vectorial Maxwell Equations with Continuous Angular Indices

T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper argues that Maxwell's equations in spherical geometry admit solutions with continuous real angular indices $(\ell,m)$ whose singular fields carry finite electromagnetic energy exactly when $\ell > -1/2$.

desk verdict A bold but unproven continuous-index Maxwell framework whose central ℓ > −1/2 energy threshold rests on an asserted regularizing function and is contradicted by the paper's own asymptotics; worth referee time but not publication as it stands. read the letter →

arxiv 2508.02675 v1 pith:U6XJYVVT submitted 2025-07-10 math.NA cs.NAmath-phmath.MPphysics.class-ph

Mustafa Bakr This is my paper · ORCID
classification math.NAcs.NAmath-phmath.MPphysics.class-ph MSC 35Q6133C5546E3535B40
keywords continuousangularspectrumsingularelectromagneticmodesfiniteenergycriterionMaxwellequationssphericalgeometryweightedSobolevspacesspectralintegralsnon-separablevector
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish that Maxwell's equations in spherical geometry still have physically meaningful solutions when the angular momentum indices $\ell$ and $m$ are allowed to be any real numbers, not just integers. Such solutions are singular at the origin, with radial scaling $E_r \sim r^{\ell-1}$, yet the paper argues their electromagnetic energy stays finite exactly when $\ell > -1/2$. If true, this gives a new continuous family of admissible singular field configurations beyond the standard integer spherical-harmonic basis, with possible applications to wedge-like cavities, field concentration, and non-periodic azimuthal domains. The argument is built on a continuous spectral integral over non-integer modes, coupled vector equations, and weighted Sobolev spaces.

What carries the argument

The load-bearing machinery is a continuous-index replacement for spherical harmonics: generalized functions $\Phi_{\ell m}(\theta,\phi)$ built from analytically continued associated Legendre functions $P^{|m|}_{\ell}(\cos\theta)$, a spectral weight $w(\ell,m) = \pi\,\Gamma(\ell+|m|+1)/(\sin(\pi(\ell-|m|))\,\Gamma(\ell-|m|+1))$, and dual functions $\tilde{\Psi}_{\ell m}$ involving $Q^{|m|}_\ell$ that give a delta-normalized biorthogonal system. Fields are expanded as spectral integrals $\mathbf{E}(r)=\int a(\ell,m)\,r^{\alpha(\ell,m)}\boldsymbol{\Phi}_{\ell m}\,d\ell\,dm$, with the singular exponent tied to the angular operator eigenvalue by $\alpha(\ell,m)=\tfrac12(\sqrt{1+4\lambda_{\ell m}}-1)$. The coupled angular components are reconstructed through Green's functions built from spherical Bessel functions of non-integer order $\nu=\sqrt{\ell(\ell+1)-1}$, and the energy-convergence argument rests on a regularizing function $F(\ell,m;\xi)=1+O(\xi^2)$ with $\xi=r/\sin\theta$. This machinery converts an apparently divergent angular integral into a convergent energy integral.

What would settle it

Compute the coupled Green's function integrals (Eqs. (153)--(159)) at $r\to 0$ for $\ell=-0.25$, $m=0.5$ without inserting the asserted regularization; if the resulting energy integral diverges, the claim that $\ell>-1/2$ suffices is false.

Watch

Extended reading notes

Core claim

The central claim is that the full vectorial Maxwell equations, not an approximate or separable version, admit solutions labelled by continuous real angular indices $(\ell,m)$, with fields of the form $\mathbf{E} = \int a(\ell,m)\, r^{\alpha(\ell,m)} \boldsymbol{\Phi}_{\ell m}(\theta,\phi)\, d\ell\, dm$. Near $r=0$ the components scale as $E_r \sim r^{\ell-1}$, $E_\theta \sim \ell/r$, and $E_\varphi \sim m/(r^3\sin^2\theta)$; the apparent divergences are controlled because the curl coupling makes the angular and radial parts non-separable. The paper proves that the integrated energy converges precisely for $\ell > -1/2$ for all $m$, and that the fields belong to weighted Sobolev spaces $H^s_w(\Omega)$ with $s < \ell + 1/2$. It constructs the spectral kernels and biorthogonal function systems that make the continuous expansion explicit, and verifies the singularity exponents numerically by Galerkin projection.

Load-bearing premise

The energy-convergence criterion $\ell>-1/2$ rests on a correction factor $F(\ell,m;\xi)$ that the paper asserts is $1+O(\xi^2)$ near the origin, with the derivation deferred to an appendix that is not present; if that factor cannot be derived from the coupled Maxwell equations, the finite-energy claim for singular modes does not follow.

Editorial extensions

If this is right

  • Fields with $\ell$ in $(-1/2,0)$ are singular at the origin yet carry finite energy, so they are admissible in the energy norm even though they are not square-integrable in the usual $L^2$ sense.
  • In a spherical cavity with azimuthal span $\Phi_0<2\pi$, allowed azimuthal indices are $m=n\Phi_0/2\pi$, and $\ell$ becomes coupled to $m$; the continuous spectrum replaces the discrete spherical-harmonic ladder.
  • Truncating the continuous spectral integral converges: $\|f-f_N\|_{L^2(S^2)} \le C\|f\|_{H^s} N^{-s+1/2+\epsilon}$ for $s>1/2$, so the expansion is computationally usable.
  • Resonances appear as poles of the spectral coefficient $a(\ell,m)$ in the complex $\ell$ plane, with lifetimes encoded by the imaginary parts, extending quasi-normal-mode ideas to continuous angular indices.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the decisive check is the regularizing function $F(\ell,m;\xi)$; until its derivation is supplied, the $\ell>-1/2$ criterion should be read as conditional on that function's existence.
  • Editorial inference: if the criterion holds, a conical or wedge cavity whose azimuthal span is not $2\pi$ should show enhanced, singular field concentration near the apex for the lowest continuous mode, which could be probed numerically or experimentally.
  • Editorial inference: the same analytic continuation may transfer to acoustic or quantum wave problems with broken rotational symmetry, where continuous 'angular momentum' singular modes could be constructed by the same biorthogonal Green's function route.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper claims to extend Maxwell's equations in spherical geometry to a continuous spectral representation with real angular indices (ℓ,m), going beyond the discrete spherical-harmonic basis. It constructs a biorthogonal system of generalized vector spherical functions with a spectral weight w(ℓ,m), derives coupled integral equations for the angular field components, and states that the electromagnetic energy converges exactly when ℓ > -1/2. The manuscript also reports Galerkin and spectral-integral numerics, including an optimization procedure for the spectral weight parameters.

Significance. If the central claim were established, it would represent a new class of admissible singular electromagnetic modes beyond the standard spherical-harmonic basis, with potential applications to wedge geometries, broken azimuthal symmetry, and singular cavities. The paper's strengths include explicit formulas for the spectral weight (Eq. (69)), the Green's function (Eq. (106)), and a detailed numerical optimization procedure (Section V). However, the manuscript does not provide proofs for the completeness of the continuous biorthogonal system, the contraction estimate for the coupled integral equations, or the regularizing function F that underlies the energy criterion; one of its own asymptotic summaries (Eq. (210)) contradicts the regularized scaling used to obtain finite energy. At present the significance of the claimed result is not assessable from the submitted text.

major comments (4)
  1. [Sec. IV.C, Eqs. (189)-(192) and (210)] The energy criterion ℓ > -1/2 is derived by inserting the regularizing function F(ℓ,m;ξ) with F = 1 + O(ξ^2) into the angular-component scaling. The derivation of F is deferred to 'Appendix A', which is not present in the manuscript, and the paper's own summary Eq. (210) states Eθ ∼ ℓ/r and Eφ ∼ m/(r^3 sin^2θ). The Eθ scaling alone makes the energy integral ∫|Eθ|^2 r^2 dr diverge linearly for every nonzero ℓ, so the finite-energy claim is internally inconsistent unless Eq. (210) is explicitly identified as an incorrect separable approximation and the true non-separable asymptotics for all three components are derived. Because the ℓ > -1/2 criterion depends entirely on F, this is a load-bearing gap.
  2. [Sec. III.F, Eqs. (60) and (69)-(72)] The completeness of the continuous-index biorthogonal system is asserted without proof. The weight w(ℓ,m) = πΓ(ℓ+|m|+1)/(sin(π(ℓ-|m|))Γ(ℓ-|m|+1)) diverges at integer ℓ-|m|, and the manuscript states that one recovers usual orthonormality at integer limits but gives no limiting procedure. On the compact sphere, the standard spherical harmonics already form a complete set, so a continuous-index basis cannot simply be added without a theorem showing that the combined system is a resolution of the identity. The completeness relations (70) and (72) are used for arbitrary field expansions and must be proved.
  3. [Sec. III.G.6, Eq. (126)] The contraction mapping argument for the coupled integral equations is asserted rather than proved. The bound γ ≤ max{|Cθφ|, |Cθφ*|} · sup|G| · (integration bounds) is not a theorem: no function space is specified, and the Green's function (106) together with the source terms r'^{-2} d/dr'[r'^2 Eφ] requires control of derivatives and singular integrals that is not supplied. Since existence and uniqueness of the angular components is a prerequisite for the spectral construction, this is a load-bearing step that needs a rigorous estimate.
  4. [Sec. V.C.3, Eqs. (288)-(294)] The numerical validation of finite energy is circular. The spectral weight parameters A, p, q, β are optimized subject to the constraint p > 3/2 - ℓ_min, which is exactly the condition obtained in Eq. (292) from the same energy-convergence analysis being tested. The finite-energy result shown in Figure 1(d) is therefore generated under the very constraint derived from the claim under validation, and it does not provide an independent confirmation of the ℓ > -1/2 criterion.
minor comments (4)
  1. [Sec. IV.D, Eq. (202)] The conversion of ℓ > -1/2 into λ > -3/4 is not correct as stated: with α(α+1) = λ and α = ℓ - 1, the positive branch α = (-1 + sqrt(1+4λ))/2 is ≥ -1/2 for every λ ≥ -1/4, so the inequality λ > -3/4 is automatically satisfied and does not encode the admissibility threshold. The branch analysis must be redone.
  2. [Sec. IV.E, Eq. (211)] The summary bullet 'Energy convergence requires ℓ > -1/2 and m = 0' contradicts the criterion in Eq. (195) (ℓ > -1/2 for all m) and the preceding analysis that |m| > 0 is needed for the angular integral to converge. Please correct the summary to match the derived condition.
  3. [Abstract and Sec. I] There are typos ('behavoiur', 'To address this, We develop...') and several equation references are imprecise (e.g., the discussion around Eq. (107) refers to 'equation (28)' rather than the displayed numbering). A careful editorial pass is needed.
  4. [Sec. V.D and Figure 1] The numerical experiments are not reproducible from the text: the Galerkin basis size, quadrature order, and the value of the regularization ε in Eq. (C10) used for the displayed results are not stated, and the claimed boundary-condition error (1.2%) is not defined precisely.

Circularity Check

2 steps flagged · score 6.0 of 10

The ℓ > −1/2 finite-energy criterion is read off from an assumed regularizing scaling (Eqs. 189–190) whose derivation is deferred to a missing appendix, while the numerical validation optimizes spectral weights under the same convergence constraint (Eq. 294).

  1. self definitional [Section IV.C.3–IV.E, Eqs. (189)–(192) and (195); derivation deferred to absent Appendix A]
    "Eϕ ∼ rℓ−1 sin|m|−1 θ · F(ℓ, m; ξ), (189) ... F (ℓ, m; ξ) = 1 + O(ξ2) for ξ → 0. (190) ... Wϕ ∼ Φ0 Z ϵ 0 r2ℓ dr Z π 0 sin2|m|−1 θ dθ. (192) ... the radial integral converges for ℓ >−1/2"

    The headline criterion ℓ > −1/2 is obtained by inserting the assumed scaling Eφ ∼ r^{ℓ−1} sin^{|m|−1}θ with F = 1 + O(ξ²) into the energy integral; the integrand r^{2ℓ} and threshold 2ℓ > −1 are algebraic consequences of the chosen exponent, not outputs of Maxwell's equations. The existence and asymptotics of F are the entire load, but they are only asserted, with the derivation deferred to an Appendix A that is not present. The paper's own Eq. (210) still gives Eφ ∼ m/(r^3 sin^2 θ) for the same modes, which is not energy-integrable for m ≠ 0; Eq. (189) is the mechanism that converts divergence into convergence, so the proof reduces to assuming the regularizing factor.

  2. fitted input called prediction [Section V.C.3.d, Eqs. (292)–(294); Figure 1(d)]
    "This establishes the constraint p > 3/2 − ℓ for energy convergence. (292) ... subject to p > 3/2 − ℓmin, q > 0, β > 0, (294) ... For the test case with ℓmin = 0.1, the resulting parameters were A = 0.832, p = 1.74, q = 0.92, and β = 0.437."

    The spectral weight parameters are optimized under the constraint p > 3/2 − ℓ_min, which is exactly the energy-convergence condition derived in Eq. (292) from the assumed scaling. The finite-energy fields displayed in Figure 1(d) are then generated with these optimized weights, so the numerical 'confirmation' of singular-mode energy convergence enforces the target condition by construction rather than testing it against an independent prediction.

full rationale

Score 6: the circularity is partial but central. The main theorem, finite electromagnetic energy for ℓ > −1/2, is carried by Eqs. (189)–(190), where the regularizing function F is asserted to make Eφ scale as r^{ℓ−1} sin^{|m|−1}θ with F = 1 + O(ξ²). Inserting this assumed scaling yields 2ℓ > −1 immediately; the paper's own unregularized asymptotics in Eq. (210) give a different, non-integrable Eφ scaling for m ≠ 0 and are never reconciled. The derivation of F is referenced to Appendix A, which is absent from the manuscript, so the presented proof has no independent content beyond the assumed exponent. The numerical validation in Figure 1(d) is also partly confirmation by construction: the parameters A, p, q, β are optimized subject to p > 3/2 − ℓ_min (Eq. 294), which is the paper's own convergence bound. I found no load-bearing self-citation chain: refs. [8,9] are invoked for the cylindrical analogue and for consistency, but the spherical criterion is not reduced to them. The reduction is through the asserted regularizing function and the fitted spectral weights, so the score is 6 rather than higher. This is a correctness and completeness gap as much as a circularity: a proof of existence of F with the stated properties, if present, could break the circular step.

Assumptions & free parameters 4 free parameters · 4 assumptions · 4 invented entities

The framework is held together by an assumed continuous-index basis and an asserted regularizing function. The only explicit free parameters are the fitted spectral weight coefficients, and the fit is constrained by the energy condition the paper claims to validate.

free parameters (4)
  • A (spectral weight amplitude) = 0.832
    Fitted by SQP to minimize boundary-condition error in Eq. (293); used to construct the finite-energy field shown in Figure 1(d).
  • p (spectral weight exponent) = 1.74
    Fitted subject to constraint p > 3/2 - ℓ_min (Eq. 294), which is the energy-convergence condition; this makes the finite-energy validation partly circular.
  • q (spectral weight exponent) = 0.92
    Fitted with p and β in the constrained optimization for the spectral weight function.
  • β (spectral weight exponent) = 0.437
    Fitted in the same optimization; controls decay of the spectral weight.
assumptions (4)
  • standard math Spectral theorem for self-adjoint operators: a vector Laplacian or curl-curl operator on a suitable function space has discrete modes plus a continuum of improper modes.
    Invoked in Section III.G.1 by analogy to justify the continuous spectral decomposition; no proof is given for the curl-curl operator with continuous angular indices.
  • domain assumption Analytic continuation of spherical harmonics to non-integer degree and order on appropriate Riemann surfaces is legitimate and yields admissible solutions.
    Section II relaxes 2π azimuthal periodicity and regularity at poles to allow real ℓ,m; the Riemann surface and boundary conditions are not rigorously defined.
  • ad hoc to paper The biorthogonal system Ψ_ℓm with weight w(ℓ,m) satisfies the resolution of identity in Eqs. (70) and (72).
    The paper says this is a logical extension and appeals to analogy with hyperbolic spectral theory; completeness is not proved and the weight has poles at integer ℓ-|m|.
  • ad hoc to paper The regularizing function F(ℓ,m;ξ) exists with F=1+O(ξ^2), making all angular-component energy integrals converge.
    Introduced in Eq. (189)-(190) specifically to remove the apparent 1/r divergence; no derivation is given.
invented entities (4)
  • Generalized continuous-index vector spherical functions Φ_ℓm
    purpose: Basis for spectral decomposition of vector fields over real ℓ,m
    No external falsifiable handle; orthogonality and completeness are assumed, not proved.
  • Spectral weight function w(ℓ,m) = π Γ(ℓ+|m|+1)/(sin(π(ℓ-|m|))Γ(ℓ-|m|+1))
    purpose: Normalizes the continuous-index basis and enforces the resolution of identity
    Its validity as a measure is assumed; it diverges at integer ℓ-|m| and no limit procedure is specified.
  • Angular operator L_ang with potential V(θ,φ)
    purpose: Encodes coupling between vector components for non-separable modes
    The eigenvalue problem is stated but the operator's domain, self-adjointness, and potential are not specified.
  • Regularizing function F(ℓ,m;ξ)
    purpose: Suppresses divergent Eθ and Eφ behavior near the origin so energy converges
    Its leading behavior is chosen to make the energy integral converge; no independent evidence or derivation is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Full Vectorial Maxwell Equations with Continuous Angular Indices." pith.science (2026). https://pith.science/paper/U6XJYVVT

@misc{pith2026250802675,
  author       = {Pith},
  title        = {Pith review of: Full Vectorial Maxwell Equations with Continuous Angular Indices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/U6XJYVVT}},
  note         = {Machine review of arXiv:2508.02675}
}
abstract

This article presents a mathematical framework for solving Maxwell's equations in cylindrical and spherical geometries with continuous angular indices. We extend beyond standard discrete harmonic decomposition to a continuous spectral representation using generalized spectral integrals, capturing electromagnetic solutions that exhibit singular behavoiur yet yield finite-energy fields at the geometric center. For continuous angular indices $\ell, m \in \mathbb{R}$, we study existence and uniqueness of solutions in weighted Sobolev spaces $H^s_{\alpha(\ell,m)}(\Omega)$ following the framework established in ~\cite{adams2003, reed1975}, prove finite energy for $\ell > -\frac{1}{2}$, and construct explicit spectral kernels via biorthogonal function systems. The framework encompasses both separable cylindrical modes with continuous azimuthal index $\nu \in (0,1)$ and non-separable spherical modes where field components couple through vectorial curl operations. We present asymptotic analysis of singular field behavior, investigate convergence rates for spectral approximations, and validate the theoretical framework through Galerkin projection methods and numerical spectral integration.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages

  1. [1]

    Characterises the singularity structure analytically

    work page

  2. [2]

    Establishes convergence and regularity in weighted Sobolev spaces

    work page

  3. [3]

    Full Vectorial Maxwell Equations with Continuous Angular Indices

    Avoids separable ansatz by solving the full coupled angular PDE system. The analysis of these singular structures in cylindrical geometry, where the variables remain separable and continuous azimuthal modes naturally arise from broken 2 π periodicity, is presented in 8,9. In this simpler case, we showed that for continuous azimuthal index ν ∈ (0, 1), the ...

    work page Pith review arXiv 2025

  4. [4]

    In classical spherical harmonic analysis, projection onto the ( ℓ, m) mode is achieved by an integral of the field against Y ∗ ℓm(θ, ϕ) over the solid angle

    Derivations and Projection Operators We introduce projection operators that project the fields onto these generalized spherical harmonics with continuous indices. In classical spherical harmonic analysis, projection onto the ( ℓ, m) mode is achieved by an integral of the field against Y ∗ ℓm(θ, ϕ) over the solid angle. For non-integer and continuous ℓ, we...

    work page

  5. [5]

    In spherical coordinates, even in standard integer-ℓ analysis, one often expands the vector field in vector spherical harmonic basis functions (which mix components Er, Eθ, Eϕ)

    Coupling of Vector Components A significant mathematical challenge in full Maxwell equations (in contrast to scalar wave equations) is the coupling introduced by the vector nature of the fields. In spherical coordinates, even in standard integer-ℓ analysis, one often expands the vector field in vector spherical harmonic basis functions (which mix componen...

    work page

  6. [6]

    Near r → 0, the solutions are singular

    Asymptotic Behavior One noteworthy result is our analysis of the asymptotic behavior of these continuous modes. Near r → 0, the solutions are singular. We find solutions behaving like Er ∼ rα(ℓ,m) with some exponent α that may be negative (indicating a singularity at the origin). Crucially, we show these singular solutions are integrable in the energy sen...

    work page

  7. [7]

    Instead, it is given by ∇2A = ∇(∇ ·A) − ∇ ×(∇ ×A)

    Vector Laplacian in Curvilinear Coordinates In spherical coordinates, the Laplacian of a vector field is not the component-wise application of the scalar Laplacian. Instead, it is given by ∇2A = ∇(∇ ·A) − ∇ ×(∇ ×A). (6) This identity ensures compatibility with Maxwell’s equations in vacuum. Importantly, in spherical coor- dinates ( r, θ, φ), the coordinat...

    work page

  8. [8]

    These satisfy orthogonality relations over the sphere and serve as the basis for studying wave equations and gauge theories on curved manifolds

    Expansion in Vector Spherical Harmonics To analyse angular structure systematically, we expand A in the vector spherical harmonics (VSH), which form a complete orthogonal basis for square-integrable vector fields on the sphere A(r, θ, ϕ) = ∞X ℓ=0 ℓX m=−ℓ h a(r) ℓm(r) Y (r) ℓm (θ, ϕ) + a(1) ℓm(r) Y (1) ℓm (θ, ϕ) + a(2) ℓm(r) Y (2) ℓm (θ, ϕ) i , (8) where t...

    work page

Show all 46 references
  1. [9]

    In particular the term ∇(∇ ·A) produces both radial and angular derivatives acting on Aθ, Aϕ

    Coupling under the Laplacian When the vector Laplacian acts on A, the radial and angular parts are no longer separable in general. In particular the term ∇(∇ ·A) produces both radial and angular derivatives acting on Aθ, Aϕ. Similarly, the term ∇ ×(∇ ×A) introduces mixing betw...

  2. [10]

    While scalar harmonics Y m ℓ diagonalise the Laplacian for scalar fields, the vector Laplacian does not commute with projection onto the radial or angular directions

    Limitations of Standard Angular Separation This coupling presents an obstacle to conventional separation of variables. While scalar harmonics Y m ℓ diagonalise the Laplacian for scalar fields, the vector Laplacian does not commute with projection onto the radial or angular dir...

  3. [11]

    Vector Expansion of Electromagnetic Fields We consider time-harmonic fields with dependence e−iωt in source-free, homogeneous, isotropic media. We expand the electric field in a basis of generalized vector spherical functions ⃗E(⃗ r) = X ℓ,m E(ℓ,m) r (r)Y m ℓ (θ, φ)ˆr + E(ℓ,m)...

  4. [12]

    Spectral Theory on Non-Compact Domains In spectral theory, the decomposition of functions on domains with broken symmetries generally involves both discrete and continuous components. For a vector function f defined on such a domain, the spectral representation takes the form ...

  5. [13]

    Vector-Valued generalized integral For vector fields satisfying Maxwell’s equations, we require a vector-valued extension of the generalized integral formalism. The electric field is represented as a spectral integral E(r) = Z C� Z C� a(ℓ, m)E(ℓ, m; r), dℓ, dm, (78) where E(ℓ,...

  6. [14]

    The coefficients α(ℓ, m) and β(ℓ, m) are determined by boundary conditions and regularity requirements

    Radial Functions and Green’s Function Approach For the radial component, we have Er(ℓ, m; r) = α(ℓ, m)h(1)ℓ(kr) + β(ℓ, m)h(2)ℓ(kr), (80) where h(1)ℓ and h(2)ℓ are the spherical Hankel functions of the first and second kind, analytically continued to non-integer order ℓ. The co...

  7. [15]

    Construction of Complete Solutions The full spectral representation of the electromagnetic field solutions requires a systematic approach that accounts for the coupled nature of the vector components. The radial component serves as the fundamental building block, with the gene...

  8. [16]

    generalized integral Interpretation We can express the spectral representation in terms of the parameters rather than ℓ. Using the relationship in equation 160, we can rewrite the spectral integral as E(r) = Z C� Z C� ˜a(s, m)E(s, m; r) ds dm, (96) where ˜a(s, m) = a(ℓ(s), m) ...

  9. [17]

    Green’s Function Method for Angular Components We begin by treating the equations for the angular components as inhomogeneous linear ODEs. For the θ-component, we have LθEθ(ℓ, m; r) = Sθ(ℓ, m; r), (97) where the differential operator Lθ is Lθ = d2 dr2 + 2 r d dr + ℓ(ℓ + 1) − 1...

  10. [18]

    (104) These solutions satisfy the appropriate boundary conditions

    Construction of Homogeneous Solutions For the operatorLθ, the appropriate homogeneous solutions are the spherical Bessel functions and spherical Neumann functions, analytically continued to non-integer order ν = p ℓ(ℓ + 1) − 1 uθ(r) = jν(kr), v θ(r) = yν(kr). (104) These solut...

  11. [19]

    For physical fields, we typically set ε → 0 and R → ∞, with additional conditions to ensure the solution remains well-behaved at these limits

    Integral Representation for � � Using this Green’s function, the solution to equation (28) can be written as Eθ(ℓ, m; r) = C1(ℓ, m)jν(kr) + C2(ℓ, m)yν(kr) + Z R ε Gθ(r, r′)Sθ(ℓ, m; r′)dr′, (107) where C1(ℓ, m) and C2(ℓ, m) are constants determined by boundary conditions, and t...

  12. [20]

    Integral Representation for � � Similarly, for the φ-component, we have LφEφ(ℓ, m; r) = Sφ(ℓ, m; r), (109) where Lφ = Lθ and Sφ(ℓ, m; r) = 1 r2 Er(ℓ, m; r)Bφ(ℓ, m) + 1 r2 d dr [r2Eθ(ℓ, m; r)]Cθφ(ℓ, m)∗. (110) The Green’s function solution is Eφ(ℓ, m; r) = D1(ℓ, m)jν(kr) + D2(ℓ...

  13. [21]

    Green’s Function Analysis and Contour Selection The integral representations derived above rely on the tensor Green’s function for Maxwell’s equations, which satisfies ∇ × ∇ ×G(r, r′) − k2G(r, r′) = Iδ (r − r′), (113) Full Vectorial Maxwell Equations with Continuous Angular In...

  14. [22]

    Coupled Integral System The integral representations derived in the previous sections form a coupled system of integral equations, since Eθ depends on Eφ and vice versa through the source terms. This coupling is fundamental to the non- separable nature of Maxwell’s equations w...

  15. [23]

    Spectral Representation of Angular Components For the complete field representation, we express Er in its spectral form Er(r, θ, φ) = Z C� Z C� a(ℓ, m)Er(ℓ, m; r)Y ℓ m(θ, φ)dℓdm, (127) Full Vectorial Maxwell Equations with Continuous Angular Indices 19 where Er(ℓ, m; r) = α(ℓ,...

  16. [24]

    Explicit Form of the Spectral Kernels To derive explicit expressions for the spectral kernels, we solve the coupled integral equations to first order in the coupling parameters. This provides an analytic approximation that captures the essential behavior Kθ(ℓ, m; r) = C1(ℓ, m)...

  17. [25]

    Rigorous treatment of both regular and singular field configurations

  18. [26]

    When an excitation frequency ω approaches a resonant frequency corresponding to a singular mode, the response amplitude scales as A(ω) ∼ 1 |ω − ωn|

    Natural incorporation of boundary conditions through the spectral weight function Full Vectorial Maxwell Equations with Continuous Angular Indices 21 The strength of resonant excitation depends on how closely an external driving matches the pole structure of a(ℓ, m). When an e...

  19. [27]

    (193) For the special case m = 0, we need a more careful analysis

    Combined Energy Convergence Conditions Combining the constraints from all components, we find that the total electromagnetic energy in the spherical case is finite if ℓ >− 1 2 and m ̸= 0. (193) For the special case m = 0, we need a more careful analysis. When m = 0, the ϕ comp...

  20. [28]

    For our field solutions with continuous indices, we have ⃗E ∈ H s w(Ω) for s < ℓ+ 1 2

    Function Space Characterization In terms of function spaces, our solutions belong to the weighted Sobolev space H s w(Ω) = f |f |H�� < ∞ , (196) where the weighted norm is defined as |f |2 H�� = Z Ω (1 + |ξ|2)s| ˆf (ξ)|2w(ξ) dξ, (197) Full Vectorial Maxwell Equations with Cont...

  21. [29]

    The field components have the following asymptotic behavior near r = 0 Er ∼ rℓ−1, E θ ∼ ℓ r , E φ ∼ m r3 sin2 θ , (210)

  22. [30]

    Energy convergence requires ℓ >− 1 2 and m = 0, (211)

  23. [31]

    In terms of eigenvalues, the condition is λℓm > − 3 4 , (212)

  24. [32]

    The singular field is correctly described by Esing(r, θ, φ) = Z 1 0 Z C� a(ℓ, m)rα(ℓ,m)Φℓm(θ, φ) dℓ dm, (213) with α(ℓ, m) = 1 2 √1 + 4λℓm − 1

  25. [33]

    The singularity occurs only in the radial coordinate r → 0, while the angular dependence remains regular. F. Non-Separable Electromagnetic Modes in Spherical Cavities with Continuous Angular Indices To provide a complete and rigorous treatment analogous to our cylindrical anal...

  26. [34]

    The non-separable nature of these equa- tions for continuous indices necessitates a more general expansion framework

    Coupled System of Equations To properly account for the coupling between angular indices, we write Maxwell’s curl equations explicitly in spherical coordinates 1 r sin θ ∂ ∂θ (Eϕ sin θ) − ∂Eθ ∂ϕ = iωµ0Hr, (224) 1 r 1 sin θ ∂Er ∂ϕ − ∂ ∂r (rEϕ) = iωµ0Hθ, (225) 1 r ∂ ∂r (rEθ) − ∂...

  27. [35]

    We discretize this operator by computing its matrix elements in our chosen basis [Lang]ij = ⟨ ⃗ψi, Lang ⃗ψj⟩

    Matrix Formulation of the Eigenvalue Problem The angular operator Lang derived from Maxwell’s equations takes the form Lang[⃗Φ] = − 1 sin θ ∂ ∂θ sin θ ∂⃗Φ ∂θ ! − 1 sin2 θ ∂2⃗Φ ∂ϕ2 + V (θ, ϕ)⃗Φ, (264) where V (θ, ϕ) is a tensorial potential term that encodes the coupling betwee...

  28. [36]

    Numerical Implementation and Convergence The matrix eigenvalue problem is solved using the implicitly restarted Arnoldi method for sparse matrices, which efficiently computes selected eigenvalues and eigenvectors. Our implementation proceeds as follows �� ������ Range of ℓ ∈ [...

  29. [37]

    �������������� We evaluate the objective function at Nθ × Nφ = 32 × 64 points on the sphere

  30. [38]

    ������� ����� We set initial values A0 = 1.0, p0 = 2.0, q0 = 1.0, β0 = 0.5

  31. [39]

    ���������� �������� We implement the constraint p >3 2 − ℓmin using a logarithmic barrier function B(p) = −µ log(p − ( 3 2 − ℓmin)), (295) where µ is a small positive parameter (typically µ = 10−3) that decreases during optimization

  32. [40]

    �������� ����������� We compute gradients of the objective function using automatic differentia- tion to enhance accuracy

  33. [41]

    Full Vectorial Maxwell Equations with Continuous Angular Indices 39

    ���� ������ We employ a backtracking line search with Armijo conditions f (x + αd) ≤ f (x) + c1α∇f (x)T d, (296) with typical values c1 = 10−4 and initial step size α = 1.0. Full Vectorial Maxwell Equations with Continuous Angular Indices 39

  34. [42]

    This optimization approach consistently produces spectral weight functions that accurately satisfy boundary conditions while maintaining physical constraints

    ����������� �������� We terminate the optimization when ∥∇L∥ < ϵgrad or |fk − fk−1| |fk| < ϵrel, (297) where L is the Lagrangian, fk is the objective function value at iteration k, and typical values are ϵgrad = 10−5 and ϵrel = 10−7. This optimization approach consistently pro...

  35. [43]

    �������� ��������� ������� We dynamically scale parameters to balance their numerical influence during optimization

  36. [44]

    ���������������� �������� We start with a coarse grid and progressively refine it as optimization proceeds

  37. [45]

    ���������������The objective function evaluation is parallelized across angular points for efficiency

  38. [46]

    �������������� ���������� The parameter λ is adjusted adaptively based on the current divergence magnitude. This detailed implementation ensures robust determination of the spectral weight function, balancing accu- racy at boundaries with physical constraints including energy ...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.