Abstract:In this paper, we propose and analyze a second-order time-stepping numerical scheme for the inhomogeneous backward fractional Feynman-Kac equation with nonsmooth initial data. The complex parameters and time-space coupled Riemann-Liouville fractional substantial integral and derivative in the equation bring challenges on numerical analysis and computations. The nonlocal operators are approximated by using the weighted and shifted Grünwald difference (WSGD) formula. Then a second-order WSGD scheme is obtained after making some initial corrections. Moreover, the error estimates of the proposed time-stepping scheme are rigorously established without the regularity requirement on the exact solution. Finally, some numerical experiments are performed to validate the efficiency and accuracy of the proposed numerical scheme.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to propose and analyze a new second - order time - weighted shifted Grünwald difference (WSGD) numerical scheme for the inhomogeneous backward fractional Feynman - Kac equation with non - smooth initial data. Specifically, this research focuses on the following aspects: 1. **Numerical scheme design**: Based on the weighted shifted Grünwald formula, a modified time - discretization method is designed for the inhomogeneous term and non - smooth initial data. 2. **Error estimation**: Strict error estimates are established for the proposed discrete scheme, and these estimates do not require the regularity of the exact solution. ### Background and problem description The backward fractional Feynman - Kac equation involved in this research can be expressed as: \[ \begin{cases} \frac{\partial G(x,\rho,t)}{\partial t}={}_0D_{t,x}^{1 - \alpha}\Delta G(x,\rho,t)-\rho U(x)G(x,\rho,t)+f(x,\rho,t),&(x,t)\in\Omega\times(0,T],\\ G(x,\rho,0) = G_0(x),&x\in\Omega,\\ G(x,\rho,t)=0,&(x,t)\in\partial\Omega\times(0,T], \end{cases} \] where ${}_0D_{t,x}^{1 - \alpha}$ represents the Riemann - Liouville fractional substantial derivative, $\alpha\in(0,1)$, $\Delta$ is the Laplace operator, $\Omega$ is a bounded convex polygonal region, and $U(x)$ and $f(x,\rho,t)$ are functions and source terms defined on $\Omega$, respectively. ### Main contributions 1. **New time - discretization method**: Based on the weighted shifted Grünwald formula, a new modified second - order time - discretization method is proposed, which is suitable for cases with inhomogeneous terms and non - smooth initial data. 2. **Strict error estimates**: A strict error estimate analysis is carried out for the proposed discrete scheme, ensuring the effectiveness and accuracy of the numerical method. ### Numerical experiments To verify the effectiveness and convergence rate of the proposed method, the author carried out several numerical experiments. For example, in the one - dimensional case, for the homogeneous problem with non - smooth initial data, the theoretical analysis results are verified by calculating the errors and convergence rates at different time steps. In summary, this paper aims to provide an effective numerical method for solving the challenging backward fractional Feynman - Kac equation, and the accuracy and reliability of the method are verified through strict mathematical analysis and numerical experiments.

Analysis of a WSGD scheme for backward fractional Feynman-Kac equation with nonsmooth data

Generalized Finite Difference Method for a Class of Multidimensional Space-Fractional Diffusion Equations

Second-Order Time Stepping Scheme Combined with a Mixed Element Method for a 2D Nonlinear Fourth-Order Fractional Integro-Differential Equations

A fast implicit difference scheme for solving the generalized time-space fractional diffusion equations with variable coefficients

Analysis of the L1 scheme for fractional wave equations with nonsmooth data

Fast second-order implicit difference schemes for time distributed-order and Riesz space fractional diffusion-wave equations

Stability analysis of a second-order difference scheme for the time-fractional mixed sub-diffusion and diffusion-wave equation

High‐order finite difference/spectral‐Galerkin approximations for the nonlinear time–space fractional Ginzburg–Landau equation

A second-order L 2 - 1 σ difference scheme for the nonlinear time–space fractional Schrödinger equation

Second‐order difference scheme for the time fractional mixed diffusion‐wave equation with initial weak regularity

Numerical methods for forward fractional Feynman–Kac equation

High-Order Dissipation-Preserving Methods for Nonlinear Fractional Generalized Wave Equations

An efficient second order stabilized scheme for the two dimensional time fractional Allen-Cahn equation

Analysis of a local discontinuous Galerkin scheme for fractional Korteweg-de Vries equation

Galerkin finite element approximation for semilinear stochastic time-tempered fractional wave equations with multiplicative white noise and fractional Gaussian noise

Dissipative WENO stabilization of high-order discontinuous Galerkin methods for hyperbolic problems

Mathematical Analysis and a Second-Order Compact Scheme for Nonlinear Caputo–Hadamard Fractional Sub-diffusion Equations

A Discontinuous Galerkin Method and Its Error Estimate for Nonlinear Fourth-Order Wave Equations.

Semi-implicit Galerkin–Legendre Spectral Schemes for Nonlinear Time-Space Fractional Diffusion–Reaction Equations with Smooth and Nonsmooth Solutions

Improved Uniform Error Bounds of Exponential Wave Integrator Method for Long-Time Dynamics of the Space Fractional Klein-Gordon Equation with Weak Nonlinearity

Convergence analysis of a fast second‐order time‐stepping numerical method for two‐dimensional nonlinear time–space fractional Schrödinger equation