---
title: Learnable Fractional Reaction-Diffusion (LFRD2)
url: https://www.emergentmind.com/topics/learnable-fractional-reaction-diffusion-dynamics-lfrd2
type: topic
---

# Learnable Fractional Reaction-Diffusion (LFRD2)

Searching arXiv for LFRD2 and closely related fractional reaction–diffusion work.
Learnable Fractional Reaction-Diffusion Dynamics (LFRD2) denotes a hybrid, physics-informed framework that embeds a time-fractional reaction-diffusion evolution law inside a learnable depth-restoration pipeline. In its defining formulation, introduced for under-display Time-of-Flight imaging, LFRD2 combines a deep initial estimator with a Deep Fractional Reaction-Diffusion module that performs iterative refinement through a Caputo time-fractional reaction-diffusion update with learnable differential orders and an efficient continuous convolution operator [2511.01704]. More broadly, LFRD2 sits at the intersection of fractional PDE modeling, operator-based image restoration, and learnable nonlocal dynamics: it inherits long-memory temporal behavior from fractional calculus, nonlocal transport structure from fractional diffusion theory, and data adaptivity from neural parameterization [2511.01704; 1502.06127; 1211.0063].

## 1. Origin and defining formulation

The term LFRD2 originates in “Learnable Fractional Reaction-Diffusion Dynamics for Under-Display ToF Imaging and Beyond” [2511.01704]. The framework was introduced to address under-display ToF degradation caused by transparent OLED layers, specifically signal attenuation, multi-path interference, and temporal noise, which significantly compromise depth quality [2511.01704]. The proposed remedy is a two-stage architecture composed of a Deep Initial State Builder and a Deep Fractional Reaction-Diffusion module [2511.01704].

In this formulation, the initial depth estimate \(u_0\) is produced by a backbone network, and refinement proceeds iteratively through
\[
u_{n+1}
= \sum_{t=0}^n w_t u_t + \mathcal{D}_t(u_n) + \mathcal{R}_t(u_n, u_0),
\]
where \(\mathcal{D}_t\) is a diffusion term, \(\mathcal{R}_t\) is a reaction term, and \(w_t\) are memory weights induced by the fractional derivative [2511.01704]. The underlying PDE is written as
\[
^C_0 D_t^{\alpha} u_{n+1} = \mathrm{div}\bigl(g(|\nabla u_n|)\nabla u_n\bigr) + \lambda (u_0 - u_n),
\]
with \( ^C_0 D_t^{\alpha}\) the Caputo time-fractional derivative, \(g(\cdot)\) a learnable diffusivity, and \(\lambda=0.01\) following TNRD [2511.01704].

This structure places LFRD2 in a longer lineage of fractional reaction-diffusion models. Earlier work studied computable solutions for fractional reaction-diffusion equations with generalized Riemann-Liouville time derivatives and Riesz-Feller space derivatives, using Laplace-Fourier methods and generalized Mittag-Leffler functions [1502.06127]. Related work also derived closed-form solutions for distributed-order systems with Riemann-Liouville time derivatives and Riesz-Feller space operators, with solutions expressed through Mittag-Leffler and Srivastava-Daoust functions [1211.0063]. LFRD2 differs in that it does not primarily seek closed-form PDE solutions; instead, it makes the PDE dynamics learnable and task-adaptive [2511.01704].

## 2. Fractional temporal memory and reaction-diffusion structure

The central mathematical distinction of LFRD2 is the replacement of the standard first-order time derivative by a Caputo derivative of order \(\alpha\in(0,1)\), defined as
\[
^C_0D_t^{\alpha} u(t) = \frac{1}{\Gamma(1-\alpha)} \int_0^t (t-\tau)^{-\alpha} u'(\tau)\,d\tau.
\]
This kernel weights all past times, so the refinement at iteration \(n+1\) depends on the entire history of previous states rather than only on \(u_n\) [2511.01704]. The paper explicitly contrasts this with integer-order dynamics, where the update is Markovian [2511.01704].

