pith. sign in

arxiv: 2509.06183 · v2 · submitted 2025-09-07 · 🧮 math.AP

Forward and inverse problems of a semilinear transport equation

Kui Ren , Yimin Zhong This is my paper

Pith reviewed 2026-05-18 17:55 UTC · model grok-4.3

classification 🧮 math.AP
keywords semilinear transport equationwell-posednessinverse problemsstability estimatesradiative transportphotoacoustic imagingnonlinear absorptionweighted norms
0
0 comments X

The pith

A semilinear radiative transport model is well-posed for arbitrary boundary data and yields stable recovery of the nonlinear absorption coefficient from internal measurements.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proves that the forward problem for a semilinear transport equation admits solutions when boundary data are not restricted to be small, extending earlier local results to a global theory. It then derives stability estimates for the inverse problem of recovering the absorption coefficient that depends on the angular average of the solution, using internal data. These estimates introduce a weighted norm that controls boundary contributions and thereby produces a single L1 stability theory that covers both the transport and diffusion regimes. The setup is motivated by the need to image multi-photon absorption in heterogeneous media via photoacoustic techniques. A reader would care because the results remove a practical restriction on input size and allow the same reconstruction framework to apply across different physical scales.

Core claim

The absorption coefficient is a nonlinear function of the angular average of the transport solution; with this modeling choice the forward problem possesses a unique solution for general boundary data, while the inverse problem admits a stable reconstruction whose L1 error is controlled by a weighted norm that down-weights the boundary layer and thereby unifies the stability theory for both the pure transport and the diffusion limits.

What carries the argument

The weighted norm that penalizes the contribution from the boundary region, which converts separate L1 stability statements for the transport and diffusion regimes into a single uniform estimate.

If this is right

  • Larger boundary inputs become admissible in photoacoustic experiments without losing well-posedness.
  • A single reconstruction algorithm and error bound can be used for both ballistic and diffusive regimes.
  • The weighted norm supplies a concrete way to quantify how boundary truncation affects stability.
  • Multi-photon absorption maps can be recovered directly from measured internal fluence data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same weighted-norm technique may apply to other nonlinearities in transport models that are not strictly angular averages.
  • Numerical schemes that enforce the weighted norm could improve conditioning when boundary data are noisy or truncated.
  • The unified stability result suggests that hybrid diffusion-transport forward solvers can share the same inverse framework.

Load-bearing premise

The absorption coefficient is assumed to be a function of the angular average of the transport solution, and internal data for that average are available.

What would settle it

An explicit large-boundary-data example in which either existence or uniqueness of the forward solution fails, or a family of coefficients whose internal-data difference is small yet whose L1 difference is large even after weighting.

read the original abstract

We study forward and inverse problems for a semilinear radiative transport model where the absorption coefficient depends on the angular average of the transport solution. Our first result is the well-posedness theory for the transport model with general boundary data, which significantly improves previous theories for small boundary data. For the inverse problem of reconstructing the nonlinear absorption coefficient from internal data, we develop stability results for the reconstructions and unify an $L^1$ stability theory for both the diffusion and transport regimes by introducing a weighted norm that penalizes the contribution from the boundary region. The problems studied here are motivated by applications such as photoacoustic imaging of multi-photon absorption of heterogeneous media.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript develops well-posedness theory for a semilinear radiative transport equation in which the absorption coefficient depends on the angular average of the solution, extending prior results from small to general boundary data. It further establishes stability estimates for the inverse problem of recovering the nonlinear absorption coefficient from internal data, unifying an L¹ stability theory across diffusion and transport regimes through a weighted norm that penalizes boundary contributions. The work is motivated by multi-photon photoacoustic imaging applications.

Significance. If the central claims hold, the results advance the analysis of semilinear transport models by removing restrictive smallness conditions on boundary data and by providing a unified stability framework that bridges regimes. The introduction of the weighted norm for L¹ estimates is a technically useful device that could inform related inverse problems in imaging.

major comments (2)
  1. [§3] §3, Theorem 3.2: the fixed-point argument for global well-posedness with general boundary data requires a uniform bound on the Lipschitz constant of the nonlinearity; the proof sketch does not explicitly verify that this constant remains controlled independently of the L¹ norm of the boundary flux, which is the key improvement over prior small-data results.
  2. [§5] §5, Eq. (5.7): the weighted norm used to unify the L¹ stability estimates across regimes is defined by adding a boundary penalty term; it is not shown whether this norm is equivalent (with constants independent of the absorption coefficient) to the standard L¹ norm in the pure transport regime, which would be needed to confirm the unification claim.
