REVIEW 2 major objections 2 minor 1 cited by
Inverse initial data reconstruction for a memory convection-diffusion equation via Legendre spatial reduction and Tikhonov regularization
T0 review · 2 major / 2 minor · reviewed 2026-06-26 · grok-4.3
Pith's one-line read Legendre spatial reduction combined with Tikhonov regularization recovers the initial condition from final-time data for a memory convection-diffusion equation.
desk verdict Standard Legendre reduction plus Tikhonov for a memory PDE inverse problem; convergence proof holds only under an external parameter choice. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Finite tensor-product Legendre expansion in space that reduces the inverse problem to a terminal-value system for time-dependent coefficients, solved by Tikhonov-regularized least squares with an H^2 penalty.
What would settle it
A concrete sequence of noisy final-time data and corresponding regularization parameters, both tending to zero, for which the regularized Legendre-coefficient vectors fail to approach the minimum-norm solution of the reduced system at a fixed truncation order.
Extended reading notes
Core claim
For a fixed truncation order, the Tikhonov-regularized minimizers converge to the finite-dimensional minimum-norm solution as the noise level and the regularization parameter vanish, under a suitable choice of the regularization parameter.
Load-bearing premise
A regularization parameter can be chosen that vanishes with the noise level in a way that does not require knowledge of the exact solution.
Editorial extensions
If this is right
- Uniqueness of the initial datum holds when coefficients are independent of space, proved via Fourier transform and analytic continuation of a scalar Volterra equation.
- The original infinite-dimensional ill-posed problem is replaced by a finite-dimensional, computable least-squares problem whose regularized solutions converge under the stated parameter rule.
- Two-dimensional numerical tests confirm that the reduced regularized reconstructions remain stable for moderate noise levels.
- The memory term is retained exactly inside the reduced system of ordinary differential equations for the Legendre coefficients.
Reading between the lines
- Separate analysis of the truncation error is needed before the method can be applied to data with arbitrarily fine spatial features.
- The same reduction-plus-regularization pattern may apply directly to other linear parabolic equations that contain nonlocal time operators.
- Practical implementation would still require an a-posteriori rule such as the discrepancy principle to select the regularization parameter from data alone.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies an inverse initial-data problem for a convection-diffusion equation with memory term. For spatially constant coefficients it proves uniqueness via Fourier transform reducing the problem to a scalar Volterra equation whose analyticity yields uniqueness. For variable coefficients it introduces a Legendre-Galerkin spatial reduction that converts the PDE inverse problem into a finite-dimensional terminal-value ODE system for the coefficient vector; this reduced problem is solved by Tikhonov regularization with an H² penalty. The central theoretical result states that, for any fixed truncation order N, the regularized minimizers converge to the minimum-norm solution of the finite-dimensional problem as the noise level δ and the regularization parameter α both tend to zero, provided α is chosen suitably (i.e., α→0 and δ²/α→0). Two-dimensional numerical illustrations are presented.
Significance. If the convergence statement can be made constructive, the combination of exact spatial reduction with a provably convergent regularization scheme would supply a theoretically supported numerical method for a class of severely ill-posed inverse problems that include memory. The explicit reduction to a finite terminal-value system and the convergence proof for that reduced system are genuine strengths; the numerical examples provide at least preliminary evidence of practical behavior.
major comments (2)
- [Abstract / convergence theorem] Abstract and convergence statement: the claim that the regularized minimizers converge to the finite-dimensional minimum-norm solution 'under a suitable choice of the regularization parameter' is load-bearing for the method, yet the manuscript gives no indication that an a-posteriori rule (discrepancy principle, balancing principle, etc.) is derived or that the proof supplies a constructive selection of α(δ) independent of the unknown solution. The standard source-condition argument requires α→0 and δ²/α→0, both of which remain external to the data.
- [Legendre spatial reduction section] Reduction step: the passage from the original memory PDE to the finite tensor-product Legendre system for the time-dependent coefficients must preserve the memory integral exactly; any truncation or quadrature error introduced at this stage would propagate into the subsequent Tikhonov analysis. The manuscript should state explicitly whether the memory term is discretized exactly or approximated and, if the latter, quantify the additional consistency error.
minor comments (2)
- [Abstract] The abstract states that uniqueness holds 'in a spatially independent coefficient setting' but does not clarify whether the same Fourier-Volterra argument extends, even formally, to the variable-coefficient case before the numerical reduction is applied.
- [Notation / preliminaries] Notation for the memory kernel and the precise form of the H² penalty should be introduced once and used consistently; several symbols appear only in the abstract and are not redefined in the main text.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for highlighting both the strengths of the spatial reduction and the limitations of the convergence result. We respond to each major comment below.
read point-by-point responses
-
Referee: [Abstract / convergence theorem] Abstract and convergence statement: the claim that the regularized minimizers converge to the finite-dimensional minimum-norm solution 'under a suitable choice of the regularization parameter' is load-bearing for the method, yet the manuscript gives no indication that an a-posteriori rule (discrepancy principle, balancing principle, etc.) is derived or that the proof supplies a constructive selection of α(δ) independent of the unknown solution. The standard source-condition argument requires α→0 and δ²/α→0, both of which remain external to the data.
Authors: We agree that the convergence theorem is stated under the external conditions α→0 and δ²/α→0 and that no a-posteriori selection rule is derived or analyzed in the manuscript. The result therefore remains non-constructive with respect to practical parameter choice. In the revision we will modify the abstract and the theorem statement to make this assumption explicit and will add a short remark in the conclusions noting that discrepancy-principle or balancing-principle rules could be considered in applications, while acknowledging that their rigorous justification would require additional source conditions not treated here. revision: partial
-
Referee: [Legendre spatial reduction section] Reduction step: the passage from the original memory PDE to the finite tensor-product Legendre system for the time-dependent coefficients must preserve the memory integral exactly; any truncation or quadrature error introduced at this stage would propagate into the subsequent Tikhonov analysis. The manuscript should state explicitly whether the memory term is discretized exactly or approximated and, if the latter, quantify the additional consistency error.
Authors: The reduction is obtained by applying the orthogonal projection onto the tensor-product Legendre space directly to the entire PDE, including the memory integral. Because the projection is linear, the memory term is transferred exactly into the finite-dimensional system; no truncation or quadrature approximation is introduced at the reduction stage. The resulting terminal-value ODE system is therefore the exact Galerkin projection of the original problem. We will insert an explicit clarifying paragraph in the Legendre spatial reduction section stating that the memory integral is preserved exactly and that no additional consistency error arises from the reduction itself. revision: yes
Circularity Check
No significant circularity; derivation is self-contained
full rationale
Uniqueness is established via external Fourier transform plus analyticity of a scalar Volterra equation. The convergence statement is a standard Tikhonov result proved directly on the explicitly constructed finite-dimensional reduced system after Legendre truncation; it does not reduce by the paper's own equations to any fitted quantity, self-citation chain, or ansatz smuggled from prior work. The a-priori parameter choice is external and does not create a definitional loop. No load-bearing self-citations or renamings appear in the derivation chain.
Assumptions & free parameters
free parameters (2)
- truncation order
- regularization parameter
assumptions (2)
- domain assumption Fourier transform converts the PDE into a scalar Volterra equation whose solution is analytic
- standard math Legendre polynomials form a complete orthogonal basis allowing exact spatial reduction for the chosen truncation
Cite this review
Pith. "Pith review of Inverse initial data reconstruction for a memory convection-diffusion equation via Legendre spatial reduction and Tikhonov regularization." pith.science (2026). https://pith.science/paper/FEB3PUAL
@misc{pith2026260620875,
author = {Pith},
title = {Pith review of: Inverse initial data reconstruction for a memory convection-diffusion equation via Legendre spatial reduction and Tikhonov regularization},
year = {2026},
howpublished = {\url{https://pith.science/paper/FEB3PUAL}},
note = {Machine review of arXiv:2606.20875}
}
abstract
We study an inverse initial data problem for a convection-diffusion equation with memory, where the goal is to recover the unknown initial condition from final-time data. The model includes convection, an instantaneous Laplacian term, and a nonlocal-in-time memory term involving the Laplacian of the past states, which leads to a severely ill-posed backward problem. We prove uniqueness in a spatially independent coefficient setting by applying the Fourier transform and using an analyticity argument for a scalar Volterra equation. For the variable-coefficient case, we develop a computational method based on Legendre spatial dimensional reduction and Tikhonov regularization. The solution is approximated by a finite tensor-product Legendre expansion, thereby reducing the inverse problem to a finite-dimensional terminal-value system for the time-dependent coefficients. We solve the reduced problem by a Tikhonov-regularized least-squares method with an $H^2$ penalty. For a fixed truncation order, we prove that the regularized minimizers converge to the finite-dimensional minimum-norm solution as the noise level and the regularization parameter vanish, under a suitable choice of the regularization parameter. Some two-dimensional numerical examples are presented to illustrate the performance of the proposed method.
Figures
Forward citations
Cited by 1 Pith paper
-
Carleman--Picard and time-dimensional reduction for inverse initial-data problems in nonlinear transport with memory
A Carleman-Picard iteration with Legendre-exponential time reduction globally converges, within a truncated reduced model, for reconstructing initial data of quasilinear transport with memory from outflow measurements.
Reference graph
Works this paper leans on
-
[1]
K. A. Ames, G. W. Clark, J. F. Epperson, and S. F. Oppenheimer. A comparison of regular- izations for an ill-posed problem.Mathematics of Computation, 67(224):1451–1471, 1998
1998
-
[2]
K. A. Ames and J. F. Epperson. A kernel-based method for the approximate solution of backward parabolic problems.SIAM Journal on Numerical Analysis, 34(4):1357–1390, 1997
1997
-
[3]
Atifi, E.-H
K. Atifi, E.-H. Essoufi, and B. Khouiti. An inverse backward problem for degenerate two- dimensional parabolic equation.Opuscula Mathematica, 40(4):427–448, 2020
2020
-
[4]
Canuto and A
C. Canuto and A. Quarteroni. Approximation results for orthogonal polynomials in Sobolev spaces.Mathematics of Computation, 38(157):67–86, 1982
1982
-
[5]
A. S. Carasso. Determining surface temperature from interior observations.SIAM Journal on Applied Mathematics, 42(3):558–574, 1982
1982
-
[6]
Chen and J
Q. Chen and J. Liu. Solving an inverse parabolic problem by optimization from final measure- ment data.Journal of Computational and Applied Mathematics, 193:183–203, 2006
2006
-
[7]
D. K. Durdiev and A. Sh. Rashidov. Inverse problem of determining the kernel in an integro- differential equation of parabolic type.Differential Equations, 50(1):110–116, 2014. 26
2014
-
[8]
D. K. Durdiev and Zh. Zh. Zhumaev. Memory kernel reconstruction problems in the integro- differential equation of rigid heat conductor.Mathematical Methods in the Applied Sciences, 45(14):8374–8388, 2022
2022
Show all 53 references
-
[9]
D. K. Durdiev and Zh. Zh. Zhumaev. One-dimensional inverse problems of finding the kernel of integro-differential heat equation in a bounded domain.Ukrainian Mathematical Journal, 73(3):1723–1740, 2022
2022
-
[10]
D. K. Durdiev and Zh. Zh. Zhumaev. On determination of the coefficient and kernel in an integro-differential equation of parabolic type.Eurasian Journal of Mathematical and Com- puter Applications, 11(1):49–65, 2023
2023
-
[11]
Efendiev, W
Y. Efendiev, W. T. Leung, W. Li, S.-M. Pun, and P. N. Vabishchevich. Nonlocal transport equations in multiscale media: Modeling, dememorization, and discretizations.Journal of Computational Physics, 472:111555, 2023
2023
-
[12]
R. E. Ewing. The approximation of certain parabolic equations backward in time by Sobolev equations.SIAM Journal on Mathematical Analysis, 6(2):283–294, 1975
1975
-
[13]
C.-L. Fu. Simplified Tikhonov and Fourier regularization methods on a general sideways parabolic equation.Journal of Computational and Applied Mathematics, 167(2):449–463, 2004
2004
-
[14]
H. Han, D. Ingham, and Y. Yuan. The boundary element method for the solution of the backward heat conduction equation.Journal of Computational Physics, 116:292–299, 1995
1995
-
[15]
P. C. Hansen and D. P. O’Leary. The use of the L-curve in the regularization of discrete ill-posed problems.SIAM Journal on Scientific Computing, 14:1487–1503, 1993
1993
-
[16]
D. N. Hao, N. V. Duc, N. V. Thang, and T. Nguyen. Regularization of backward time- fractional parabolic equations by Sobolev-type equations.Journal of Inverse and Ill-Posed Problems, 28(5):659–676, 2020
2020
-
[17]
D. N. Hao, V. N. Duc, and D. Lesnic. Regularization of parabolic equations backward in time by a non-local boundary value problem method.IMA Journal of Applied Mathematics, 75:291–315, 2010
2010
-
[18]
A. H. Hasanov and J. L. Mueller. A numerical method for backward parabolic problems with non-selfadjoint elliptic operators.Applied Numerical Mathematics, 37:55–78, 2001
2001
-
[19]
Isakov.Inverse Problems for Partial Differential Equations, volume 127 ofApplied Mathe- matical Sciences
V. Isakov.Inverse Problems for Partial Differential Equations, volume 127 ofApplied Mathe- matical Sciences. Springer, New York, 2 edition, 2006
2006
-
[20]
B. Wu, Y. Gao, L. Yan, S. Wu, and Z. Wang. Existence and uniqueness of an inverse memory kernel for an integro-differential parabolic equation with free boundary.Journal of Dynamical and Control Systems, 24:237–252, 2018
2018
-
[21]
Janno and A
J. Janno and A. Lorenzi. Recovering memory kernels in parabolic transmission problems. Journal of Inverse and Ill-Posed Problems, 16:239–265, 2008
2008
-
[22]
M. V. Klibanov. Carleman estimates for the regularization of ill-posed Cauchy problems. Applied Numerical Mathematics, 94:46–74, 2015
2015
-
[23]
M. V. Klibanov and A. G. Yagola. Convergent numerical methods for parabolic equations with reversed time via a new Carleman estimate.Inverse Problems, 35(11):115012, 2019. 27
2019
-
[24]
M. N. Koleva and L. G. Vulkov. Numerical solution of the retrospective inverse parabolic problem on disjoint intervals.Computation, 11(10):Article 204, 2023
2023
-
[25]
Langer, O
U. Langer, O. Steinbach, F. Tr¨ oltzsch, and H. Yang. Space-Time Finite Element Methods for Initial Temperature Reconstruction. InDomain Decomposition Methods in Science and Engineering XXVI, volume 145 ofLecture Notes in Computational Science and Engineering, pages 617–624. S...
2023
-
[26]
Latt` es and J.-L
R. Latt` es and J.-L. Lions.M´ ethode de quasi-r´ eversibilit´ e et Applications. Dunod, Paris, 1967
1967
-
[27]
M. M. Lavrent’ev, V. G. Romanov, and S. P. Shishatski˘ ı.Ill-Posed Problems of Mathematical Physics and Analysis, volume 64 ofTranslations of Mathematical Monographs. American Mathematical Society, Providence, RI, 1986
1986
-
[28]
T. T. Le. Global reconstruction of initial conditions of nonlinear parabolic equations via the Carleman-contraction method. In D.-L. Nguyen, L. H. Nguyen, and T.-P. Nguyen, editors, Advances in Inverse problems for Partial Differential Equations, volume 784 ofContemporary Math...
2023
-
[29]
T. T. Le and L. H. Nguyen. A convergent numerical method to recover the initial condition of nonlinear parabolic equations from lateral Cauchy data.Journal of Inverse and Ill-posed Problems, 30(2):265–286, 2022
2022
-
[30]
T. T. Le, C. B. Van, D. T. Dang, and L. H. Nguyen. Inverse initial data reconstruction for Maxwell’s equations via time-dimensional reduction method.Journal of Computational Physics, 559:114896, 2026
2026
-
[31]
Thuy T Le, Minh-Binh Tran, and Loc H Nguyen. A globally convergent Carleman-Picard method for an inverse initial-value problem for a nonlinear diffusive coagulation-fragmentation equation coagulation-fragmentation equation.arXiv preprint arXiv:2603.21185, 2026
2026
-
[32]
Li and L
Q. Li and L. H. Nguyen. Recovering the initial condition of parabolic equations from lateral Cauchy data via the quasi-reversibility method.Inverse Problems in Science and Engineering, 2019
2019
-
[33]
T. N. Luan and T. Q. Khanh. On the backward problem for parabolic equations with memory. Applicable Analysis, 100(7):1414–1431, 2021
2021
-
[34]
N. S. Mera. The method of fundamental solutions for the backward heat conduction problem. Inverse Problems in Science and Engineering, 13(1):65–78, 2005
2005
-
[35]
K. Miller. Stabilized quasi-reversibility and other nearly-best-possible methods for non-well- posed problems. InSymposium on Non-Well-Posed Problems and Logarithmic Convexity, volume 316 ofLecture Notes in Mathematics, pages 161–176. Springer-Verlag, 1973
1973
-
[36]
Inverse initial data for nonlinear Schr¨ odinger equation via Carleman estimates and the contraction principle.arXiv preprint arXiv:2605.11409, 2026
Navaraj Neupane and Loc Nguyen. Inverse initial data for nonlinear Schr¨ odinger equation via Carleman estimates and the contraction principle.arXiv preprint arXiv:2605.11409, 2026
2026 arXiv
-
[37]
P. M. Nguyen, M. Nguyen, and L. H. Nguyen. Forward-time Black–Scholes reconstruction via regularized Legendre reduction.preprint arxiv:2606.12450, 2026
2026 arXiv
-
[38]
A Carleman contraction method for inverse initial data recovery in the Navier-Stokes equations with unknown body force.arXiv e-prints, pages arXiv–2604, 2026
Phuong M Nguyen and Loc H Nguyen. A Carleman contraction method for inverse initial data recovery in the Navier-Stokes equations with unknown body force.arXiv e-prints, pages arXiv–2604, 2026. 28
2026
-
[39]
L. E. Payne.Improperly Posed Problems in Partial Differential Equations. SIAM, Philadelphia, 1975
1975
-
[40]
Prilepko, D
A. Prilepko, D. Orlovsky, and I. Vasin.Methods for Solving Inverse Problems in Mathematical Physics. Marcel Dekker, New York, 2000
2000
-
[41]
Pham Hoang Quan, Le Minh Triet, and Dang Duc Trong. A regularization of the backward problem for nonlinear parabolic equation with Time-dependent coefficient.International Jour- nal of Mathematics and Mathematical Sciences, 2012:Article 943621, 2012
2012
-
[42]
A. A. Samarskii, P. N. Vabishchevich, and V. I. Vasil’ev. Iterative solution of a retrospec- tive inverse problem of heat conduction.Mathematical Models and Computer Simulations / Matematicheskoe Modelirovanie, 9:119–127, 1997
1997
-
[43]
R. E. Showalter. The final value problem for evolution equations.Journal of Mathematical Analysis and Applications, 47:563–572, 1974
1974
-
[44]
L. Su, J. Huang, V. I. Vasil’ev, A. Li, and A. M. Kardashevsky. A numerical method for solving retrospective inverse problem of fractional parabolic equation.Journal of Computational and Applied Mathematics, 413:Article 114366, 2022
2022
-
[45]
A. N. Tikhonov and V. Y. Arsenin.Solutions of Ill-Posed Problems. Winston & Sons, Wash- ington, D.C., 1977
1977
-
[46]
Trong and N
D. Trong and N. H. Tuan. Regularization and error estimates for nonhomogeneous backward heat problems.Electronic Journal of Differential Equations, pages 1–10, 2006
2006
-
[47]
P. H. Quan, D. D. Trong, L. M. Triet, and N. H. Tuan. A modified quasi-boundary value method for regularizing a backward problem with time-dependent coefficient.Inverse Problems in Science and Engineering, 19(3):409–423, 2011
2011
-
[48]
D. D. Trong, C. V. Le, K. D. Luu, and L. H. Nguyen. Recovery of initial displacement and velocity in anisotropic elastic systems by the time dimensional reduction method.Journal of Computational Physics, 542:114371, 2025
2025
-
[49]
C. B. Van, T. T. Le, and L. H. Nguyen. The inverse initial data problem for anisotropic Navier–Stokes equations via Legendre time reduction method.Communications in Nonlinear Science and Numerical Simulation,, 161:110074, 2026
2026
-
[50]
Yang and Z.-C
L. Yang and Z.-C. Deng. An inverse backward problem for degenerate parabolic equations. Numerical Methods for Partial Differential Equations, 33(6):1900–1923, 2017
1900
-
[51]
Cheng and C
W. Cheng and C. L. Fu. A modified Tikhonov regularization method for an axisymmetric backward heat equation.Acta Mathematica Sinica, English Series, 26(11):2157–2164, 2010
2010
-
[52]
W.-Q. Hou, F. Yang, X.-X. Li, and Z.-J. Tian. The modified Tikhonov regularization method for backward heat conduction problem with a complete parabolic equation inR n.Numerical Algorithms, 102:981–1010, 2026
2026
-
[53]
Zhang and Z
Z. Zhang and Z. Zhou. Backward diffusion-wave problem: stability, regularization and approx- imation.SIAM Journal on Scientific Computing, 44(5):A3183–A3216, 2022. 29
2022
Reviewed June 26, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.