The diffusion term in LFRD2 has Perona-Malik form,
\[
\mathrm{div}(g(|\nabla u_n|)\nabla u_n),
\]
while the reaction term is
\[
\lambda(u_0-u_n).
\]
The reaction component anchors the evolution to the initial estimate, thereby counteracting excessive smoothing, whereas the diffusion component performs nonlinear spatial regularization [2511.01704]. This split is consistent with broader fractional reaction-diffusion literature, where diffusion is represented by nonlocal or fractional operators and reaction contributes local decay, growth, or forcing terms [1502.06127; 1904.07651; 1806.07980].

A key conceptual point is that fractional time derivatives encode memory through power-law kernels rather than exponential forgetting. Earlier analytical studies showed that time-fractional reaction and diffusion equations are naturally expressed in terms of generalized Mittag-Leffler functions, which generalize exponential relaxation and yield heavy-tailed temporal behavior [1001.2293; 1502.06127; 1211.0063]. This suggests that LFRD2’s use of a learnable \(\alpha\) is not merely a numerical variation on iterative refinement; it situates the model within a class of non-Markovian evolution laws whose effective memory spectrum can be adapted from data [2511.01704].

## 3. Discretization, learnability, and continuous convolution

LFRD2 discretizes the Caputo derivative with the L1 approximation
\[
^C_0 D_t^{\alpha} u_{n+1} \approx (\Delta t)^{-\alpha} \Gamma(2-\alpha) \sum_{k=0}^{n} a_k^{(\alpha)} \bigl(u_{n+1-k} - u_{n-k}\bigr),
\]
where
\[
a_k^{(\alpha)} = (k+1)^{1-\alpha}-k^{1-\alpha}.
\]
With \(\Delta t=1\), this leads to the explicit update
\[
\begin{aligned}
u_{n+1} &= u_n + S\bigl[ \mathrm{div}(g(|\nabla u_n|)\nabla u_n) + \lambda (u_0 - u_n) \bigr] \\
&\quad - \Bigl[\sum_{k=1}^{n} a_k^{(\alpha)}\bigl(u_{n+1-k}-u_{n-k}\bigr)\Bigr],
\end{aligned}
\]
with
\[
S = \frac{\Gamma(2-\alpha)}{a_0^{(\alpha)}}.
\]
The second line is the fractional history correction, and it is the mechanism by which all previous increments influence the current refinement [2511.01704].

The differentiability order \(\alpha\) is not fixed. The paper states that the differentiation orders are dynamically generated by a neural network, and ablations comparing fixed \(\alpha\in\{0.1,0.3,0.5,0.7,0.9\}\) against the variable-order model show that the variable-order model performs best [2511.01704]. This is the defining “learnable” aspect of LFRD2: the model does not merely tune coefficients in a preset PDE but adapts the effective temporal memory law itself [2511.01704].

The spatial diffusion operator is implemented through a continuous convolution construction rather than a standard finite-difference stencil. Building on a representation of convolution via repeated integration and differentiation,
\[
u * \mathcal{K}
=
u^{(-n)} * \mathcal{K}^{(n)},
\]
LFRD2 predicts estimated Dirac-delta-like kernel quantities directly from network features and approximates the repeated antiderivative by a local linear form \(u^{(-2)} \approx A u(x_0,y_0) + B\) [2511.01704]. The resulting operator is designed to be more efficient than Neural Field Convolutions while preserving a continuous-operator interpretation [2511.01704].

This learnable-operator view has clear precedents in the fractional PDE literature. Spectral formulations of the fractional Laplacian show that in Fourier space the operator acts by multiplication with \((k^2+l^2)^{\alpha/2}\) on periodic domains [1904.07651]. Riesz and Riesz-Feller derivatives similarly admit explicit Fourier symbols such as \(-|k|^\alpha\) in the symmetric case or \(-\psi_{\alpha,\theta}(k)\) in the skew case [1502.06127; 1211.0063]. LFRD2 does not explicitly parameterize these classical symbols, but it occupies the same operator-theoretic niche: a learnable approximation to nonlocal smoothing constrained by PDE form.

## 4. Relation to broader fractional reaction-diffusion theory

Fractional reaction-diffusion systems predate LFRD2 by many years. Earlier analytical work established explicit solution formulas for time-fractional and space-fractional equations using Laplace and Fourier transforms, with generalized Mittag-Leffler and Fox \(H\)-function kernels [1001.2293; 1502.06127; 1211.0063]. These studies make precise how fractional orders control memory, anomalous transport, and relaxation tails.

