Vanishing Clutter: Fast and Accurate Shape Imaging via Active Cloaking

Abstract: Imaging a target sample in near-field scanning optical microscopy (NSOM) is fundamentally limited by measurement artifacts and data contamination from multiple probe-sample scattering. The probe, essential for subwavelength resolution, inherently perturbs the local field, degrading the signal-to-noise ratio and rendering the inverse problem for quantitative shape reconstruction highly ill-posed. We first establish the well-posedness of the corresponding forward model, providing a rigorous foundation for subsequent imaging. We then repurpose active cloaking--conventionally the antagonist of imaging--as an enabling mechanism to eliminate probe-induced interference. Rather than directly reconstructing the sample from corrupted data, we actively cloak the probe by formulating an optimal control problem and prove the existence and stability of its minimizers. Leveraging the theory of localized anomalous resonance in layered plasmonic structures, we derive an exact closed-form minimizer, thereby circumventing the computationally prohibitive iterative solution of the optimal control problem. The resulting cloaking-driven interference removal yields a virtually probe-free measurement environment, enabling fast, artifact-free shape reconstruction. Extensive numerical experiments demonstrate accurate shape reconstruction and dramatic acceleration--often by orders of magnitude--over conventional iterative methods applied directly to probe--contaminated data without cloaking--based preprocessing, validating the robustness and transformative potential of the proposed approach for high-fidelity subwavelength imaging.

• The paper demonstrates that active cloaking can suppress probe and interaction scattering through anomalous localized resonance, making measurements converge toward those from an isolated sample as loss approaches zero.
• The method combines an optimal-control formulation with a reduced-domain Levenberg–Marquardt reconstruction algorithm, which matches full-domain accuracy while achieving speedups from about 4× for single-layer probes to nearly 90× for complex multi-probe systems.
• The rigorous invisibility result applies to separated samples and concentric circular shells with specific plasmonic parameters, leaving eccentric geometries, three-dimensional settings, finite wavelengths, and theoretical noise robustness open for future work.

Problem setting and motivation

Near-field scanning optical microscopy (NSOM) achieves subwavelength resolution by bringing a probe into the extreme near field of a sample, but this proximity is precisely what corrupts the data: probe–sample multiple scattering generates standing waves, scattered surface plasmon polaritons, and spurious "clutter" that degrades the signal-to-noise ratio and makes quantitative shape reconstruction severely ill-posed. The paper by Xiao and Zheng (2608.13839) addresses this by inverting the usual role of cloaking. Rather than hiding a target from an observer, the authors cloak the probe itself, so that its scattering signature is suppressed before reconstruction is attempted. The central claim is that active cloaking can serve as a physics-driven preprocessor that yields measurements equivalent to those of an ideal, invisible sensor.

The mathematical model consists of three cylindrical (effectively two-dimensional) domains of class $C^{1,\eta}$: a penetrable sample $D_1$ with permittivity $\varepsilon_1>0$, a probe core $D_2$ with $\varepsilon_2>0$, and a cloaking shell $D_3\setminus\overline{D_2}$ filled with a Drude-like metamaterial of permittivity $\varepsilon_3=-s+i\delta$, where $s>0$ is the control parameter and $\delta>0$ models plasmonic loss. Under quasi-static illumination $u^i(x)=x\cdot d$, the total potential satisfies a transmission problem with continuity of potential and weighted normal flux across all interfaces, together with the decay condition $D_1$0. The inverse problem is to recover $D_1$1 from total-field measurements on a circle $D_1$2, assuming the probe geometry and parameters are known.

Well-posedness of the forward problem

The authors first establish uniqueness via an energy method: taking the imaginary part of a Green's identity over $D_1$3 forces $D_1$4 in the lossy shell (since only $D_1$5 has nonzero imaginary part), and the real part then annihilates all gradients, yielding $D_1$6 for homogeneous data. This argument depends essentially on $D_1$7; the lossless case is not covered.

