{"id":"1630078e-c253-414a-88eb-2c99b894e8d3","arxiv_id":"1908.08683","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper proves compactness and non-uniqueness for the inverse eddy current conductivity problem and develops a Sobolev-gradient nonlinear conjugate gradient method that reconstructs inclusions in 3D.","lead":"The paper proves that recovering the electrical conductivity of materials from low-frequency eddy current measurements is mathematically ill-posed and non-unique, then builds a regularized optimization algorithm that reconstructs inclusions in 3D box domains. A smart generalist might read it because it supplies rigorous foundations for a widely used nondestructive testing technique.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 3.2 proves non-uniqueness only for the induced source iωσE(σ), not for the conductivity σ; no pair of distinct conductivities with identical boundary data is constructed, so the abstract's non-uniqueness claim is unproven.","rationale":"The reader's weakest assumption identifies exactly the load-bearing gap: Theorem 3.2 addresses the induced source σE(σ), not the conductivity σ. My reading confirms that the proof decomposes a single source into radiating and non-radiating parts and concludes that the source is non-unique, but never constructs two conductivities with identical tangential boundary data. This is not a mere presentational issue: the abstract and Section 3.2 title claim non-uniqueness of the recovery process for σ, and the inverse-source result cannot be transferred without additional structure. The compactness lemma, regularity theorem, and regularized minimization analysis appear plausible and are independent of this gap, so the appropriate disposition remains CONDITIONAL rather than REJECT. The proposed construction test directly probes whether the missing transfer can be made; if it fails, the paper should either weaken the advertised claim or provide a separate non-uniqueness proof for σ.","tokens_in":20635,"tokens_out":11825,"duration_ms":121908,"concrete_test":"Test the transfer step in Theorem 3.2 by attempting an explicit construction: choose a nonzero J2∈W⊥ with support in Ωc, let E=E0+E(J2) solve (3.3) with source J2, and define σ1=J2/(iωE) wherever E≠0 (and zero elsewhere). Verify whether σ1 is real-valued, nonnegative, belongs to H1_0(Ωc), and satisfies E=E(σ1) for the original source Js. If no such admissible σ1 can be produced from any J2∈W⊥, then the inverse-source non-uniqueness does not imply non-uniqueness of conductivity recovery, and the advertised claim requires a separate proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2 is titled 'Non-uniqueness of the recovery of the conductivity', and the abstract advertises non-uniqueness of the recovery process. The actual Theorem 3.2 proves only that for a fixed σ satisfying n×(E(σ)-E0)≠0 on Γ, the secondary source iωσE(σ) cannot be determined uniquely from n×E on Γ. This is an inverse-source non-uniqueness statement. It does not yield two distinct conductivities with the same boundary data. To close the gap one would need, given the decomposition iωσE(σ)=J1+J2 with nonzero J2∈W⊥, to produce σ1≠σ such that iωσ1E(σ1)=J1 and n×E(σ1)=n×E(σ) on Γ. In particular σ1 would have to satisfy the fixed-point equation σ1=J1/(iωE(σ1)), be real-valued, nonnegative, compactly supported in Ωc, and lie in H1_0(Ωc). None of these properties is shown. The compactness result and the regularized optimization analysis are not affected, but the advertised non-uniqueness of conductivity recovery is not established by the paper's arguments.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the inverse eddy current problem of recovering the conductivity σ in a bounded three-dimensional conductor from tangential electric field measurements n×E on a boundary part Γ. The authors first prove well-posedness of the forward saddle-point formulation, establish a regularity result for E(σ)−Eobs in the air subdomain Ω0, and prove compactness of the map σ↦n×E(σ)|Γ. They then claim non-uniqueness of the conductivity recovery, formulate a Tikhonov-regularized constrained optimization problem, prove existence and stability of minimizers, derive the adjoint-based gradient, discretize the problem with edge elements, and propose a nonlinear conjugate gradient method with a Sobolev gradient. Numerical experiments on one and two inclusions are reported.","tokens_in":20834,"tokens_out":15307,"duration_ms":154003,"significance":"If the stated results are taken at face value, the paper contributes a fairly complete analysis pipeline for an eddy-current inversion problem: forward well-posedness, a regularity result that justifies using an L2 misfit on Γ, compactness of the forward map, existence/stability of regularized minimizers, and a concrete adjoint/NLCG algorithm with numerical demonstrations. The main advertised theoretical novelty, however, is the non-uniqueness of the conductivity recovery, and that claim is not established by the arguments in Section 3.2. The paper actually proves non-uniqueness of the equivalent source iωσE(σ) in the auxiliary inverse source problem, not non-uniqueness of σ itself. The optimization analysis and numerical framework are valuable and do not depend on the overclaimed conductivity non-uniqueness, so the paper is salvageable with a precise reformulation of the ill-posedness statement.","major_comments":[{"comment":"Theorem 3.2 overreaches in a load-bearing way. The proof shows only that the secondary source Je = iωσE(σ) can be decomposed as J1 + J2 with J2 ∈ W⊥, J2 ≠ 0, and that by Theorem 3.1 the component J2 produces no tangential field on Γ. This means the inverse source problem (3.4) does not have a unique source. It does not produce two distinct conductivities σ and σ1 with n×E(σ) = n×E(σ1) on Γ. To obtain conductivity non-uniqueness one would need to construct σ1 ≠ σ such that iωσ1E(σ1) = J1 and n×E(σ1) = n×E(σ) on Γ, and one would need σ1 to be real-valued, nonnegative, compactly supported in Ωc, and in H1_0(Ωc). None of these properties is shown, and the fixed-point equation σ1 = J1/(iωE(σ1)) is not analyzed. Consequently the title of Section 3.2, the abstract's claim of non-uniqueness of the recovery process, and the concluding remarks overstate the result. The authors should either prove conductivity non-uniqueness or explicitly restrict the claim to non-uniqueness of the equivalent source in the inverse source problem.","section":"§3.2, Theorem 3.2"},{"comment":"The proof of Theorem 2.2 contains a missing argument that is load-bearing for the regularity result. It is asserted that 'with the arguments in Theorem 3.1 of Section 3, we know that ∇φ is a non-radiating source, then Eφ|Ω0 = 0.' To apply the non-radiating criterion of Theorem 3.1, one needs to prove that ∇φ ∈ W⊥, where W is defined only later in Section 3.2. This is not shown in the proof, and it is not immediate from the displayed decomposition Je = J0 + ∇φ. A short argument using that every u ∈ W satisfies ∇·u = 0 in Ωc (from the equation in the definition of W) would close the gap, but as written the proof of Theorem 2.2 is incomplete. Since Theorem 2.2 is used to justify the L2 misfit functional and is subsequently invoked in Lemmas 2.3, 3.1, 3.2, and Theorem 3.3, this gap should be repaired explicitly.","section":"§2.2, proof of Theorem 2.2"}],"minor_comments":[{"comment":"The statement says σn → σ* in L2(Γ), but the proof and the intended compactness argument require convergence in L2(Ωc); this should be corrected.","section":"Lemma 3.1"},{"comment":"The proof cites 'Lemma 2.2' for the regularity of E(σn)−E(σ*) in Ω0, but Lemma 2.2 is the uniqueness result; the regularity bound comes from Theorem 2.2 and Lemma 2.3. This citation should be fixed.","section":"Lemma 3.2"},{"comment":"The statement writes H1_Γ(Ωc) for the multiplier φ; the correct space is H1_Γ(Ω0) as used in the proof.","section":"Theorem 2.1"},{"comment":"The phrase 'homogeneous eigenfunction corresponding to imaginary eigenvalue iωµσ0' is sign-inconsistent with the definition of W, which has +iωσ0u = 0; this should be corrected or clarified.","section":"Theorem 3.2, proof"},{"comment":"The expression Je = iωµσE(σ) is inconsistent with the earlier definition Je = iωσE(σ) in Section 3.2; the factor µ should be removed or the convention explained.","section":"Remark 3.2"},{"comment":"The numerical section reports only visual comparisons and iteration counts; a quantitative reconstruction error with respect to the true σ and a clear stopping criterion would strengthen the feasibility claim.","section":"§5, Numerical experiments"}],"recommendation":"major_revision","confidential_remarks":"The paper's optimization and numerical content is solid, but the advertised non-uniqueness of conductivity recovery is not proven. The fix is feasible within the manuscript's scope: either prove conductivity non-uniqueness under additional assumptions or reframe the contribution as non-uniqueness of the equivalent inverse source problem. I recommend major revision rather than rejection, since the regularized optimization results and numerical framework remain valid regardless of the outcome of that correction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper and the stress-test note. I think the note is right, and it names the main problem: the advertised non-uniqueness of conductivity recovery is not established. Theorem 3.2 proves non-uniqueness of the induced source iωσE(σ) in the inverse source problem, not non-uniqueness of σ. The reduction in Section 3.2 drops the constraint that the source must have the form iωσE(σ) for some admissible conductivity. To close the gap you would need to construct σ1≠σ with iωσ1E(σ1)=J1 and the same trace on Γ, preserving real-valuedness, support, and H1 regularity. None of that is shown. So the abstract's \"non-uniqueness of the recovery process\" overstates what the proof gives. This is a genuine gap in one advertised contribution, but it is not a fatal one.\n\nWhat is good: the compactness of σ→n×E(σ) from H1_0(Ωc) into L2(Γ) is a real result for this eddy current model, and the H1/2 regularity for E(σ)-Eobs on Ω0 is a useful justification for using the L2 boundary misfit rather than the intractable H^{-1/2}(Div;Γ) norm. The regularized constrained optimization and its existence/stability analysis are standard but carefully adapted to the complex-valued saddle-point structure. The real/imaginary Lagrangian, the adjoint gradient, and the Sobolev-gradient NLCG are competently put together. The numerics demonstrate feasibility on nontrivial 3D examples, especially the faster convergence with the Sobolev gradient. The citation pattern is fine: [26], [12], and [10] are genuinely used for non-radiating sources, singularities, and augmented Lagrangian ideas; self-citations [7], [8], [18] support auxiliary well-posedness and preconditioning.\n\nSofter spots: Theorem 2.2 asserts that ∇φ is a non-radiating source with a pointer to the arguments in Theorem 3.1; that step needs a real proof or a precise citation. The numerical section is feasibility-only—no quantitative reconstruction errors or systematic regularization parameter choice—fine for a methods paper, but not evidence of practical accuracy.\n\nWho this is for: researchers working on inverse eddy current or PDE-constrained optimization for nondestructive testing. It deserves a serious referee. The non-uniqueness claim needs to be either proved with a genuine two-conductivity construction or rewritten as inverse-source non-uniqueness with the conductivity question left open.","headline":"The paper has solid compactness, regularity, and optimization work, but its headline non-uniqueness claim only proves non-uniqueness of the induced source, not of the conductivity itself.","tokens_in":21412,"tokens_out":4643,"would_cite":true,"duration_ms":46415,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R30","35B30"],"pacs":[],"model":"deepseek-v4-flash","headline":"The inverse eddy current problem is ill-posed and non-unique, and the paper builds a regularized optimization that recovers inclusions numerically.","keywords":["inverse eddy current","ill-posedness","non-uniqueness","regularity","stability","Lagrangian","adjoint problem","nonlinear conjugate gradient"],"falsifier":"A concrete test: pick two distinct admissible conductivities $\\sigma_1,\\sigma_2$ (for instance two different inclusions in $\\Omega_c$), solve the forward problem (2.6) for each, and compare $n\\times E(\\sigma_1)$ and $n\\times E(\\sigma_2)$ on $\\Gamma$. Exact equality would confirm true non-uniqueness of conductivity; a proof that $\\sigma \\mapsto \\sigma E(\\sigma)$ is injective on the admissible class would show only the source is ambiguous and leave conductivity uniqueness open.","tokens_in":20368,"feed_emoji":"🧲","tokens_out":7241,"duration_ms":61252,"temperature":0.7,"pith_summary":"This paper studies whether the spatial distribution of electrical conductivity inside a conductor can be recovered from measurements of the tangential electric field on the boundary, the setting of eddy-current nondestructive testing. It establishes two structural obstructions: the map sending conductivity to the boundary tangential field is compact (Lemma 3.2), and the recovery is non-unique because the induced source term $i\\omega\\sigma E(\\sigma)$ always contains a non-radiating part invisible to boundary data (Theorem 3.2). It then shows that adding a Tikhonov-type regularization turns the problem into a constrained minimization with provable existence and stability of minimizers. A Lagrangian/adjoint formulation yields the gradient, and a nonlinear conjugate gradient method with a Sobolev gradient recovers separated inclusions in numerical tests, with or without noise.","feed_headline":"Surface fields alone can't pin down conductivity in eddy current scans","feed_subtitle":"Even noise-free surface measurements leave the conductivity ambiguous, so recovery is recast as a regularized optimization.","key_machinery":"The load-bearing object is the orthogonal decomposition of the secondary source space $L^2(\\Omega_c)^3 = W \\oplus W^\\perp$, where $W$ consists of fields solving the homogeneous eddy-current equation inside the conductor with a vanishing tangential trace condition. Sources lying in $W^\\perp$ are non-radiating: they produce $n\\times E=0$ on the interface and measurement surface, so any source with a nonzero $W^\\perp$ component is invisible to boundary data. The compactness of the forward map comes from the compact embedding $H^1(\\Omega_c)\\hookrightarrow L^2(\\Omega_c)$ plus the continuity estimate linking conductivity differences to $H^{1/2}$ differences of the field in the air region, a regularity gain supplied by non-radiating-source arguments. Together these identify exactly why the inverse problem needs regularization.","core_discovery":"On the paper's own terms, the central discovery is that the inverse eddy current problem based on tangential electric boundary data is fundamentally ill-posed: the forward operator $\\sigma \\mapsto n\\times E(\\sigma)$ is compact from $H^1_0(\\Omega_c)$ into $L^2(\\Gamma)$, so bounded sets of conductivities produce precompact sets of data, and the secondary source $i\\omega\\sigma E(\\sigma)$ is never uniquely determined when the measured field differs from the background field. The proof decomposes the source space as $L^2(\\Omega_c)^3 = W \\oplus W^\\perp$; sources in $W^\\perp$ are non-radiating (Theorem 3.1), and Theorem 3.2 shows both summands are nonzero, so the invisible component prevents unique recovery. The rest of the paper builds a well-posed surrogate: minimize $\\Phi_\\alpha(\\sigma)=\\frac12\\|n\\times(E(\\sigma)-E_{\\rm obs})\\|^2_{L^2(\\Gamma)}+\\frac{\\alpha}{2}\\|\\nabla\\sigma\\|^2_{L^2(\\Omega_c)}$, prove existence and stability of its minimizers, and solve it numerically by a feasible Lagrangian, adjoint gradient, and NLCG with Sobolev gradient, with numerical reconstructions of separated inclusions.","pith_inferences":["The stated non-uniqueness is for the source $\\sigma E(\\sigma)$, not for $\\sigma$ itself; whether two distinct conductivities can produce identical boundary data is left open, since the map $\\sigma\\mapsto \\sigma E(\\sigma)$ could conceivably be injective even though the source is not unique.","The $W/W^\\perp$ splitting suggests a natural definition of 'detectable' conductivity perturbations: a perturbation is invisible to first order if its secondary source lies mostly in $W^\\perp$; this could be tested numerically by optimizing a perturbation to minimize the boundary trace.","Because the splitting depends on frequency $\\omega$, multi-frequency measurements may shrink the invisible component and could restore uniqueness even when single-frequency data cannot; this is an extension the paper does not pursue.","The same Lagrangian/adjoint machinery would apply to recovering other parameters such as magnetic permeability or to multi-frequency data, since the structure of the state and adjoint systems is unchanged."],"forward_implications":["Because the data-to-conductivity map is compact, small measurement errors can be amplified arbitrarily; any numerical inversion must be regularized, not just sampled finely.","The non-radiating part of the induced source means perfect noise-free boundary data still cannot single out a unique conductivity distribution in general.","The regularized functional has at least one minimizer and minimizers depend stably on the data, so the optimization formulation is a legitimate surrogate for the ill-posed inversion.","The $H^{1/2}$ regularity of the field difference in the air region justifies using the computable $L^2(\\Gamma)$ misfit instead of the harder $H^{-1/2}(\\mathrm{Div};\\Gamma)$ trace norm.","The NLCG algorithm with the Sobolev gradient reconstructs locations and sizes of separated inclusions, and its convergence is markedly faster than with the plain $L^2$ gradient."],"supporting_citations":[{"why":"Supplies the non-radiating source decomposition and inverse-source uniqueness technique used in Section 3.2.","marker":"[26]"},{"why":"Provides the regularity and imbedding result (Theorem 6.1) used to prove the $H^{1/2}$ regularity of the field difference in the air region.","marker":"[12]"},{"why":"Defines the $H^{-1/2}(\\mathrm{Div};\\Gamma)$ tangential trace space and the trace theorems used throughout the formulation.","marker":"[6]"},{"why":"Establishes coercivity and well-posedness of the forward eddy-current saddle-point problem and supplies the edge-element background.","marker":"[7]"},{"why":"Gives the standard regularization existence and stability arguments adapted in Theorems 3.3 and 3.4.","marker":"[10]"},{"why":"Provides the analogous regularized-minimizer analysis for inverse Maxwell and elliptic problems.","marker":"[13]"},{"why":"Introduces the Sobolev gradient preconditioning idea used in Section 4.6.","marker":"[18]"},{"why":"Supplies the edge-element and trace-space background for the finite element discretization.","marker":"[20]"}],"fun_headline_variants":["Eddy current inversion: conductivity invisible to surface data","Non-unique eddy current inverse problem solved via optimization","Hidden sources break uniqueness in eddy current conductivity recovery","Ill-posed eddy current recovery recast as regularized optimization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reduction from conductivity to the secondary source $\\sigma E(\\sigma)$ is assumed lossless: Theorem 3.2 proves only that the source is not unique, and the paper takes that to mean the conductivity itself cannot be recovered, without exhibiting two distinct conductivities with identical boundary data.","fun_headline_variants_meta":{"raw":{"variants":["Eddy current inversion: conductivity invisible to surface data","Non-unique eddy current inverse problem solved via optimization","Hidden sources break uniqueness in eddy current conductivity recovery","Ill-posed eddy current recovery recast as regularized optimization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000689,"raw_usage":{"total_tokens":3145,"prompt_tokens":989,"completion_tokens":2156,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":2089}},"tokens_in":605,"tokens_out":2156,"duration_ms":15526,"temperature":1.0,"reasoning_tokens":2089,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:31:59.558903+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: pick two distinct admissible conductivities $\\sigma_1,\\sigma_2$ (for instance two different inclusions in $\\Omega_c$), solve the forward problem (2.6) for each, and compare $n\\times E(\\sigma_1)$ and $n\\times E(\\sigma_2)$ on $\\Gamma$. Exact equality would confirm true non-uniqueness of conductivity; a proof that $\\sigma \\mapsto \\sigma E(\\sigma)$ is injective on the admissible class would show only the source is ambiguous and leave conductivity uniqueness open.","supporting_citations":[{"cited_title":"Rodriguez, J","cited_arxiv_id":null,"evidence_quote":"Supplies the non-radiating source decomposition and inverse-source uniqueness technique used in Section 3.2."},{"cited_title":"Costabel, M","cited_arxiv_id":null,"evidence_quote":"Provides the regularity and imbedding result (Theorem 6.1) used to prove the $H^{1/2}$ regularity of the field difference in the air region."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the $H^{-1/2}(\\mathrm{Div};\\Gamma)$ tangential trace space and the trace theorems used throughout the formulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes coercivity and well-posedness of the forward eddy-current saddle-point problem and supplies the edge-element background."},{"cited_title":"Chen and J","cited_arxiv_id":null,"evidence_quote":"Gives the standard regularization existence and stability arguments adapted in Theorems 3.3 and 3.4."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the analogous regularized-minimizer analysis for inverse Maxwell and elliptic problems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Sobolev gradient preconditioning idea used in Section 4.6."},{"cited_title":"Monk, Finite element methods for Maxwell’s equations, Oxford University Press, 2003","cited_arxiv_id":null,"evidence_quote":"Supplies the edge-element and trace-space background for the finite element discretization."}],"review_version":1}