Distributed-order and multi-term time operators are especially relevant. One line of work considered equations of the form
\[
D_t^{\gamma_1,\delta_1}N + a D_t^{\gamma_2,\delta_2}N
=
\eta\,{}_xD_{\theta}^{\alpha}N - \omega N + U,
\]
with generalized Riemann-Liouville time derivatives and Riesz-Feller space derivatives, interpreting the two-term time operator as a discrete distributed-order model [1502.06127]. Another line studied
\[
{}_0D_t^{\alpha} N + a\, {}_0D_t^{\beta} N
=
\sum_{j=1}^{n} \mu_j\, {}_{x}D_{\theta_j}^{\gamma_j} N + \phi,
\]
and derived closed forms through Prabhakar and Srivastava-Daoust functions [1211.0063]. These results show that learnable fractional dynamics can, in principle, target not only a single time order but a mixture or distribution of temporal exponents.

The stochastic derivation of fractional transport from continuous-time random walks further clarifies the modeling role of these operators. Heavy-tailed waiting times lead to Caputo time derivatives, while heavy-tailed jumps lead to Riesz space derivatives, yielding equations such as
\[
\partial_t^\beta u(x,t)=\partial_{|x|}^\alpha u(x,t),
\]
or fractional reaction-diffusion equations where the reaction term itself is convolved with the same memory kernel [2503.13096]. This provides a mechanistic interpretation of why a learnable \(\alpha\) can be meaningful: it may encode effective waiting-time or transport heterogeneity rather than serving only as a generic tuning parameter [2503.13096].

On bounded domains, regularity theory shows that solution smoothness is limited by endpoint singular behavior determined by the fractional operator parameters [1911.03261]. This suggests a practical implication for learnable models: naive Euclidean smoothness assumptions may be mismatched to the actual function spaces of fractional PDE solutions. A plausible implication is that LFRD2-style models could benefit from architecture or loss designs that better reflect weighted Sobolev or nonlocal regularity structures, especially outside image-restoration settings [1911.03261].

## 5. Stability, well-posedness, and numerical dynamics

Fractional reaction-diffusion systems require a distinct stability theory because temporal decay is governed by sectorial spectral conditions rather than only by \(\Re\lambda<0\). Recent work proved a linearization principle for abstract fractional reaction-diffusion equations with Caputo derivative of order \(\alpha\in(0,1)\), showing that stability is associated with the spectral condition \(|\arg\lambda|>\alpha\pi/2\), while instability arises when some eigenvalue satisfies \(|\arg\lambda|<\alpha\pi/2\) [2507.02094]. The same work established a counterpart of classical Turing instability for fractional systems [2507.02094]. This provides a rigorous backdrop for any learnable fractional dynamics whose iterates approximate a time-fractional PDE.

Well-posedness on bounded domains with spectral fractional Laplacians has also been analyzed using analytic semigroup theory. For semilinear problems of the form
\[
\partial_t u + \varepsilon^2(-\Delta)_g^\alpha u = \mathcal{N}(u),
\]
one can reformulate inhomogeneous boundary conditions via harmonic lifting, derive Duhamel formulas, obtain local wellposedness for locally Lipschitz reactions, and establish positivity preservation and invariant sets for specific systems such as bistable equations and Gray-Scott dynamics [2601.21422]. This is relevant because it shows that learnable reaction terms can, at least in principle, be constrained so that the overall fractional dynamics remain well posed and positivity preserving [2601.21422].

Numerical studies further clarify what fractional order changes in practice. Stabilized semi-implicit Fourier spectral schemes for nonlinear space-fractional reaction-diffusion systems have shown that varying \(\alpha\) changes morphology in fractional Allen-Cahn, Gray-Scott, and FitzHugh-Nagumo models, and that the effect of changing \(\alpha\) is not equivalent to simply rescaling the diffusion coefficient [1904.07651]. Related Gray-Scott simulations on bounded domains showed that decreasing \(\alpha\) yields different steady patterns and a scaling law for radial distribution functions in terms of fractional order [1806.07980]. These results directly rebut a common oversimplification: fractional order is not merely an alternative parameterization of diffusion strength. It changes the operator class and therefore the spatiotemporal dynamics themselves [1904.07651; 1806.07980].