Existence proceeds through layer potentials. The field in the presence of $D_1$8 alone admits the explicit representation

$D_1$9

with $\varepsilon_1>0$0, and the full solution is expanded using Green's functions modified by $\varepsilon_1>0$1. The densities on $\varepsilon_1>0$2 and $\varepsilon_1>0$3 satisfy a compactly perturbed diagonal system $\varepsilon_1>0$4, where invertibility of $\varepsilon_1>0$5 follows from $\varepsilon_1>0$6 under the assumed permittivity signs; Fredholm theory then gives existence. A Lipschitz stability result shows that the density pair $\varepsilon_1>0$7 depends on $\varepsilon_1>0$8 with constant $\varepsilon_1>0$9 uniformly on compact subsets $D_2$0, which is the key regularity input for everything that follows.

Optimal control formulation of probe cloaking

The probe's contribution to the measured field is decomposed as $D_2$1, where $D_2$2 collects the shell's own scattering plus all sample–probe interaction terms. Cloaking is formulated as minimizing

$D_2$3

over $D_2$4. Using the stability theorem, the authors prove existence of minimizers, Lipschitz stability of $D_2$5 in $D_2$6, and convergence of approximate minimizers computed from perturbed fields. These results give the optimal-control route to cloaking a rigorous footing, but they do not identify the minimizer explicitly, nor do they show that $D_2$7 as $D_2$8. That gap is closed by the resonance analysis below.

Closed-form invisibility via anomalous localized resonance

For concentric disks $D_2$9, $\varepsilon_2>0$0 with $\varepsilon_2>0$1, permittivities $\varepsilon_2>0$2, $\varepsilon_2>0$3, one obtains $\varepsilon_2>0$4, reducing the density system to $\varepsilon_2>0$5. The critical radius is $\varepsilon_2>0$6, matching the classical CALR threshold. Two structural estimates drive the analysis:

• Sample-probe coupling is weak at distance: if the sample center satisfies $\varepsilon_2>0$7 and is sufficiently far, the interaction operator $\varepsilon_2>0$8 satisfies $\varepsilon_2>0$9 with $D_3\setminus\overline{D_2}$0, so it is a small compact perturbation.
• Boundedness off the critical region: for sources supported outside $D_3\setminus\overline{D_2}$1 whose Newtonian potential is harmonic there, the Fourier-mode solution of $D_3\setminus\overline{D_2}$2 remains uniformly bounded as $D_3\setminus\overline{D_2}$3, because $D_3\setminus\overline{D_2}$4 for small $D_3\setminus\overline{D_2}$5.

Combining these via an asymptotic expansion in powers of $D_3\setminus\overline{D_2}$6, the authors prove their main analytical result: for $D_3\setminus\overline{D_2}$7 and $D_3\setminus\overline{D_2}$8,

$D_3\setminus\overline{D_2}$9

i.e., the entire probe-plus-interaction contribution vanishes outside the critical annulus, leaving $\varepsilon_3=-s+i\delta$0 — the field of the isolated sample. This connects the optimal-control formulation to CALR-based cloaking: the resonant choice $\varepsilon_3=-s+i\delta$1 is effectively the minimizer of $\varepsilon_3=-s+i\delta$2, obtained in closed form rather than iteratively. The implication is immediate: once the probe is cloaked, inversion may proceed against the much simpler single-inclusion forward model.

Two caveats attach to this result. It requires the concentric-disk geometry, the specific permittivity balance $\varepsilon_3=-s+i\delta$3, $\varepsilon_3=-s+i\delta$4, and sufficient probe–sample separation; and the proof treats the coupled term formally, relying on uniform boundedness of $\varepsilon_3=-s+i\delta$5 rather than an explicit expression for the coupled Fourier coefficients.

Reconstruction algorithm and numerical performance

