New type of the unique continuation property for a fractional diffusion equation and an inverse source problem

In this work, a new type of the unique continuation property for time-fractional diﬀusion equations is studied. The proof is mainly based on the Laplace transform and the properties of Bessel functions. As an application, the uniqueness of the inverse problem in the simultaneous determination of spatially dependent source terms and fractional order from sparse boundary observation data is established


Introduction
Classical unique continuation, which shows that the local information of the solution can uniquely determine the global information of the solution, is an important property of elliptic and parabolic equations.More specifically, the vanishment of the solution of the homogeneous problem in the open subset results in its vanishment in the entire domain, see Saut and Scheurer [1].There are many important applications of the unique continuation properties.For the inverse source problem, we refer the readers to El Badia and Ha-Duong [2], Hu, Kian, and Zhao [3].For approximate controllability, see Cannarsa, Tort, and Yamamoto [4], Dou and Lu [5].
Major advances have been made in fractional calculus in the last few decades.As a generalization of the classical diffusion equation, the fractional diffusion equation has become a research hotspot in mathematics.The extension of the unique continuation properties to fractional diffusion equations has attracted the attention of many researchers.Li and Yamamoto [6] investigated the lateral Cauchy problem for the one-dimensional time-fractional diffusion equation.As a direct conclusion of the uniqueness of the Cauchy problem, they proved that the classical unique continuation property is valid.Other researchers in recent years have investigated the unique continuation property for fractional differential equations via Carleman estimates, see Xu, Cheng, and Yamamoto [7], Cheng, Lin, and Nakamura [8].Sakamoto and Yamamoto [9] indicated the weak unique continuation property of the time-fractional diffusion equations with the homogeneous Dirichlet boundary condition on the whole boundary.Jiang et al. [10] generalized the result in [9] to the multiterm case.The literature mentioned above focuses on the unique continuation properties in the inner open subset.However, as far as we know, there are few published results concerning the properties from sparse data.We wonder whether the global information of the solution can be determined from a small finite number N of measurement points.The primary goal of this paper is to investigate the unique continuation property of the following problem.
Let ⊂ R 2 be a unit disc and T > 0. Consider the following initial-boundary value problem of the time-fractional diffusion equation: ( c ∂ α t u denotes the Caputo derivative of u at time t > 0, which is defined by Here 0 < α < 1. (•) denotes the gamma function.For the definition and properties of fractional derivatives, see Podlubny [11].
Problem 1.1 Let u be the solution of (1) and Here, the boundary flux data ∂u ∂ν are used, and ν is the unit outward normal vector of ∂ .
In practice, environmental authorities often need to determine the intensity and location of pollution sources based on monitoring data.The study of the inverse source problem has become popular due to the aforementioned issues.Anomalous transport poses significant challenges for accurate prediction and remediation of groundwater contamination.Fractional calculus has attracted more and more attention in anomalous diffusion due to its heritability and memorability.We refer the readers to Zhang, Meerschaert, and Baeumer [12], Sun et al. [13], and Yin et al. [14].Furthermore, it is known that the fractional order is related to the inhomogeneity of the media, but it is not clear which physical law can relate the inhomogeneity to the fractional order.So we are also required to consider the inverse fractional order problem.There exists a large and rapidly growing number of publications related to the inverse problems in determining sources, fractional orders, and other unknown coefficients.For some early work on the determination of the source terms, we refer the readers to Zhang and Xu [15], Kirane and Malik [16], Chi, Li, and Jia [17], Liu, Rundell, and Yamamoto [18]; and for more recent works, we refer to Liu and Zhang [19], Rundell and Zhang [20], Li and Zhang [21], Phuong, Kumar, and Binh [22], Binh and Long [23], Phuong, Thi, and Luc [24] and the references therein.In particular, we mention the reference [10] by Jiang et al., where the uniqueness of an inverse problem in determining the spatial component in the source term by interior measurements utilizing the weak unique continuation property was proven.For the determination of the fractional orders and other unknown coefficients, refer to Cheng et al. [25], Li et al. [26], Li et al. [27], Kian et al. [28], Ozbilge and Demir [29], Jday and Mdimagh [30], Phuong et al. [31], Phuong et al. [32], Long and Saadati [33].It is worth noting that some of the above inverse problems are studied in disc or rectangular domains mathematically.Inversion studies in these special domains have also received widespread attention in engineering defect identification (see [34][35][36] for instance).The study of inverse problems in these special regions can provide theoretical support for engineering numerical simulations.
This paper considers the inverse problem of simultaneously determining spatially dependent source terms and fractional orders from two point Neumann boundary data.More precisely, we consider an initial-boundary value problem: ( The above model explains the evolution of density u(x, t) at location x and time t of some substances, such as contaminants, where g(t)f (x) is a source term in .The source term can often be assumed to be modeled in the form of the separation of variables.g(t) and f (x) describe the spatial distribution of the source and the time evolution pattern, respectively.We are dedicated to using sparse boundary data to simultaneously determine the spatially dependent source term f (x) and fractional order α.
A precise mathematical statement of this inverse problem is provided below.
Problem 1.2 Let be the unit disc in R 2 and the fractional derivative be defined as above.
Assume that the temporal component g(t) is known in equation ( 3).The boundary flux data are given at two points: Can we uniquely determine f (x) and α simultaneously?Here, Z ob is defined as in Problem 1.1.
It is worth noting that the observations required in this work are only for a limited period (i.e., t ∈ (0, T)), which is in contrast to Li and Zhang [21,Theorem 1.1].Our aim is to find the connection between the above inverse problem and Problem 1.1, and to apply the conclusion of Problem 1.1 to prove the uniqueness of this inverse problem.
We give the following uniqueness theorem to answer Problem 1.2.
Denote the two sets of unknown solutions to equation (3) as u 1 and u 2 .If The rest of this paper is organized as follows.In Sect.2, some prior knowledge is listed, such as the eigensystem of the Laplacian operator, the properties of Bessel functions, and the Mittag-Leffler function.In Sect.3, we give the proof of Theorem 1.1.In Sect.4, we present the proof of Theorem 1.2 as an application of Theorem 1.1.Finally, concluding remarks are provided in Sect. 5.

Preliminary information
In this section, we first set up notations and introduce the Dirichlet eigensystem of the Laplacian operator.Let L 2 ( ) be a usual L 2 -space with the inner product •, • , and let H 1 0 ( ) denote the usual Sobolev spaces.We introduce the eigensystem {(λ n , ϕ n )} ∞ n=-∞ of the Laplacian operator -on with the Dirichlet boundary condition, i.e., Since is the unit disc in R 2 as mentioned above, we consider the eigensystem in polar coordinates for convenience.According to the Bessel function and its related properties, { λ n , ϕ n } ∞ n=0 is given as follows: and ϕ n denotes the corresponding eigenfunction which forms a complete orthonormal basis of L 2 ( ).Here, (r, θ ) are the polar coordinates on , demonstrates the dependence of m on n such that m(n) = -m(-n), and ω n is the normalized coefficient and allows the form See [21, Sect.2.2] for further details.Sometimes, we also write λ n = λ |n| for simplicity.Therefore, for a given eigenvalue λ n 0 , in the case of m(n 0 ) = 0, the corresponding eigenpairs are given as For the purpose of the latter proof, we list some properties of Bessel functions.
Lemma 2.1 [37] Bessel functions have the following differential relations and recurrence relations: Next, we introduce the Mittag-Leffler function and some of its properties.The generalized Mittag-Leffler function plays an important role in fractional calculus [38].The function is defined as follows: We provide the asymptotic property of the Mittag-Leffler function below.
Then, we introduce the lemma for the Laplace transform of the Mittag-Leffler function.First, we define s ∈ C + := {s ∈ C : Re(s) > 0} to ensure analyticity.Here, s is the Laplace transform parameter.

A new type of unique continuation
In this section, we present the proof of Theorem 1.1.From the above lemma, we are ready to give the proof of the main result.
Proof of Theorem 1.1 We first obtain the eigenfunction expansions of the solution to (1) by the Fourier method.We multiply both sides of (1) by ϕ n (x) and integrate the equation with respect to x.Using integration by parts for the second term and (11) where u n (t) = u(•, t), ϕ n .By using the Mittag-Leffler function, we can formally obtain the expansion Together with (6), the following holds: Next, we assert that u(r, θ , t) is t-analytic.To validate this claim, since E α,1 (-λ n t α ) is analytic in C, we denote Thus, for any K ⊂⊂ {z ∈ C; Re z > 0}, is analytic in K .According to Lemma 2.2 and Lemma 2.3, the following holds: Thus, lim n→∞ u +u + N L ∞ (K;L 2 ( )) = 0 so that u + (r, θ , t) is analytic in K .Analogously, we can obtain that u -is also analytic in K .As a result of the arbitrariness of K , the above assertion is valid.A similar argument can be found in Sakamoto and Yamamoto [9].
From the analyticity of u(r, θ , t), the solution u to the initial-boundary value problem (1) can be analytically extended from (0, T) to (0, ∞).For simplicity, we still denote the extension by u.Therefore, we arrive at the following initial-boundary value problem: and the condition ∂u ∂ν (z 1 , t) Next, we perform the Laplace transform in ( 13) and use the formula to derive the transformed equation From the eigensystem { λ n , ϕ n } ∞ n=-∞ of -on with the Dirichlet boundary condition, the following holds: By using the formula •; s) and Lemma 2.1, the following holds: and together with ( 14), we have Note that is a unit disk, as mentioned in Sect. 1. Thus z 1 and z 2 are boundary observation points implying r = 1, there holds and We set η = s α , which yields It is readily seen that ( 16) holds for η ∈ C\{-λ n } ∞ n=1 .Then, for any n = 1, 2, . . ., we can take a sufficiently small circle centered at -λ n that does not include distinct eigenvalues.Integrating (16) on this circle yields Moreover, from the properties of Bessel functions in Lemma 2.1, it follows that is the Bessel function of order |m(n)| with √ λ n as its zero point, the above formula can be reduced to By noting the definition of ω n , we further see that By combining the above formulas, we finally get We divide the next proofs into the following two cases.Case 1. Provided that m(n) = 0, that is, n = 0, u 0 , ϕ 0 = 0 holds.Case 2. Provided that m(n) = 0, the following hold: and We conclude that From the matrix determinant and Euler's formula, we conclude Recalling Theorem 1.1, θ 1 and θ 2 are fulfilled ∈ πQ, Q is the set of rational numbers.

Fractional order and source term identification
In this section, our goal is to give the proof for the uniqueness of the inverse problem for the determination of the fractional order and source term.
Proof of Theorem 1.2 To prove Theorem 1.2, we first prove that the observation To achieve this, we extend the function g from (0, T) to (0, ∞) by letting g = 0 outside of (0, T).For simplicity, we still denote the extension by g and obtain the following initial-boundary value problem: The solution u j to the initial-boundary value problem ( 19) is an extension of the solution u to problem (1.2).For simplicity, we still denote the extension as u.By combining equation (18) and the series representation of the solution u (e.g., see Sakamoto and Yamamoto [9]), we have We perform Laplace transform on both sides of the above formula to obtain which implies Moreover, from the observation By choosing ε ∈ (0, T) and multiplying se sε on both sides of the above formula, we obtain We note from the initial value theorem of the Laplace transforms that lim Re s→∞ sĝ(s) = g(0) > 0, and then If α 1 = α 2 , then from Lemma 3.7 in Li and Zhang [21], it follows that f 1 , ϕ n = f 2 , ϕ n = 0 for all n ∈ Z, which contradicts f = 0. Therefore, we must have α 1 = α 2 .Now, it remains to prove that if ∂u ∂ν (z 1 , t) = ∂u ∂ν (z 2 , t) = 0 holds, then f (x) = 0 can be derived.Let u satisfy the initial-boundary value problem (3), where g ∈ C 1 [0, T] and f ∈ L 2 ( ).According to Lemma 2.4, u takes the form of (9).Therefore, the following holds: In particular, then the following holds: Applying the operator J 1-α to both sides of the above formula yields According to Gronwall's inequality, we obtain ∂v ∂ν (z , t) = 0, z ∈ ∂ , and = 1, 2. Lastly, we apply Theorem 1.1 to the homogeneous problem (10) to derive v = 0 in × (0, T), in which we can conclude f = 0 in .We therefore finish the proof of the Theorem 1.2.

Concluding remarks and future work
In this paper, a new type of the unique continuation property for time-fractional diffusion equations is studied.The proof of this unique continuation principle was concerned, which leads to uniqueness of the corresponding inverse problems.We gave a representation formula for the solution in polar coordinates based on eigenfunction expansion and the Mittag-Leffler function.With the help of the Laplace transformation and the properties of the Bessel function, we proved the unique continuation property from two point Neumann boundary data.As an application, we considered the uniqueness of an inverse problem for simultaneously determining the spatial component in the source term and the fractional order from two point finite time observation data.The uniqueness result is slightly sensitive to the geometry of the domain.We will consider the case in other domains such as rectangular and elliptic in the future.Otherwise, it would be interesting to investigate the unique continuation properties of the solution in a high dimensional case.