{"id":"8a701fc4-aa16-46be-a352-c2bc4edeb406","arxiv_id":"1908.05800","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The orthogonality sampling method is rigorously justified for electromagnetic inverse scattering of anisotropic media, with resolution, stability, and equivalence to the direct sampling method.","lead":"This paper proves a theoretical foundation for the orthogonality sampling method, a fast non-iterative imaging technique that locates anisotropic electromagnetic scatterers from far-field data at one frequency. It establishes resolution and stability estimates and shows the method is equivalent to the direct sampling method, with three-dimensional numerical examples.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 6's decay proof contains an incorrect Funk-Hecke computation; the displayed v_j formula is quantitatively wrong, so the O(1/dist^2) decay of IOSM is not established as written, though the rate appears salvageable.","rationale":"The reader's weakest_assumption points to the coercivity of T on Range(H) imported from Haddar [10]. That is a genuine condition, but it is an explicitly stated, standard hypothesis, not an internal defect. My stress-test pass found a more immediate, concrete internal problem: the Funk-Hecke computation in the proof of Theorem 6 is quantitatively wrong, so the decay statement of the central theorem is not proved as written. This is load-bearing because the separation between inside and outside values of the imaging functional relies on both the lower bound and the O(1/dist^2) decay; if the decay proof fails, the plotting criterion is not justified by the argument given. However, the correct computation preserves the O(1/|z|) decay of H phi_ys, so the theorem's conclusion is very likely true and only the proof needs correction. The unproved Theorem 10 is a minor completeness issue, and the numerical experiments are supportive but not independently reproducible. Overall, the paper's central claims appear sound in substance but need a corrected derivation before the proof is fully convincing; this is consistent with the reader's CONDITIONAL verdict, so I recommend no change to that verdict. The disagreement with the reader's weakest_assumption is partial: the coercivity is a limitation, but the incorrect v_j formula is the sharper and more actionable concern.","tokens_in":12070,"tokens_out":31891,"duration_ms":311106,"concrete_test":"Recompute the integral I(z)=∫_{S2}(d×p)×d e^{-ik d·z} ds(d) for z=(1,0,0), p=(1,0,0), k=1 by direct numerical quadrature and compare with the paper's v_j formula in the proof of Theorem 6. Also verify the corrected Funk-Hecke identity I(z)=4π[p(j0(kr)-j1(kr)/(kr)) - (p·z)z/r^2(j0(kr)-3j1(kr)/(kr))] and check that it still satisfies |I(z)|=O(1/|z|) as |z|→∞, which is the only property needed for the decay part of Theorem 6.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central imaging claim in Theorem 6 depends on the decay estimate IOSM(ys) = O(1/dist(ys,Omega)^2), which is proved by bounding ||H phi_ys||^2 via an explicit Funk-Hecke computation. The displayed formula for the integral I(z)=∫_{S2}(d×p)×d e^{-ik d·z} ds(d) is incorrect. Taking z=(1,0,0), p=(1,0,0), k=1, direct quadrature or the standard identity I(z)=4π[p(j0(kr)-j1(kr)/(kr)) - (p·z)z/r^2(j0(kr)-3j1(kr)/(kr))] gives I(z)=8π j1(1)e1 ≈ 7.57 e1. The paper's v_j formula gives v_1=0.9036, and since I=-4π v/k^2, this yields I≈-11.35 e1: wrong sign and magnitude. Thus the proof of the decay statement is not valid as written. The decay rate itself is plausibly correct (stationary phase gives H phi_ys = O(1/r) uniformly for x in Omega), so the theorem may be repairable, but the current paper does not supply a correct derivation. The lower-bound half of Theorem 6 also inherits the external coercivity condition of Theorem 4 (cited to Haddar [10]); that is an explicit assumption rather than an internal error, but it is a further reason the central claim is conditional on material hypotheses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper treats the electromagnetic inverse shape problem for penetrable anisotropic scatterers at fixed frequency, using multi-static far-field data. It defines an orthogonality sampling imaging functional IOSM(ys), proves (Theorem 6) that IOSM is bounded below by a positive constant for sampling points inside the scatterer and decays as O(dist(ys,Ω)^{-2}) outside, establishes a noise-stability estimate (Theorem 7), and shows an equivalence between IOSM and a direct sampling functional IDSM (Theorem 9), with a corresponding resolution statement (Theorem 10). The analytical machinery is the factorization F=H*TH of the far-field operator, the range characterization of H*, and a coercivity result for T taken from previous work; numerical examples for three-dimensional anisotropic scatterers are also presented.","tokens_in":12380,"tokens_out":14888,"duration_ms":138461,"significance":"If the proofs are repaired, the paper would give a rigorous resolution analysis for two simple, fast, and robust sampling methods in three-dimensional electromagnetic inverse scattering, including anisotropic media. The use of external benchmarks is transparent and not circular: the range characterization of the factorization method and the coercivity of the middle operator are cited from Kirsch [14] and Haddar [10], respectively, and no ad-hoc parameters are introduced. The stability estimate is a useful quantitative addition. However, the central exterior-decay proof in Theorem 6 contains a concrete computational error, and Theorem 10 is stated without a proof; these issues are load-bearing for the main claims, although they appear repairable.","major_comments":[{"comment":"The decay estimate IOSM(ys)=O(dist(ys,Ω)^{-2}) is not established by the computation given. In the proof, the identity curl_z curl_z(p e^{-ikd·z}) = -k^2 (d×p)×d e^{-ikd·z} has the wrong sign; the correct identity has +k^2 (d×p)×d e^{-ikd·z}. The subsequent evaluation of I(z)=∫_{S^2}(d×p)×d e^{-ikd·z} ds(d) is also quantitatively incorrect. For example, with z=(1,0,0), p=(1,0,0), and k=1, direct computation gives I(z)=8π j_1(1)e_1 ≈ 7.57 e_1, whereas the displayed v_j formula gives v_1=0.9036 and the relation I=-4π v/k^2 yields I≈-11.35 e_1, so both sign and magnitude are wrong. Because this computation is the only support for the exterior decay assertion, the proof of Theorem 6 is invalid as written. The decay rate is plausible and likely repairable, for instance from the exact identity I(z)=4π[p(j0(kr)-j1(kr)/(kr))-(p·z)z/r^2(j0(kr)-3j1(kr)/(kr))] or by a stationary-phase argument, but it must be corrected.","section":"§4, proof of Theorem 6"},{"comment":"Theorem 10 is stated without a proof. The sentence \"Following the method for proving Theorem 6 we can obtain similar resolution analysis\" is not a proof, and the decay half would inherit the incorrect Funk-Hecke computation identified above. The authors should either supply the proof or state explicitly why the DSM functional inherits the results, since the theorem is one of the paper's main claims.","section":"§4, Theorem 10"}],"minor_comments":[{"comment":"In the Cauchy-Schwarz step, |⟨Fφ_ys,φ_ys⟩|^2 ≤ |S^2|^2 ‖p·Fφ_ys‖^2 should be |⟨Fφ_ys,φ_ys⟩|^2 ≤ |S^2| ‖p·Fφ_ys‖^2; the extra factor of |S^2| changes the constants in Theorem 6 and Theorem 9, although it does not affect the positivity or the decay rate.","section":"§4, proof of Theorem 6"},{"comment":"The stability estimate contains the notation |S|^2, which should presumably be |S^2|, and the displayed constant appears to be too large by a factor of |S^2|; checking with ‖φ_ys‖^2=(8π/3)|p|^2 gives a smaller prefactor. The qualitative stability statement is unaffected.","section":"§4, Theorem 7"},{"comment":"The indexing in the Funk-Hecke formula is nonstandard: one usually writes Y_ℓ^m with m=−ℓ,...,ℓ, rather than the displayed \"m∈N∪{0} and ℓ=−m,...,m\". This should be corrected to avoid confusion.","section":"§4, Funk-Hecke formula"},{"comment":"The entries \"1 .5\" and \"1 .2\" in the matrix A presumably mean 1.5 and 1.2; the spacing makes the values ambiguous.","section":"§5, equation (13)"},{"comment":"There is a typo in the sentence \"This resolution analysis implies the for any sampling points ys∈Ω...\" — \"the\" should be removed.","section":"§4, paragraph after Theorem 6"},{"comment":"The hypotheses of Theorem 4 (coercivity of T on Range(H)) and the transmission-eigenvalue assumption are not repeated in the statements of Theorems 6 and 9; since these conditions are load-bearing for the positivity of the imaging functionals inside Ω, they should be stated explicitly or referenced unambiguously in each theorem.","section":"§4, Theorems 6 and 9"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has a sound overall architecture and the lower-bound and equivalence arguments are mostly correct modulo constant factors, but the invalid Funk-Hecke computation in Theorem 6 touches the main decay claim, and Theorem 10 lacks a proof. These are local and likely repairable, so I recommend major revision rather than rejection. I would also encourage the authors to have the numerical section probe at least one far-field decay prediction, since the current experiments are purely qualitative."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nMy take: this paper extends the orthogonality sampling method (OSM) and direct sampling method (DSM) to electromagnetic scattering from anisotropic media with arbitrary-shaped scatterers. That is a real step beyond prior work on Helmholtz (Liu) and small scatterers (Nguyen). The claimed resolution, stability, and equivalence results are new relative to that work, and the paper is clearly written.\n\nThe soft spots are real. The proof of the decay estimate in Theorem 6 contains an incorrect Funk-Hecke computation. The displayed formula for v_j is quantitatively wrong: for a concrete case (z=(1,0,0), p=(1,0,0), k=1) it gives about -11.4 instead of +7.6, with the wrong sign and magnitude. The standard vector identity gives the correct O(1/r) behavior, so the claimed O(1/dist^2) decay is likely true and probably repairable, but the paper as written does not prove it. This is the core of the resolution analysis, so it is a load-bearing flaw, not a minor typo.\n\nThere is also a minor typo in the Cauchy-Schwarz step in Theorem 6: the factor should be |S^2|, not |S^2|^2. More importantly, the proof of Theorem 10 (the analogous resolution result for DSM) is omitted entirely; the equivalence theorem alone does not give the claimed O(1/dist^2) decay, so a real proof is needed. The lower-bound half of Theorem 6 rests on a coercivity assumption on the middle operator T, imported from Haddar. That is an explicit material condition, which is fine, but it means the theoretical guarantee is conditional.\n\nOn the positive side, the paper builds on the factorization method and does not fit anything. The stability estimate is straightforward and correct. The numerical experiments are synthetic and lack code/data, but they show the expected behavior for several shapes.\n\nThis is a genuine within-subfield contribution. I would send it to peer review, but only with the expectation of major revision: correct the Funk-Hecke derivation or downgrade the resolution claim, and add a proof of Theorem 10. In its current form I would not cite it, because trusting the central bound would mean relying on a derivation that is demonstrably wrong.\n\nFor you: worthwhile if you care about sampling methods in inverse scattering, but do not read it for the proof of the decay estimate. I would engage with it as a reviewer rather than as a citation.","headline":"A worthy extension of OSM and DSM to anisotropic electromagnetic scattering, but the decay proof in Theorem 6 has a concrete, load-bearing computational error that must be corrected before the resolution analysis can be trusted.","tokens_in":12894,"tokens_out":10147,"would_cite":false,"duration_ms":88374,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R30","35R09","65R20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that the orthogonality sampling imaging functional—a single inner product of far-field data with a known dipole wave—is positive inside an anisotropic electromagnetic scatterer and decays quadratically outside it…","keywords":["orthogonality sampling method","direct sampling method","inverse electromagnetic scattering","Maxwell's equations","anisotropic media","far field operator","factorization method","Funk-Hecke formula"],"falsifier":"Run the OSM on synthetic far-field data for a three-dimensional anisotropic scatterer with a purely real, non-dissipative contrast $P$ at a wave number $k$ that is a transmission eigenvalue, and compare the minimum of $I_{\\mathrm{OSM}}$ over the scatterer with its values just outside; if interior values are not bounded away from the exterior background, the coercivity condition in Theorem 4 is exactly what carries the detection property.","tokens_in":11860,"feed_emoji":"📡","tokens_out":11732,"duration_ms":102903,"temperature":0.7,"pith_summary":"The paper establishes a rigorous foundation for the orthogonality sampling method (OSM) in electromagnetic inverse scattering: from fixed-frequency far-field data produced by an anisotropic penetrable scatterer, the imaging functional $I_{\\mathrm{OSM}}(y_s)$—computed as a single inner product of the data with a known electric dipole wave—is proved to be bounded below by a positive constant at every sampling point $y_s$ inside the scatterer and to decay like $O(1/\\operatorname{dist}(y_s,\\Omega)^2)$ outside it. Because this two-sided behaviour holds, plotting the functional over a sampling grid recovers the scatterer's location and shape without solving an ill-posed equation at each point. The paper also proves that the direct sampling method (DSM) functional is equivalent to the OSM functional up to explicit constants, so the DSM inherits the same resolution and decay guarantees, and it proves a stability estimate showing that noise in the far-field data changes the image only to order $\\delta$. The result matters because it converts a very cheap imaging heuristic into a quantitative reconstruction tool for three-dimensional anisotropic media, with numerical examples confirming the predicted behaviour.","feed_headline":"Simple EM imaging method is proven to find scatterers","feed_subtitle":"Inside the object the imaging functional stays positive; outside it decays as the inverse square of distance.","key_machinery":"The load-bearing object is the factorization of the far-field operator $F=H^*TH$, where $H:L^2_t(S^2)\\to L^2(\\Omega)^3$ superposes incident plane waves, $H^*$ maps a volume source to its electric dipole far-field pattern, and $T=k^2P(f+v)$ records the contrast-induced source inside the scatterer. Two ingredients make the argument move: the coercivity of $T$ on the range of $H$ (a uniform positivity condition on $\\langle Tv,v\\rangle$, guaranteed when the contrast $P$ is dissipative in a uniform sense and $k$ is not a transmission eigenvalue) and the range characterization $\\varphi_{y_s}\\in\\operatorname{Range}(H^*)$ if and only if $y_s\\in\\Omega$. The Funk-Hecke formula, applied to the spherical Bessel function $j_0$, turns these operator facts into the explicit decay rate $O(1/\\operatorname{dist}(y_s,\\Omega)^2)$ for sampling points outside the scatterer.","core_discovery":"On the paper's own terms, the central discovery is that the orthogonality sampling functional $I_{\\mathrm{OSM}}(y_s)=\\|p\\cdot F\\varphi_{y_s}\\|^2_{L^2(S^2)}$ is a reliable indicator of whether a sampling point lies in the scatterer. Theorem 6 proves the two-sided behaviour: for $y_s\\in\\Omega$ the functional is bounded below by a positive constant, while for $y_s\\notin\\Omega$ it decays as $O(1/\\operatorname{dist}(y_s,\\Omega)^2)$ as $y_s$ recedes. The proof runs through the factorization $F=H^*TH$ of the far-field operator into the plane-wave superposition operator $H$, the contrast-dependent middle operator $T=k^2P(f+v)$, and the adjoint $H^*$; when $T$ is coercive on the range of $H$, the functional inherits a positive lower bound exactly on $\\Omega$, and the Funk-Hecke formula converts the decay of the spherical Bessel function $j_0$ into the explicit quadratic decay outside. The same ingredients give the equivalence $c_1 I_{\\mathrm{OSM}}\\le I_{\\mathrm{DSM}}\\le c_2\\sqrt{I_{\\mathrm{OSM}}}$, so the direct sampling functional is positive inside and decays like $O(1/\\operatorname{dist}(y_s,\\Omega)^2)$ outside as well.","pith_inferences":["Because the lower bound depends on coercivity of $T$, the theory predicts that for lossless (purely real) anisotropic contrasts the guaranteed positivity inside the scatterer can fail; testing a real-valued contrast at a transmission eigenvalue would expose exactly where the proof breaks.","The Funk-Hecke decay mechanism is specific to full spherical aperture data; an analogous calculation with finitely many incident directions would likely replace the $O(1/\\operatorname{dist}(y_s,\\Omega)^2)$ decay by an oscillatory factor, so resolution for sparse apertures may degrade more quickly than the theorem suggests.","The equivalence theorem suggests that, up to constants, the OSM image behaves like the square of the DSM image; this is consistent with the paper's numerical observation that $I_{\\mathrm{OSM}}$ resembles $(I_{\\mathrm{DSM}})^2$ more closely than $I_{\\mathrm{DSM}}$ itself.","If the coercivity conditions hold, the same factorization proof should yield explicit decay rates for other qualitative indicators, such as the linear sampling method's norm, wherever an analogous range characterization is available."],"forward_implications":["Because $I_{\\mathrm{OSM}}$ has a uniform positive lower bound on $\\Omega$ and decays quadratically outside, thresholding a plot of the functional over a sampling grid yields an estimate of the scatterer's support without solving any ill-posed linear system.","Since $I_{\\mathrm{DSM}}$ is equivalent to $I_{\\mathrm{OSM}}$ up to constants, the direct sampling method is justified with the same resolution and decay guarantees for matrix-valued anisotropic contrasts, not only scalar ones.","The stability estimate shows that a relative noise level $\\delta$ in the far-field operator changes the imaging functional by at most $|p|^2|S^2|\\|F\\|^2(\\delta^2+2\\delta)$ uniformly over all sampling points, so the reconstruction is continuous in the data.","Because the functional is uniformly positive inside $\\Omega$ and quadratically small outside, the transition zone controls how sharply an isosurface can separate the scatterer from the background, making the method's resolution quantifiable.","The numerical examples show the method still gives reasonable reconstructions at 60% and 90% noise levels and with a reduced number of incident directions, consistent with the stability estimate."],"supporting_citations":[{"why":"Supplies the coercivity theorem for the middle operator T on Range(H), the key input for the positive lower bound in Theorem 6.","marker":"[10]"},{"why":"Supplies the factorization F=H*TH and the range characterization phi_ys in Range(H*) iff ys in Omega used throughout the proof.","marker":"[14]"},{"why":"Supplies the reciprocity relation for the far-field pattern and the Funk-Hecke formula used to derive the quadratic decay.","marker":"[7]"},{"why":"Establishes the OSM-DSM equivalence in the Helmholtz case, the pattern directly transferred here to Maxwell's equations.","marker":"[17]"},{"why":"Introduces the orthogonality sampling method whose electromagnetic version is the object of the paper.","marker":"[20]"},{"why":"Provides the spectral solver used to generate the synthetic three-dimensional anisotropic scattering data for the numerical examples.","marker":"[19]"}],"fun_headline_variants":["EM orthogonality sampling: proven scatterer indicator","Orthogonality sampling functional pinpoints scatterers","Same functional proves EM scatterer location","Scatterers stand out in new EM imaging proof","EM imaging: positive inside, decays outside scatterers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire inside-scatterer lower bound rests on the contrast $P$ being dissipative enough—either $\\operatorname{Im}(P)$ uniformly positive definite, or $\\operatorname{Re}(P)+\\alpha\\operatorname{Im}(P)$ uniformly positive definite with $\\operatorname{Im}(P)$ positive semidefinite—and on $k$ not being a transmission eigenvalue; if that material condition fails, the paper does not prove the imaging functional stays positive inside the scatterer.","fun_headline_variants_meta":{"raw":{"variants":["EM orthogonality sampling: proven scatterer indicator","Orthogonality sampling functional pinpoints scatterers","Same functional proves EM scatterer location","Scatterers stand out in new EM imaging proof","EM imaging: positive inside, decays outside scatterers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000193,"raw_usage":{"total_tokens":1355,"prompt_tokens":955,"completion_tokens":400,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":571,"completion_tokens_details":{"reasoning_tokens":328}},"tokens_in":571,"tokens_out":400,"duration_ms":4522,"temperature":1.0,"reasoning_tokens":328,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:05:22.712614+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the OSM on synthetic far-field data for a three-dimensional anisotropic scatterer with a purely real, non-dissipative contrast $P$ at a wave number $k$ that is a transmission eigenvalue, and compare the minimum of $I_{\\mathrm{OSM}}$ over the scatterer with its values just outside; if interior values are not bounded away from the exterior background, the coercivity condition in Theorem 4 is exactly what carries the detection property.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the coercivity theorem for the middle operator T on Range(H), the key input for the positive lower bound in Theorem 6."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the factorization F=H*TH and the range characterization phi_ys in Range(H*) iff ys in Omega used throughout the proof."},{"cited_title":"Colton and R","cited_arxiv_id":null,"evidence_quote":"Supplies the reciprocity relation for the far-field pattern and the Funk-Hecke formula used to derive the quadratic decay."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the OSM-DSM equivalence in the Helmholtz case, the pattern directly transferred here to Maxwell's equations."},{"cited_title":"Potthast","cited_arxiv_id":null,"evidence_quote":"Introduces the orthogonality sampling method whose electromagnetic version is the object of the paper."}],"review_version":1}