LFRD2’s empirical design choices are consistent with this broader picture. The model uses an explicit fractional update with learned order, reaction anchoring, and continuous convolution, and the paper reports that training can become unstable with poor strategies and that future work may consider implicit schemes [2511.01704]. This suggests that stability-aware numerical design remains an open issue for learnable fractional PDE modules, even when their task-level performance is strong.

## 6. Applications, empirical performance, and scope

LFRD2 was introduced for under-display ToF imaging, but its empirical scope is broader. On four benchmark datasets, the paper reports effectiveness for under-display ToF restoration, generic ToF denoising, and depth super-resolution [2511.01704]. On SUD-ToF, LFRD2 achieved MAE \(8.41\) mm and RMSE \(10.99\) mm, improving over UD-ToFnet at MAE \(8.88\) mm and RMSE \(11.50\) mm [2511.01704]. On RUD-ToF, it achieved MAE \(16.73\) mm and RMSE \(30.94\) mm, compared with UD-ToFnet at MAE \(17.29\) mm and RMSE \(31.11\) mm [2511.01704]. On FLAT, it achieved MAE \(4.13\) mm and RMSE \(7.35\) mm [2511.01704]. On NYUv2, it improved over DSR-EI across \(\times4\), \(\times8\), and \(\times16\) upsampling factors, reporting \(2.85/0.47\), \(12.8/1.16\), and \(52.3/2.58\) in MAE/RMSE, respectively [2511.01704].

The ablations are especially revealing. Removing fractional calculus degrades performance, as does removing continuous convolution, and a variable-order model outperforms fixed-order alternatives [2511.01704]. Comparisons with GRU- and LSTM-based iterative refinements also show that the fractional dynamics module achieves the best MAE/RMSE trade-off among the tested recurrent refinements while adding only \(0.18\)M parameters over the UD-ToFnet backbone [2511.01704]. These results support the claim that the benefit is not reducible to generic recurrence alone.

A broader implication is that LFRD2 should be viewed less as a ToF-specific architecture than as a learnable PDE layer for inverse problems where one has an initial estimate and benefits from memory-aware iterative refinement. The paper itself makes this point by replacing the initial-state builder according to task and using the same fractional refinement mechanism “for under-display ToF imaging and beyond” [2511.01704]. This suggests that the core transferable ingredient is the learnable fractional reaction-diffusion update, not the upstream sensor-specific encoder.

A common misconception is that LFRD2 is simply a classical reaction-diffusion network with one extra hyperparameter. The surrounding literature indicates otherwise. Fractional derivatives alter the temporal operator from local to nonlocal, fractional spatial operators alter transport from Gaussian to heavy-tailed or skewed propagation, and distributed-order variants introduce multi-scale memory unavailable to single-order models [1001.2293; 1502.06127; 1211.0063; 2503.13096]. LFRD2 adopts only part of this full fractional toolkit, but it does so in a way that makes the memory law itself learnable [2511.01704].

Another misconception is that learnability necessarily sacrifices interpretability. In LFRD2 the neural components parameterize diffusivity, order selection, and convolution kernels, but the update law remains an explicit discretization of a fractional reaction-diffusion PDE [2511.01704]. This does not make the model analytically transparent in the same sense as closed-form Mittag-Leffler solutions, but it preserves a strong operator-level interpretation absent from purely black-box regressors.

The main unresolved issues concern stability, parameter identifiability, and extension beyond the present formulation. The broader literature points to several natural directions: richer distributed-order time operators [1502.06127; 1211.0063], explicit stability regularization based on fractional linearization theory [2507.02094], better bounded-domain operator design informed by well-posedness and invariant sets [2601.21422], and architectures adapted to the limited regularity of fractional PDE solutions on bounded domains [1911.03261]. Within that larger landscape, LFRD2 can be understood as an initial, application-driven instance of a more general program: making fractional reaction-diffusion dynamics trainable while preserving enough PDE structure to remain physically interpretable [2511.01704].

Source: https://www.emergentmind.com/topics/learnable-fractional-reaction-diffusion-dynamics-lfrd2