Secondary Wave-Source Response Field (SWRF)

• Secondary Wave-Source Response Field (SWRF) is a representation that encodes observable wave responses from latent, stochastic sources for inverse analysis.
• It unifies methods in stochastic inverse problems and differentiable radar imaging by enabling recovery of source statistics and soft shadow boundaries.
• SWRF facilitates both high-frequency asymptotic analysis and gradient-based inversion, ensuring robust reconstructions in various wave phenomena.

Searching arXiv for the cited SWRF-related papers and closely related terminology. Secondary Wave-Source Response Field (SWRF) denotes a representation in which the quantity of interest is not the primary source or object itself, but a field that encodes how that source or object gives rise to an observable wave response. In the stochastic inverse-source setting of wave equations, this viewpoint appears naturally when a rough random source is inferred indirectly through its radiated far-field pattern, which functions as a secondary response field carrying recoverable statistical structure (Li et al., 2021). In forward scatter radar (FSR), SWRF is introduced explicitly as a continuous, learnable scalar field $p(x',z')\in\mathbb{R}$ on the aperture plane $x'oz'$, interpreted as the propensity of each point $Q=(x',z')$ to act as a secondary wave source under incident illumination; this field replaces a non-differentiable hard shadow boundary by a soft physical surrogate for end-to-end inversion (Lei et al., 15 Aug 2025).

1. Conceptual definition and scope

The term SWRF spans two closely related usages. In one usage, it is a viewpoint for inverse problems in which the measured quantity is a propagated wave response rather than the original source. The paper on inverse source problems for stochastic wave equations fits this viewpoint “very naturally”: the original object of interest is a rough random source, but what is actually measured is a radiated wave field, and in the high-frequency and far-field regime that radiated field becomes a secondary response field encoding statistical information about the source (Li et al., 2021).

In the other usage, SWRF is a concrete latent field used inside an inversion framework. In DPI-SPR for FSR, SWRF is defined as a continuous, learnable scalar field $p(x',z')$ whose sigmoid activation

$\sigma(p)=\frac{1}{1+\exp(-p)}$

is treated as the probability that a point belongs to the shadow-casting region or secondary source region (Lei et al., 15 Aug 2025). High $p$ implies $\sigma(p)\approx 1$ inside the shadow boundary, while low $p$ implies $\sigma(p)\approx 0$ outside.

These two usages are methodologically aligned. In both cases, the observable wave quantity is treated as a structured response field from which hidden source information can be inferred. This suggests a unifying interpretation: SWRF is a representation that mediates between an unobserved source description and an observed wave phenomenon, with the secondary field chosen so that it is both physically meaningful and invertible in an appropriate asymptotic or optimization sense.

2. SWRF in stochastic wave inverse source problems

For stochastic acoustic, biharmonic, electromagnetic, and elastic wave equations, the direct source-to-response map has the form

$f \longmapsto u \longmapsto u^\infty,$

where $x'oz'$0 is the random source, $x'oz'$1 is the wave field, and $x'oz'$2 is the far-field pattern (Li et al., 2021). In this setting, the far-field pattern is the response field.

The source is modeled as a centered generalized microlocally isotropic Gaussian random field. “Centered” means

$x'oz'$3

For a complex-valued field $x'oz'$4, the covariance and relation operators are defined by

$x'oz'$5

Here $x'oz'$6 is the covariance operator and $x'oz'$7 is the relation operator, also called the pseudo-covariance operator (Li et al., 2021).

Microlocal isotropy of order $x'oz'$8 in a bounded domain $x'oz'$9 means that both operators are classical pseudo-differential operators of order $Q=(x',z')$0, with symbols

$Q=(x',z')$1

The amplitudes $Q=(x',z')$2 and $Q=(x',z')$3 are the principal symbol amplitudes of the covariance and relation operators, respectively (Li et al., 2021). Their Fourier representation and asymptotic decomposition are given by

$Q=(x',z')$4

$Q=(x',z')$5

$Q=(x',z')$6

and

$Q=(x',z')$7

Within this framework, SWRF is not an added latent variable but an interpretive role played by the far-field pattern. The paper’s central inverse statement is that the principal symbols of the covariance and relation operators can be uniquely determined by a single realization of the far-field pattern averaged over the frequency band with probability one (Li et al., 2021).

3. Direct source-to-response mappings and far-field structure

The direct problem is developed in a unified way for four wave models: the Helmholtz equation, the biharmonic wave equation, Maxwell’s system, and the Navier equation (Li et al., 2021). In each case, the random source generates a wave field via a linear Green-function representation, and the far-field is an asymptotic transform of that field.

For the stochastic Helmholtz equation,

$Q=(x',z')$8

with Sommerfeld radiation condition, the solution is

$Q=(x',z')$9

and the far-field pattern is

$p(x',z')$0

Accordingly, the far-field is essentially a Fourier transform of the random source restricted to the frequency sphere (Li et al., 2021).

