{"id":"3b477fcb-250d-48af-ae5e-5f86853e12e0","arxiv_id":"2502.01037","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A direct, non-iterative reconstruction algorithm localizes a point fluorophore in a homogeneous scattering half-space by inverting asymptotic peak-time formulas and intersecting three spheres through a tetrahedron construction.","lead":"This paper derives explicit formulas that convert the peak time of a fluorescence signal measured on the surface of a scattering medium into the distance to an internal glowing point, then uses three such measurements to compute the point's location. It offers a fast, non-iterative way to localize a fluorescent target, a task usually done with expensive iterative reconstructions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Inversion (4.3) is derived for the approximate peak time, but its equivalence to the exact measured peak time is only spot-checked numerically, leaving the claimed error control unproven.","rationale":"The reader identified as the weakest assumption that the exact peak time of U_m is close enough to the approximate peak time used in the asymptotic inversions, and that this is only verified numerically. My read of the paper confirms this is the load-bearing gap: every subsequent step of the inversion — the sphere radius (4.6) and the tetrahedron vertex (4.10) — inherits errors from the λ(t) map, and λ(t) is derived from an asymptotic expansion of an approximate peak time, not from a theorem about the exact critical point. The numerical comparisons in Section 3 are reassuring but cover only the specific geometry in (3.1), and the reconstruction experiments in Section 5 do not isolate the exact-versus-approximate error from the effects of noise and of the tetrahedron geometry. The paper is otherwise coherent: the asymptotic derivations are internally consistent, the inversion algebra leading to (4.10) checks out, and the numerical feasibility is demonstrated. Because the missing piece is an error bound or a systematic numerical validation rather than a demonstrated contradiction, the appropriate verdict remains CONDITIONAL, unchanged from the reader's assessment.","tokens_in":13872,"tokens_out":21479,"duration_ms":195772,"concrete_test":"For a grid of targets with xc3 from 10 to 40 mm and S-D separations from 4 to 20 mm, compute the exact peak time of U_m by solving (1.1)-(1.2) with a high-accuracy quadrature or finite-element solver, and compare it with the approximate peak time obtained from (2.8) and (2.23). Then feed the exact peak times through (4.3) and compare the recovered λ with the true λ. If the worst-case relative error in λ exceeds twice the nominal O(t^{-1}) / O(t^{-1/2}) remainder, the claimed error control does not hold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central inversion formula (4.3) turns a measured peak time t into the distance parameter λ. However, Theorems 2.3 and 2.7 control the peak time of the approximate profile u_a^m, defined through (2.8) and (2.23), not the peak time of the exact solution U_m of (1.1)-(1.2). The link between those two peak times is only tested for the parameter values in Section 3; no bound is given for the difference δt := t_exact - t_approx. Since the inversion has slope dλ/dt ≈ k^{1/2} in the small-lifetime case and ≈ (k-ℓ^{-1})^{1/2} in the large-lifetime case, a systematic deviation δt shifts every reconstructed λ by roughly that slope times δt. The shifted λ then enters all three sphere radii (4.6) and the tetrahedron vertex (4.10). For ℓ ≫ 1 the nominal remainder in (4.3) is only O(t^{-1/2}), so for t ~ 10^3 ps it already allows errors of several percent; the noise-free experiments show exactly this magnitude. The paper therefore does not establish that the reported 1-5% errors are covered by the asymptotic error terms, nor that the scheme is accurate outside the tested configurations.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a direct, non-iterative inversion scheme for point-target fluorescence diffuse optical tomography (FDOT) using measured peak times of temporal response functions. The authors derive asymptotic expansions for an approximate peak time—defined as the critical point of an asymptotic profile of the emission solution—in two regimes: small fluorescence lifetime (Theorem 2.3) and large fluorescence lifetime (Theorem 2.7). These expansions are inverted to express the distance parameter λ in terms of peak time (Section 4, Eq. (4.3)), which determines a sphere on which the target lies (Eq. (4.5)). Using three source-detector pairs, the target is reconstructed as the vertex of a tetrahedron whose edges are the sphere radii (Eqs. (4.7)–(4.10)). Numerical experiments in Section 5 report noise-free reconstruction errors of roughly 1–5% and demonstrate robustness to time-jitter noise.","tokens_in":14148,"tokens_out":21204,"duration_ms":171997,"significance":"If the error between the approximate and the exact peak time could be rigorously controlled, this would be a valuable contribution: the scheme is explicit, avoids iterative optimization, uses peak-time measurements that are relatively robust to noise, and extends earlier peak-time formulas to include the extrapolation length β and nonzero fluorescence lifetime. The asymptotic formulas are closed-form and the tetrahedron reconstruction is straightforward to implement. The paper also provides numerical verification across several parameter sets. However, the central claim currently rests on an unproven link between the exact and approximate peak times and on a remainder estimate in the inversion that is not correct at the stated order, so the significance is conditional on those issues being resolved.","major_comments":[{"comment":"The asymptotic theorems are stated for an approximate peak time defined as the critical point of the profile u_a^m, not for the exact peak time t_peak = argmax_t U_m used in the Inverse Problem of Section 1. Lemma 2.1 (Eq. (2.5)) gives only a pointwise approximation u_m = u_a^m(1+O(λ^{-1})), with no control of the derivative or of the argmax. Therefore applying the inversion (4.3) to a measured t_peak has no proven error bound. The numerical comparisons in Section 3, while supportive, are spot checks and do not supply the missing estimate. Since the entire inversion scheme is driven by measured peak times, this gap is load-bearing.","section":"Section 2, Definitions 2.2/2.6 and Theorems 2.3/2.7; Section 4, Eq. (4.3)"},{"comment":"The replacement αλ = ᾶt + O(t^{-1/2}) is not correct at the stated order. From λ = (k−ℓ^{-1})^{1/2} t + O(t^{1/2}) one has λ^{1/2} = (k−ℓ^{-1})^{1/4} t^{1/2} + O(1), so αλ^2 = ᾶt^2 + O(t^{-1/2}) and hence αλ = ᾶt + O(t^{-1/2}/√log t). When the term −αλ t^{1/2} is expanded, the induced error in λ is O(t^{1/2} · t^{-1/2}/√log t) = O(1/√log t), not O(t^{-1/2}). Thus the final inversion formula (4.3) for ℓ ≫ 1 has a remainder that is only o(1), and its convergence is only logarithmic. At the finite times used in Section 5 (t ~ 10^2–10^3 ps), 1/√log t is of order 0.1–0.4, so the asymptotic analysis does not by itself explain the observed few-percent accuracy; the numerical success may rely on small constants not captured by the big-O estimates.","section":"Section 4, Eq. (4.2)"},{"comment":"The algebraic reconstruction from the three radii r, r1, r2 is presented without any error propagation or conditioning analysis. The vertex formula (4.10) involves square roots and differences of fourth powers; small perturbations in the radii can be amplified, especially for nearly degenerate tetrahedra. The numerical tables cover only a few geometries and do not report the radii or the angles θ1, θ2, so the reader cannot assess the stability of the reconstruction. A quantitative sensitivity estimate is needed to support the claim that three peak-time measurements determine the target position accurately.","section":"Section 5, Eq. (4.10)"}],"minor_comments":[{"comment":"The notation in the β-dependent term is ambiguous; the denominator should be written consistently as √k + β√(vD), and the same expression should be used in (1.5), (2.9), and (4.1).","section":"Eqs. (1.5), (2.9), (4.1)"},{"comment":"Lemma 2.5 is quoted from the submitted manuscript [2] and its proof is not included. Please make the paper self-contained by providing a proof or a detailed proof sketch, especially because the lemma is the basis for the large-lifetime equation (2.23).","section":"Section 2, Lemma 2.5"},{"comment":"The definition of ᾶt requires the argument of the logarithm to lie in (0,1) for ᾶt to be real; this follows from the condition in (4.3) but should be stated explicitly.","section":"Section 4, Eq. (4.4)"},{"comment":"Reporting the radii r, r1, r2 and the angles θ1, θ2 for the numerical examples would help the reader evaluate the conditioning of the tetrahedron reconstruction and reproduce the results.","section":"Section 5"},{"comment":"There are several typographical and grammatical issues, e.g., 'unknwon' in Section 6, 'the proof the Theorem 2.7' at the end of the proof of Theorem 2.7, and redundant phrases such as 'is as in' in Section 3.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper leans heavily on the authors' own prior work [1,2], with [2] still in submission; this is not a circularity, but the reliance on an unavailable manuscript for a key lemma should be addressed. The main unresolved issue—the missing bound between exact and approximate peak times—is not a disagreement with the field's consensus but a genuine gap in the proof. The formula (4.2) error is a concrete technical mistake that should be fixed. The paper otherwise fits the scope of math.AP, and the explicit scheme is potentially useful, so a major revision is warranted rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a genuine closed-form inversion for point-target time-domain FDOT in a homogeneous half-space. The asymptotic expansions in the small- and large-lifetime regimes are explicit, the numerics look encouraging, and the tetrahedron reconstruction is simple and clearly described. The main soft spot is that the paper does not bound the difference between the exact peak time and the approximate peak time used in the inversion.\n\nWhat's new: Theorems 2.3 and 2.7 give asymptotic expansions for the approximate peak time, generalizing the authors' earlier equation (1.8) to (1.9) with a boundary-term correction. Equations (4.1)-(4.3) invert those expansions to produce a direct formula for the distance parameter lambda(t), and (4.5)-(4.10) turn three S-D pairs into the target location. I have not seen this explicit tetrahedron construction elsewhere. The derivations are internally consistent, and the paper is honest about the approximations.\n\nWhere it is soft: the stress-test point is the one that matters. The theorems control the approximate peak time, defined as the solution to (2.8) or (2.23), not the measured peak time of U_m. The gap between them is only verified numerically in Section 3 for a finite set of parameters. Since (4.3) has slope O(1) in t, a systematic offset delta-t shifts all three radii in (4.6) and the reconstructed vertex in (4.10). For large lifetimes the nominal remainder in (4.3) is only O(t^{-1/2}), which for t ~ 10^3 ps leaves errors of a few percent; the noise-free RelErr values in Tables 1 and 3 are consistent with that. The authors do not claim rigorous control, so this is a soft spot, not a fatal flaw. The noise experiments are single-sample per level, and no forward-solver code is shipped, so the figures cannot be independently reproduced. Those are addressable in revision.\n\nWho this is for: anyone who wants a fast, non-iterative reconstruction of a single fluorescent point target in time-domain FDOT. It deserves refereeing because the core asymptotics are worked out carefully and the method is transparent. I would send it to a referee with a request to either prove a bound on the exact-approximate peak-time gap or demonstrate numerically that the gap stays small across a wider parameter range, and to add repeated noise runs. Engage with it; this is a serious paper that needs revision rather than rejection.","headline":"Direct closed-form FDOT point-target inversion with clean asymptotics, but the exact-to-approximate peak-time gap is numerically spot-checked, not bounded.","tokens_in":14676,"tokens_out":2958,"would_cite":false,"duration_ms":32229,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R30","35K20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that three peak-time measurements suffice to reconstruct a fluorescent point target by converting each peak time into a sphere radius and taking the target as the tetrahedron vertex.","keywords":["FDOT","peak time","asymptotic analysis","reconstruction algorithm","fluorescence diffuse optical tomography","inverse problem","time-domain measurement","point target"],"falsifier":"Compute the true peak time from a full numerical solution of the coupled diffusion equations for a configuration where $\\lambda$ is only moderately large, such as a target depth comparable to the source-detector separation, then run the inversion (4.3)-(4.10). If the reconstructed vertex error grows beyond the $O(t^{-1})$ or $O(t^{-1/2})$ prediction, the asymptotic peak-time approximation is breaking. A sharper test is to compare the positive solution of (2.23) with the true peak time over a grid of target depth, lifetime, and absorption values and locate where the relative difference exceeds a few percent.","tokens_in":13672,"feed_emoji":"🎯","tokens_out":5339,"duration_ms":48300,"temperature":0.7,"pith_summary":"This paper proposes a direct, non-iterative way to locate a single fluorescent point target in a highly scattering half-space from the peak times of time-domain fluorescence measurements. The key claim is that each measured peak time determines a distance parameter $\\lambda$, which puts the target on an explicit sphere centered at the midpoint of that source-detector pair. Combining three source-detector pairs gives three sphere radii, and these radii are the edge lengths of a tetrahedron whose missing vertex is the target; the paper gives explicit formulas for that vertex. If valid, this would replace iterative image reconstruction in FDOT with a closed-form calculation using only three boundary measurements.","feed_headline":"Three peak times pinpoint a hidden fluorescence target","feed_subtitle":"Each measured peak time becomes a sphere radius; three radii build a tetrahedron whose vertex is the target.","key_machinery":"The engine of the argument is the asymptotic approximate peak time. For small fluorescence lifetime $\\ell$, the approximate peak solves $P(t) = \\ell(P'(t)+P(t)^2)$ with $P(t) = -k - \\frac32 t^{-1} + \\lambda^2 t^{-2} - \\frac{2\\beta vD}{x_{c3}+\\beta vD t}$, giving $t_p = k^{-1/2}\\lambda - \\frac74 k^{-1} + \\frac{\\sqrt{k}}{\\sqrt{k}+\\beta\\sqrt{vD}} k^{-1} + \\ell + O(\\lambda^{-1})$. For large $\\ell$, the approximate peak solves equation (2.23), giving $t_p = (k-\\ell^{-1})^{-1/2}\\lambda + (k-\\ell^{-1})^{-3/4}\\alpha_\\lambda \\lambda^{1/2} + O(\\lambda^{-1/2})$. Inverting these expansions yields the explicit $\\lambda(t)$ in (4.3), and the sphere identity turns $\\lambda$ into the tetrahedron edge length. The inversion scheme reduces the inverse problem to a purely geometric construction: three known sphere radii around three known midpoints determine the target as the tetrahedron vertex.","core_discovery":"The central discovery is that the peak time of the emission signal can be inverted explicitly into the target-distance parameter $\\lambda(t)$: for small fluorescence lifetime $\\ell$, $\\lambda(t) = k^{1/2}t + \\frac{7}{4}k^{-1/2} - \\frac{\\sqrt{k}}{\\sqrt{k}+\\beta\\sqrt{vD}}k^{-1/2} - \\ell k^{1/2}$, and for large $\\ell$, $\\lambda(t) = (k-\\ell^{-1})^{1/2}t - \\tilde\\alpha_t t^{1/2} + \\frac{\\tilde\\alpha_t^2}{2}(k-\\ell^{-1})^{-1/2}$, up to remainder terms. Once $\\lambda$ is known, the identity $\\lambda^2 = (|x_d-x_c|^2 + |x_s-x_c|^2)/(2vD)$ places the unknown target $x_c$ on the sphere centered at $(x_d+x_s)/2$ with radius $r = \\sqrt{vD\\lambda^2 - |x_d-x_s|^2/4}$. Three source-detector pairs give three such radii, which form the edges of a tetrahedron, and formulas (4.10) recover the target as the fourth vertex. In noise-free numerical tests the reconstruction errors are about 1.4% to 5.3% relative error for targets at depth 20-30 mm, and the depth estimate remains acceptable at 5% time-jitter noise.","pith_inferences":["A natural testable extension would feed experimentally measured peak times from a tissue phantom into (4.3)-(4.10) and compare against known embedded fluorophore positions; the asymptotic error prediction says the error should shrink as source-detector-target distances grow.","Since the method needs only peak times rather than full temporal profiles, it might combine naturally with low-cost time-gated detectors; in practice the main experimental bottleneck is likely accurate calibration of the optical parameters $k$, $D$, $\\mu_a$, $\\beta$, and the lifetime $\\ell$, since all enter $\\lambda(t)$.","The tetrahedron construction suggests a geometric consistency check for redundant data: with more than three source-detector pairs, the reconstructed vertex should lie near the intersection of all the spheres, and disagreement among pairs could serve as a data-quality diagnostic.","If extended to non-point or multiple targets, the single-sphere argument likely breaks down unless the target distribution's peak-time signature is separated by lifetime or gating; the paper itself lists non-point targets and curved boundaries as future work."],"forward_implications":["Only three source-detector pairs are needed to recover a point target, and no iterative PDE solves or initial guesses are required.","The reconstruction is fully explicit: formulas (4.10) give the target coordinates directly from the three radii and two rotation angles.","The asymptotic analysis supplies error orders $O(t^{-1})$ for small lifetime and $O(t^{-1/2})$ for large lifetime, which the numerical experiments reflect as 1-5% relative error without noise.","Deeper targets at 30 mm depth reconstruct at least as accurately as shallower targets at 20 mm in the tests, consistent with the large-distance assumption $\\lambda \\gg 1$.","Time-jitter noise of 5% degrades the in-plane accuracy but still leaves a usable depth estimate, suggesting the method tolerates moderate measurement noise.","Because the work is a direct corollary of the paper's own claims, these consequences stand or fall with the asymptotic peak-time inversion."],"supporting_citations":[{"why":"Supplies the small-lifetime asymptotic peak time and the Gaussian profile approximation used in Definition 2.2 and Theorem 2.3.","marker":"[1]"},{"why":"Supplies the large-lifetime peak-time equation (1.8) that is generalized to (1.9), and the saddle-point integration estimate in Lemma 2.5.","marker":"[2]"},{"why":"Gives the zero-lifetime peak-time relation that the present work extends to include finite lifetime and boundary reflection terms.","marker":"[3]"},{"why":"Provides the derivation of the zero-lifetime emission solution (2.1), which is the starting point for the asymptotic profile approximation.","marker":"[4]"},{"why":"States the coupled diffusion model (1.1)-(1.2) for excitation and emission that defines the forward problem.","marker":"[8]"}],"fun_headline_variants":["Peak times reconstruct fluorescence target via tetrahedron","Sphere and vertex: peak-time inversion of fluorescence targets","Three peak times pinpoint hidden fluorescence target","Direct peak-time inversion locates fluorescence target"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the measured peak time is close enough to the peak of the asymptotic profile $u^a_m$ that the expansions in (2.9) and (2.25), and their inversions in (4.3), apply; this closeness is checked numerically in Section 3 for selected parameters, but it is not proved for all source-detector-target configurations.","fun_headline_variants_meta":{"raw":{"variants":["Peak times reconstruct fluorescence target via tetrahedron","Sphere and vertex: peak-time inversion of fluorescence targets","Three peak times pinpoint hidden fluorescence target","Direct peak-time inversion locates fluorescence target"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000282,"raw_usage":{"total_tokens":1671,"prompt_tokens":948,"completion_tokens":723,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":666}},"tokens_in":564,"tokens_out":723,"duration_ms":6978,"temperature":1.0,"reasoning_tokens":666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T16:49:18.832078+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the true peak time from a full numerical solution of the coupled diffusion equations for a configuration where $\\lambda$ is only moderately large, such as a target depth comparable to the source-detector separation, then run the inversion (4.3)-(4.10). If the reconstructed vertex error grows beyond the $O(t^{-1})$ or $O(t^{-1/2})$ prediction, the asymptotic peak-time approximation is breaking. A sharper test is to compare the positive solution of (2.23) with the true peak time over a grid of target depth, lifetime, and absorption values and locate where the relative difference exceeds a few percent.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the small-lifetime asymptotic peak time and the Gaussian profile approximation used in Definition 2.2 and Theorem 2.3."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the large-lifetime peak-time equation (1.8) that is generalized to (1.9), and the saddle-point integration estimate in Lemma 2.5."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the zero-lifetime peak-time relation that the present work extends to include finite lifetime and boundary reflection terms."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the derivation of the zero-lifetime emission solution (2.1), which is the starting point for the asymptotic profile approximation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"States the coupled diffusion model (1.1)-(1.2) for excitation and emission that defines the forward problem."}],"review_version":1}