minor comments (3)
  1. [§2] The statement of the semilinear model in §2 should explicitly record the precise functional setting (e.g., the space for the angular average) to make the dependence of absorption on the solution unambiguous.
  2. [Figure 1] Figure 1 (schematic of the imaging geometry) would benefit from a caption that distinguishes the internal data region from the boundary support.
  3. [§5] A short remark comparing the obtained stability constants with those in the diffusion approximation (e.g., in the cited works) would help readers assess the improvement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment of the significance, and recommendation for minor revision. We address each major comment below with clarifications and will incorporate explicit details into the revised manuscript to strengthen the presentation.

read point-by-point responses

  1. Referee: [§3] §3, Theorem 3.2: the fixed-point argument for global well-posedness with general boundary data requires a uniform bound on the Lipschitz constant of the nonlinearity; the proof sketch does not explicitly verify that this constant remains controlled independently of the L¹ norm of the boundary flux, which is the key improvement over prior small-data results.

    Authors: We appreciate the referee's observation on the need for explicit verification. In the proof of Theorem 3.2, the fixed-point map is constructed on a ball whose radius scales with the L¹ norm of the boundary flux, and the contraction relies on the Lipschitz constant of the nonlinearity being bounded uniformly. This bound follows from the global Lipschitz assumption on the nonlinearity combined with the a priori L¹ estimates obtained from the linear transport operator along characteristics, which are independent of the specific size of the boundary data. To address the lack of explicit verification in the sketch, we will expand the argument in the revised version by inserting the intermediate estimate that isolates the Lipschitz factor and shows its independence from the boundary L¹ norm using the integral form of the solution. revision: yes

  2. Referee: [§5] §5, Eq. (5.7): the weighted norm used to unify the L¹ stability estimates across regimes is defined by adding a boundary penalty term; it is not shown whether this norm is equivalent (with constants independent of the absorption coefficient) to the standard L¹ norm in the pure transport regime, which would be needed to confirm the unification claim.

    Authors: We thank the referee for raising this point regarding the unification claim. The weighted norm in Eq. (5.7) is designed to absorb boundary contributions while recovering the standard L¹ stability in the interior. In the pure transport regime, equivalence to the usual L¹ norm holds with constants depending only on the domain geometry, the transport speed, and the time horizon, but independent of the absorption coefficient; this follows from integrating along characteristics and controlling the boundary penalty term via the incoming flux data. We will add a short lemma in Section 5 establishing this equivalence explicitly, thereby confirming that the unified L¹ theory applies without additional restrictions on the absorption coefficient. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via standard analysis

full rationale

The paper establishes well-posedness for the semilinear transport model with general boundary data and L1 stability results unifying diffusion and transport regimes through a weighted norm. These rest on standard functional analysis techniques (e.g., fixed-point arguments or energy estimates) applied to the stated semilinear structure and internal data assumptions, which are introduced explicitly as modeling choices motivated by photoacoustic imaging rather than derived from the results themselves. No step reduces a claimed prediction or uniqueness result to a fitted input, self-citation chain, or definitional tautology by the paper's own equations. The improvement over prior small-data theories is presented as an extension of existing methods without load-bearing reliance on unverified self-references.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review prevents full enumeration; the model assumes a semilinear dependence of absorption on angular average and availability of internal data, which are domain assumptions rather than derived.

axioms (1)
  • domain assumption The semilinear radiative transport model is well-defined with absorption depending on angular average of the solution
    Stated directly in the abstract as the model under study