For the biharmonic equation

$p(x',z')$1

the far-field pattern is

$p(x',z')$2

For Maxwell’s equations,

$p(x',z')$3

the electric far-field is written as

$p(x',z')$4

up to the vector/tensor structure coming from Maxwell’s Green tensor.

For the Navier equation,

$p(x',z')$5

the displacement has compressional and shear far-field parts: $p(x',z')$6 with

$p(x',z')$7

$p(x',z')$8

Here the far-field is matrix-valued, with polarization projectors selecting compressional and shear content (Li et al., 2021).

The direct problems are shown to be well posed in the sense of distributions, and the resulting fields possess almost sure Sobolev regularity sufficient for far-field asymptotics and frequency averaging. For the acoustic and elastic cases the paper assumes $p(x',z')$9; for the biharmonic case $\sigma(p)=\frac{1}{1+\exp(-p)}$0; and for Maxwell $\sigma(p)=\frac{1}{1+\exp(-p)}$1 under the divergence-free distributional condition $\sigma(p)=\frac{1}{1+\exp(-p)}$2 (Li et al., 2021).

4. Inverse recovery, ergodicity, and identifiability

The inverse problem asks what can be recovered about the random source from the far-field pattern. The answer is formulated through both ensemble correlations and frequency-band averages from a single realization (Li et al., 2021).

In the acoustic case, the high-frequency expectation identities are

$\sigma(p)=\frac{1}{1+\exp(-p)}$3

and

$\sigma(p)=\frac{1}{1+\exp(-p)}$4

Thus the leading behavior of far-field correlations is determined by the Fourier transforms of the principal symbol amplitudes of the covariance and relation operators (Li et al., 2021).

The practical replacement of ensemble expectation by a single-realization frequency average is achieved through ergodicity. Almost surely,

$\sigma(p)=\frac{1}{1+\exp(-p)}$5

and

$\sigma(p)=\frac{1}{1+\exp(-p)}$6

The proof uses covariance decay estimates such as

$\sigma(p)=\frac{1}{1+\exp(-p)}$7

with analogous bounds for the other models (Li et al., 2021).

The same basic structure is established for biharmonic, Maxwell, and elastic waves, with modified scaling exponents and tensorial or projected symbol structure. For Maxwell, for example,

$\sigma(p)=\frac{1}{1+\exp(-p)}$8

$\sigma(p)=\frac{1}{1+\exp(-p)}$9

The main identifiability result is that the principal symbol amplitudes $p$0 and $p$1, or their matrix-valued analogues, are uniquely determined by the far-field data. The paper states this as uniqueness on any open subset $p$2, with analyticity providing extension from open-set knowledge to uniqueness of the amplitudes themselves (Li et al., 2021). A plausible implication is that, within this asymptotic regime, SWRF is not merely descriptive but operational: the secondary field preserves exactly the source statistics that the inverse problem can recover.

5. SWRF as a learnable field in DPI-SPR for forward scatter radar

In DPI-SPR, SWRF is introduced to remove the non-differentiability of classical forward scattering models in FSR (Lei et al., 15 Aug 2025). The classical Fresnel-Kirchhoff diffraction model (FKDM) integrates over a hard shadow region $p$3 defined by the indicator

$p$4

Because this is a step function at the boundary, the forward model is not differentiable with respect to shape parameters.

SWRF replaces this hard geometric boundary with a soft differentiable surrogate. The classical FKDM form is written as

$p$5

The Differentiable Forward Scatter Model (DFSM) replaces the hard boundary by the SWRF probability map: $p$6 Here $p$7 is a computational domain enclosing the true shadow boundary, $p$8 are the learnable SWRF parameters, and $p$9 (Lei et al., 15 Aug 2025).

The geometry-dependent distances are given in the paper as

$\sigma(p)\approx 1$0

$\sigma(p)\approx 1$1

with explicit dependence on target center and baseline geometry. The paper states that DFSM is fully differentiable with respect to SWRF parameters, uses the exact spherical-wave physical model rather than simplifying paraxial approximations, and supports near-field, far-field, and large/small diffraction-angle cases (Lei et al., 15 Aug 2025).

SWRF is not the final image. Reconstruction proceeds in two stages: continuous optimization of $\sigma(p)\approx 1$2, then binary extraction from the converged probability map. The recovered shadow profile is

$\sigma(p)\approx 1$3

The method therefore implements a soft-to-hard transition: differentiable physical inversion first, thresholded profile extraction afterward (Lei et al., 15 Aug 2025).

6. Optimization, regularization, and empirical behavior in DPI-SPR

The SWRF parameters in DPI-SPR are learned by minimizing the composite loss

$\sigma(p)\approx 1$4

where $\sigma(p)\approx 1$5 is the data consistency loss, $\sigma(p)\approx 1$6 is binary entropy regularization, and $\sigma(p)\approx 1$7 is geometric regularization (Lei et al., 15 Aug 2025).