Reconstruction parameterizes star-shaped boundaries via truncated Fourier series $\varepsilon_3=-s+i\delta$6 and applies Levenberg–Marquardt with $\varepsilon_3=-s+i\delta$7-seminorm regularization, solved by conjugate gradient. Two variants are compared: the full-domain method, which solves the complete probe–sample system ($\varepsilon_3=-s+i\delta$8 densities for an $\varepsilon_3=-s+i\delta$9-layer probe) per iteration, and the reduced-domain method, which solves only the isolated-sample problem (one density), valid when the probe is cloaked.

The numerical study covers small probes, large probes, multi-layer probes ($s>0$0), and multi-probe configurations, with synthetic data generated at higher discretization than used in inversion to avoid inverse crime. Key findings:

• Small-probe regime: with bare probes the reduced-domain method produces visibly distorted reconstructions; with cloaked probes it recovers disk targets accurately even with three simultaneous probes.
• Large-probe regime: both methods achieve comparable accuracy on circle, ellipse, and kite targets, but the reduced-domain method is dramatically faster. With $s>0$1 collocation points and $s>0$2, speedups grow with probe complexity:

Configuration	Speedup range
Single-layer probe ($s>0$3)	4.25–4.61
Three-layer probe ($s>0$4)	20.13–22.75
Five-layer probe ($s>0$5)	54.18–62.02
Multi-probe Case 1 (three single-layer)	23.78
Multi-probe Case 2	41.66
Multi-probe Case 3 (mixed layers up to 3)	89.85

The scaling is explained by linear-algebra cost: the reduced method solves an $s>0$6 system per forward solve versus $s>0$7 for the full method, and $s>0$8 stays nearly constant (about 0.6 s for one illumination, about 2.1 s for four) while $s>0$9 grows rapidly with layer count. The numerical experiments also confirm the theoretical predictions directly: $\delta>0$0 attains its minimum at $\delta>0$1 and blows up at the analytically predicted values $\delta>0$2 and $\delta>0$3, where $\delta>0$4 hits the NP eigenvalue structure; and the residual norm $\delta>0$5 tends to zero as $\delta>0$6 for $\delta>0$7 and $\delta>0$8 multi-layer probes and for all three multi-probe cases, indicating empirically that the cloaking property extends beyond the single-shell geometry covered by the theorem.

Limitations and open questions

Several restrictions should be noted plainly. The rigorous invisibility theorem holds only for concentric circular shells under the exact plasmonic balance $\delta>0$9, $u^i(x)=x\cdot d$0, and requires the sample to lie outside the critical radius $u^i(x)=x\cdot d$1 and sufficiently far from the probe; whether CALR-based cloaking persists for eccentric, non-circular, or slit-type NSOM probes is left open. The two-dimensional quasi-static model idealizes real three-dimensional probes, and the extension to finite-wavelength regimes is not addressed. The stability of the reconstruction under measurement noise is demonstrated only numerically, not theoretically, and the paper does not quantify how large $u^i(x)=x\cdot d$2 can be before the reduced-domain method loses accuracy. Finally, the equivalence between the optimal-control minimizer and the resonant value $u^i(x)=x\cdot d$3 is established analytically only in the symmetric configuration; a general proof connecting the two formulations remains open.

Conclusion

This work establishes a coherent framework in which active cloaking acts as a preprocessing step for inverse problems: the forward transmission problem is shown to be well-posed, probe cloaking is posed as an optimal control problem with provably existing and stable minimizers, and localized anomalous resonance supplies an explicit near-minimizer rendering the coated probe invisible as $u^i(x)=x\cdot d$4. The practical payoff is a reduced-domain reconstruction algorithm that matches full-domain accuracy while delivering speedups ranging from roughly 4× for single-layer probes to nearly 90× for mixed multi-probe configurations. The main open issues are the extension of the invisibility guarantee beyond concentric circular geometries and a theoretical noise-robustness analysis of the resulting inverse problem.