pith-pipeline@v0.9.0 · 5625 in / 1161 out tokens · 37857 ms · 2026-05-18T17:55:36.294947+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

  1. [1]

    Bal,Inverse transport theory and applications, Inverse Problems, 25 (2009)

    [1]G. Bal,Inverse transport theory and applications, Inverse Problems, 25 (2009). 053001. [2]G. Bal, A. Jollivet, and V. Jugnon,Inverse transport theory of photoacoustics, Inverse Problems, 26 (2010). 025011. [3]P. Bardsley, K. Ren, and R. Zhang,Quantitative photoacoustic imaging of two- photon absorption, J. Biomed. Opt., 23 (2018), p. 016002. [4]M. S. B...

    work page 2009

  2. [2]

    [5]P. Chen, I. M. Gamba, Q. Li, and L. Wang,Reconstruction of heat relaxation index in phonon transport equation, arXiv:2502.19533, (2025). [6]M. Choulli and P. Stefanov,Reconstruction of the coefficients of the stationary transport equation from boundary measurements, Inverse Probl., 12 (1996), pp. L19– L23. [7]F. J. Chung and J. C. Schotland,Inverse tra...

    work page arXiv 2025

  3. [3]

    [9]W. Denk, J. Strickler, and W. Webb,Two photon laser scanning fluorescence microscopy, Science, 248 (1990), pp. 73–76. [10]T. Ding, K. Ren, and S. Vall ´elian,A one-step reconstruction algorithm for quan- titative photoacoustic imaging, Inverse Problems, 31 (2015), p. 095005. [11]I. M. Gamba, Q. Li, and A. Nair,Reconstructing the thermal phonon transmis...

    work page 1990

  4. [4]

    Golse, P.-L

    21 [13]F. Golse, P.-L. Lions, B. Perthame, and R. Sentis,Regularity of the moments of the solution of a transport equation, Journal of functional analysis, 76 (1988), pp. 110–

    work page 1988

  5. [5]

    Hellmuth, C

    [14]K. Hellmuth, C. Klingenberg, Q. Li, and M. Tang,Kinetic chemotaxis tum- bling kernel determined from macroscopic quantities, SIAM Journal on Mathematical Analysis, 56 (2024), pp. 568–587. [15]R. B. Kellogg,Uniqueness in the Schauder fixed point theorem, Proceedings of the American Mathematical Society, 60 (1976), pp. 207–210. [16]M. V. Klibanov, J. Li...

    work page 2024

  6. [6]

    Lai and L

    [22]R.-Y. Lai and L. Yan,Stable determination of time-dependent collision kernel in the nonlinear boltzmann equation, SIAM Journal on Applied Mathematics, 84 (2024), pp. 1937–1956. [23]R.-Y. Lai and H. Zhou,Inverse problems for time-dependent nonlinear transport equations, arXiv preprint arXiv:2410.00369, (2024). [24]L. Li and Z. Ouyang,Determining the co...

    work page arXiv 2024

  7. [7]

    Machida and M

    [27]M. Machida and M. Yamamoto,Global lipschitz stability in determining coefficients of the radiative transport equation, Inverse Problems, 30 (2014), p. 035010. [28]H. Mahr,Two-photon absorption spectroscopy, in Quantum Electronics: A Treatise, H. Rabin and C. L. Tang, eds., Academic Press,

    work page 2014

  8. [8]

    [29]A. V. Mamonov and K. Ren,Quantitative photoacoustic imaging in radiative trans- port regime, arXiv preprint arXiv:1207.4664, (2012). [30]P. P. Mondal, G. Vicidomini, and A. Diaspro,Image reconstruction for multi- photon fluorescence microscopy, Appl. Phys. Lett., 92 (2008). 103902. [31]M. Reed and B. Simon,IV: Analysis of Operators, vol. 4, Elsevier,

    work page arXiv 2012

  9. [9]

    [32]K. Ren, H. Gao, and H. Zhao,A hybrid reconstruction method for quantitative PAT, SIAM Journal on Imaging Sciences, 6 (2013), pp. 32–55. [33]K. Ren and R. Zhang,Nonlinear quantitative photoacoustic tomography with two- photon absorption, SIAM J. Appl. Math., 78 (2018), pp. 479–503. [34]K. Ren, R. Zhang, and Y. Zhong,A fast algorithm for radiative trans...

    work page 2013

  10. [10]

    Sun and Q

    [41]W. Sun and Q. Li,Unique reconstruction of the heat-reflection indices at solid inter- faces, SIAM Journal on Mathematical Analysis, 54 (2022), pp. 5742–5759. [42]M. E. Taylor et al.,Partial differential equations III, vol. 2, Springer,

    work page 2022

  11. [11]

    Wang,Stability estimates of an inverse problem for the stationary transport equation, Ann

    23 [43]J.-N. Wang,Stability estimates of an inverse problem for the stationary transport equation, Ann. Inst. Henri Poincar´ e, 70 (1999), pp. 473–495. [44]D. Xu, X. Wang, Z. Xu, W. Zhou, and J. Yin,A super-resolution reconstruction algorithm for two-photon fluorescence polarization microscopy, Optics Communications, 499 (2021). 127116. [45]Y. Yamaoka, M....

    work page 1999

  12. [12]

    [47]J. Ying, F. Liu, and R. R. Alfano,Spatial distribution of two-photon-excited fluorescence in scattering media, Appl. Opt., 38 (1999), pp. 224–229. [48]H. Zhao and Y. Zhong,Instability of an inverse problem for the stationary radiative transport near the diffusion limit, SIAM Journal on Mathematical Analysis, 51 (2019), pp. 3750–3768. [49]H. Zhao and Y...

    work page 1999