The data term uses a complex Huber loss combined with a logarithmic transform: $\sigma(p)\approx 1$8 with residual $\sigma(p)\approx 1$9. The gradient scaling is

$p$0

The factor $p$1 suppresses overly large gradients from outliers and is reported to improve convergence speed and robustness under impulsive noise (Lei et al., 15 Aug 2025).

The entropy regularizer is

$p$2

with

$p$3

It is introduced because sigmoid activation can saturate near $p$4 or $p$5, causing gradients to vanish too early; the regularizer is used to keep optimization active and avoid premature saturation around trivial solutions (Lei et al., 15 Aug 2025).

The geometric regularizer is

$p$6

with smoothness, connectivity, and compactness priors. These are described as suppressing jagged high-frequency boundary artifacts, preventing fragmented objects, and favoring spatially concentrated solutions consistent with target silhouettes (Lei et al., 15 Aug 2025).

The end-to-end gradient path is

$p$7

and the update step is

$p$8

The paper specifies gradient descent / Adam, initial learning rate $p$9, decay by $\sigma(p)\approx 0$0 every 300 iterations, up to 500 iterations for simple targets, and up to 1000 iterations for noisy complex targets (Lei et al., 15 Aug 2025).

A physically motivated initialization is also given: $\sigma(p)\approx 0$1 where $\sigma(p)\approx 0$2 is the distance from the point to the center of the aperture plane, $\sigma(p)\approx 0$3 controls peak intensity, and $\sigma(p)\approx 0$4 controls the soft radius or boundary sharpness (Lei et al., 15 Aug 2025).

The reported experiments evaluate triangle, square, trapezoid, and hexagon targets in both TZ and FF settings. The paper gives the following IoU values: triangle $\sigma(p)\approx 0$5 (TZ), $\sigma(p)\approx 0$6 (FF); square $\sigma(p)\approx 0$7 (TZ), $\sigma(p)\approx 0$8 (FF); trapezoid $\sigma(p)\approx 0$9 (TZ), $f \longmapsto u \longmapsto u^\infty,$0 (FF); and hexagon $f \longmapsto u \longmapsto u^\infty,$1 (TZ), $f \longmapsto u \longmapsto u^\infty,$2 (FF). It also reports low $f \longmapsto u \longmapsto u^\infty,$3 HD values, stable reconstruction for non-convex targets under 15 dB, 8 dB, and 8 dB with mixed Gaussian plus impulse noise, and effective reconstruction down to 8 dB (Lei et al., 15 Aug 2025).

The ablation study attributes distinct roles to the main components: replacing the logarithmic loss with plain $f \longmapsto u \longmapsto u^\infty,$4 caused catastrophic failure under noise; removing regularization led to unstable or implausible reconstruction; removing entropy regularization increased susceptibility to premature saturation and local minima; and removing connectivity, compactness, or smoothness degraded boundary coherence through fragmented contours and distorted shapes (Lei et al., 15 Aug 2025).

7. Relation between the two SWRF formulations

The two papers use SWRF at different levels of abstraction. In the stochastic wave-equation paper, the “response field” is the far-field pattern of the radiated wave, and the inverse problem recovers the principal symbols of source covariance and relation operators from asymptotic correlation identities and ergodic frequency averaging (Li et al., 2021). In DPI-SPR, SWRF is a learnable latent field inside the forward model itself, designed to replace a non-differentiable boundary representation and make inversion tractable by automatic differentiation (Lei et al., 15 Aug 2025).

Their common structure is nonetheless clear. In both formulations, a primary object that is difficult to access directly—a rough generalized random source in one case, an unknown shadow boundary in the other—is not inverted in raw form. Instead, inversion proceeds through a secondary field that is wave-mediated, physically structured, and amenable either to high-frequency asymptotics or to gradient-based optimization. This suggests that SWRF is best understood not as a single equation or model class, but as a representation principle for wave-based inverse problems.

A common misconception would be to treat SWRF as necessarily identical across domains. The available evidence does not support that. In (Li et al., 2021), SWRF is an interpretive viewpoint centered on the far-field pattern as the measured secondary response field. In (Lei et al., 15 Aug 2025), SWRF is explicitly defined as a continuous, learnable scalar field $f \longmapsto u \longmapsto u^\infty,$5. The shared terminology therefore indicates a family resemblance rather than a single canonical mathematical object.

Another plausible implication is that SWRF-based formulations are especially useful when the direct object representation is either distributional and statistically structured, or geometrically sharp and non-differentiable. In the first case, asymptotic response statistics recover microlocal information; in the second, soft activation and differentiable physics permit end-to-end inversion. Both usages preserve a strong connection between latent representation and wave propagation, which is the central technical feature of the term as currently evidenced in the arXiv literature (Li et al., 2021, Lei et al., 15 Aug 2025).