On the Nonlocal Problems in Time for Time-fractional Subdiffusion Equations

The nonlocal boundary value problem, dtρu(t)+Au(t)=f(t) (0&lt;ρ&lt;1, 0&lt;t≤T), u(ξ)=αu(0)+φ (α is a constant and 0&lt;ξ≤T), in an arbitrary separable Hilbert space H with the strongly positive selfadjoint operator A, is considered. The operator dt on the left hand side of the equation expresses either the Caputo derivative or the Riemann-Liouville derivative; naturally, in the case of the Riemann - Liouville derivatives, the nonlocal boundary condition should be slightly changed. Existence and uniqueness theorems for solutions of the problems under consideration are proved. The influence of the constant α on the existence of a solution to problems is investigated. Inequalities of coercivity type are obtained and it is shown that these inequalities differ depending on the considered type of fractional derivatives. The inverse problems of determining the right-hand side of the equation and the function φ in the boundary conditions are investigated.


Introduction
Let H be a separable Hilbert space with the scalar product (·, ·) and the norm || · || and A : H → H be an arbitrary unbounded positive selfadjoint operator in H. Suppose that A has a complete in H system of orthonormal eigenfunctions {v k } and a countable set of positive eigenvalues λ k . It is convenient to assume that the eigenvalues do not decrease as their number increases, i.e., 0 < λ 1 ≤ λ 2 · ·· → +∞.
Using the definitions of a strong integral and a strong derivative, fractional analogues of integrals and derivatives can be determined for vector-valued functions (or simply functions) h : R + → H, while the well-known formulae and properties are preserved (see, e.g., [1]). Recall that the fractional integration of order σ < 0 of the function h(t) defined on [0, ∞) has the form If in this definition we interchange differentiation and fractional integration, then we obtain the definition of the regularized derivative, that is, the definition of the fractional derivative in the sense of Caputo: Note that if ρ = 1, then fractional derivatives coincides with the ordinary classical derivative of the first order: ∂ t h(t) = D t h(t) = d dt h(t). Let ρ ∈ (0, 1) be a fixed number and let C((a, b); H) stand for a set of continuous functions u(t) of t ∈ (a, b) with values in H.
The subject of this work is the following two nonlocal boundary value problems: and where f (t), g(t) ∈ C((0, T]; H), ϕ, φ ∈ H and α is a constant, ξ-fixed point. These problems are also called the forward problems.
The definition of the solution to the nonlocal problem (3) is introduced in a similar way.
If α = 0 (and ξ = T), then these problems are called the backward problems. The backward problems in case (2) were studied in detail, for example, in [2][3][4]. The work [5] is devoted to the study of the backward problem in case (3). Therefore, in what follows we only consider the case α = 0.
The backward problems for the diffusion process are of great importance in engineering fields and are aimed at determining the previous state of a physical field (for example, at t = 0) based on its current information (see, e.g., [3] and for the classical head equation see [6]). However, regardless of the fact that the Riemann-Liouville or the Caputo derivative is taken into the equation, this problem is ill-possed in the sense of Hadamard. In other words, a small change of u(T) in the norm of space H leads to large changes in the initial data. As can be seen from the main results of papers [2][3][4][5] (note, in these works 0 < ρ < 1), the situation changes if we take the sufficiently smooth function u(T). Since the problem is ill-posed, many authors have considered various regularization options for finding the initial condition (see, for one-dimensional elliptical part, Liu and Yamamoto [2], for the nonlinear case, Tuan, Huynh, Ngoc, and Zhou [7]). In particular, as for numerical approaches, see Tuan, Long and Tatar [8], Wang and Liu [9] and the references therein.
In the case ρ = 1 these problems are also called (see, e.g., [6], p. 214) the inverse heat conduction problem with inverse time (retrospective inverse problem). It should also be noted that, in this case, even the smoothness of the function u(T) does not guarantee the stability of the solution (see, e.g., Chapter 8.2 of [6]).
As we know, in most models described by differential (and pseudodifferential, see e.g., [10]) equations the initial condition is used. However, in practice, some other models have to use nonlocal conditions, for example, including integrals over time intervals (see, e.g., [11] for reaction-diffusion equations or [12] for fractional equations), or connecting the solution at different times, for instance at the initial time and at the terminal time (see, e.g., [13,14]).
Note, nonlocal conditions express and explain some full details about natural events because they consider additional information in the initial conditions.
The following nonlocal boundary value problem for the classical diffusion equation in an arbitrary Banach space E with the strongly positive operator A, has been extensively studied by numerous researchers (see, e.g., A. O. Ashyralyev et al. [13,14]). As shown in these papers, in contrast to the retrospective inverse problem, the problem (5) is coercively solvable in some spaces of differentiable functions. It should also be noted that various nonlocal boundary value problems for parabolic equations reduce to the boundary value problem (5) (see, e.g., [15], Chapter 1).
In the present paper we prove the existence and uniqueness theorems for solutions of problems (2) and (3). Next, we will study the dependence of the existence of a solution on the value of the parameter α. We will also prove, in contrast to the backward problems, that the solutions of problems (2) and (3) continuously depend on the right-hand side of the equation and on the function ϕ. Inequalities of coercivity type are obtained and it is shown that these inequalities differ depending on the considered type of fractional derivatives. The inverse problems of determining the right-hand side of the equation and function ϕ in the boundary conditions are investigated.
The inverse problem of determining the source function f with the final time observation has been well studied and much theoretical research has been published for classical partial differential equations. As a monograph, we should refer to Kabanikhin [6] and Prilepko, Orlovsky, and Vasin [16]. As for the fractional differential equations, one can construct theories parallel to [6,16], and the works are now made continuously. We give a brief overview of work on this inverse problem at the beginning of Section 4.
The inverse problem of defining the function ϕ arises from real-life processes. For example, when the initial temperature and final temperature for the heat equation are not indicated immediately, and it is not required to find, but information about the difference between the initial and final temperatures is sought. To the best of our knowledge, such an inverse problem was discussed only in the paper [17]. The authors considered this problem for the subdiffusion equation including the Caputo fractional derivative, the elliptical part of which is a two-variable differential expression with constant coefficients.
The remainder of this paper is composed of four sections and the Conclusion. In the next section, we introduce the Hilbert space associated with the degree of operator A and recall some properties of the Mittag-Leffler functions. Section 3 is devoted to the study of the nonlocal problem (2). Here, we first investigate problem (2) for the homogeneous equation, and then move on to the main problem. In Section 4, we study the inverse problem of determining the right-hand side of Equation (2). In this case, we assume that the unknown function f does not depend on t. The next section is devoted to the study of the inverse problem for the determination of the boundary function ϕ. Since problems (2) and (3) are studied in a similar way, in Section 6 we present only the main points of the proof of the theorem on the existence and uniqueness of the solution to problem (3). Inverse problems for Equation (3) are considered in the same way as inverse problems for Equation (2). Therefore, we omit these details.

Preliminaries
In this section, we introduce the Hilbert space of "smooth" functions related to the degree of operator A and recall some properties of the Mittag-Leffler functions, which we will use in what follows.
Let τ be an arbitrary real number. We introduce the power of operator A, acting in H according to the rule where h k are the Fourier coefficients of a function h ∈ H: h k = (h, v k ). Obviously, the domain of this operator has the form For elements of D(A τ ) we introduce the norm and together with this norm D(A τ ) turns into a Hilbert space. For 0 < ρ < 1 and an arbitrary complex number µ, by E ρ,µ (z) we denote the Mittag-Leffler function with two parameters (see, e.g., [18], p. 12): If the parameter µ = 1, then we have the classical Mittag-Leffler function: In what follows we need the asymptotic estimate of the Mittag-Leffler function with a sufficiently large negative argument. The well known estimate has the form (see, e.g., [19], p. 136) where µ is an arbitrary complex number. This estimate essentially follows from the following asymptotic estimate (see, e.g., [19], p. 134): For the Mittag-Leffler function with two parameters E ρ,ρ (−t) one can obtain a better estimate than (7). Indeed, using the asymptotic estimate (see, e.g., [19], p. 134) and the fact that E ρ,ρ (t) is real analytic, we can obtain the following inequality [5] We will also use a coarser estimate with positive number λ and 0 < ε < 1: which is easy to verify. Indeed, let t ρ λ < 1, then t < λ −1/ρ and If t ρ λ ≥ 1, then λ −1 ≤ t ρ and Proof. For x ≥ 0 this is obvious; estimates (12) follow from definition (6). For x < 0 we use the integral representation (see, e.g., [20], p. 54) Using Proposition 1, by virtue of estimates (12) and equality E ρ (0) = 1, we arrive at (see [20], p. 47).

Proposition 2.
The Mittag-Leffler function of the negative argument E ρ (−x) is monotonically decreasing function for all 0 < ρ < 1 and Proposition 3. Let ρ > 0 and λ > 0. Then for all positive t > 0 one has [3] t 0 Proof. First, we calculate the derivative of the Mittag-Leffler function Note that here the series is termwise differentiable in R. Now, by virtue of the equality we obtain the required result.

Proposition 5.
Let 0 < ρ < 1 and λ > 0. Then for all positive t one has On the other hand, using the properties of Euler's beta function B(a, b), we obtain By virtue of the definition of the Mittag-Leffler function E ρ (z) this implies the statement of the proposition.

Well-Posedness of Problem (2)
To solve problem (2), we divide it into two auxiliary problems: where ψ ∈ H is a given function.
Problem (16) is a special case of problem (2), and the solution to problem (15) is defined similarly to Definition 1.
If ψ = ϕ − ω(ξ) and ω(t) and w(t) are the corresponding solutions, then it is easy to verify that function u(t) = ω(t) + w(t) is a solution to problem (2). Therefore, it is sufficient to solve the auxiliary problems.
For problem (15) we have the following statement.
. Then problem (15) has a unique solution and this solution has the representation Moreover, there is a constant C ε > 0 such that the following coercive type inequality holds: Proof of Theorem 1. It is not hard to verify that the series (17) is a formal solution to problem (15) (see, e.g., [20], p. 173, [21]). In order to prove that function (17) is actually a solution to the problem, it remains to substantiate this formal statement, i.e., to show that the operators A and D ρ t can be applied term-by-term to series (17). Let S j (t) be the partial sum of series (17). Then Due to the Parseval equality we may write Then, by inequality (11) for 0 < ε < 1 one has or, by virtue of the generalized Minkowski inequality, Hence, we obtain Aω(t) ∈ C([0, T]; H) and in particular ω(t) ∈ C([0, T]; H).
Therefore, from the above reasoning, we have D ρ t ω(t) ∈ C((0, T]; H) and Thus, we have completed the rationale that (17) is a solution to problem (15). The last two inequalities imply estimate (18).
The uniqueness of the solution can be proved by the standard technique based on completeness of the set of eigenfunctions {v k } in H (see, e.g., [5]).
Theorem 1 is completely proved.
If f does not depend on t, then the statement of Theorem 1 is true for all f ∈ H.
Corollary 1. Let f ∈ H. Then problem (15) has a unique solution and this solution has the representation Moreover, there is a positive constant C such that the following coercive type inequality holds: Proof. Since f does not depend on t, then we have the following form for the Fourier coefficients of ω (see (17)) Application of Formula (14) to the integral shows that the formal solution to problem (15) has the form (20).
Let S j (t) be the partial sum of series (20). Then by virtue of estimate (7), we obtain Now, using this estimate and repeating the arguments similar to the proof of Theorem 1, it is easy to check that (20) is indeed a solution to problem (15) and estimate (21) holds true.
We now turn to problem (16). In accordance with the Fourier method, we will look for a solution to problem (16) in the form of a series: where T k (t), k ≥ 1, are solutions of the nonlocal problems: where ψ k are the Fourier coefficients of function ψ ∈ H. Let us denote T k (0) = b k . Then the unique solution to the differential Equation (23) with this initial condition has the form T k (t) = b k E ρ (−λ k t ρ ) (see, e.g., [20], p. 174). From the nonlocal conditions of (23) we obtain the following equation to find the unknown By virtue of property (13) of the Mittag-Leffler function, E ρ (−λ k ξ ρ ) = α for all α ≥ 1 and α < 0 (note, ξ > 0 and λ k > 0). Therefore, from (24) we have here and below, by C β we will denote a constant depending on β, not necessarily the same one. If α = 0, then E ρ (−λ k ξ ρ ) = 0, but the Mittag-Leffler function can asymptotically tend towards zero (see (8)). Therefore, in this case one has: This case, as noted above (see (4)), has been studied in detail in [2][3][4]. Let 0 < α < 1. Then according to Proposition 2, there is a unique λ 0 > 0 such that Thus, if α / ∈ (0, 1) or α ∈ (0, 1), but λ k = λ 0 for all k ≥ 1, then the formal solution of problem (16) has the form Finally, let 0 < α < 1 and λ k = λ 0 for k = k 0 , k 0 + 1, . . . , k 0 + p 0 − 1, where p 0 is the multiplicity of the eigenvalue λ k 0 . Then the nonlocal problem (23) has a solution if the boundary function ψ(x) satisfies the following orthogonality conditions and for these k ∈ K 0 arbitrary numbers b k are solutions of Equation (24). For all other k Thus, the formal solution of problem (16) in this case has the form Throughout what follows we will assume that whenever 0 < α < 1 and λ k = λ 0 , then orthogonality condition (27) is satisfied.
Let us show that the operators A and D ρ t can be applied term-by-term to series (26); for series (29) this question is considered in a completely similar way.
Let S j (t) be the partial sum of series (26). Then Due to the Parseval equality, we may write Using estimates (7), (25) and (28) we obtain Therefore if ψ ∈ H, then Au(t) ∈ C((0, T]; H). From Equation (16) one has D ρ t u(t) = −Au(t), t > 0, and the above estimates imply which means D ρ t w(t) ∈ C((0, T]; H). For S j (t), taking into account estimate (7), we obtain Hence w(t) ∈ C([0, T]; H), which was required by the definition of the solution to problem (16).
Let us investigate the uniqueness of the solution to problem (16). Suppose we have two solutions: w 1 (t), w 2 (t) and set w(t) = w 1 (t) − w 2 (t). Then, we have Let w k (t) = (w(t), v k ). Since the operator A is self-adjoint, one has and the nonlocal condition implies Let us denote w k (0) = b k . Then the unique solution to the differential Equation (35) with this initial condition has the form w k (t) = b k E ρ (−λ k t ρ ) (see, e.g., [20], p. 174). From the nonlocal conditions of (36) we obtain the following equation to find the unknown Let first α / ∈ (0, 1) or α ∈ (0, 1), but λ k = λ 0 for all k ≥ 1. Then E ρ (−λ k ξ ρ ) = α for all k. Consequently, in this case all b k are equal to zero (therefore w k (t) = 0), and by virtue of completeness of the set of eigenfunctions {v k }, we conclude that w(t) ≡ 0. Thus, problem (16) in this case has a unique solution. Now suppose that α ∈ (0, 1) and λ k = λ 0 , k ∈ K 0 . Then E ρ (−λ k ξ ρ ) = α, k ∈ K 0 and therefore Equation (37) has the following solution: b k = 0 if k / ∈ K 0 and b k is an arbitrary number for k ∈ K 0 . Thus, in this case, there is no uniqueness of the solution to problem (16).
If α ∈ (0, 1) and λ k = λ 0 , k ∈ K 0 , then we assume that the orthogonality conditions (42) are satisfied. The solution of problem (16) has the form (40) with arbitrary coefficients b k , k ∈ K 0 . Moreover, there are constants C α > 0 and C ε > 0 such that the following coercive type inequality holds: The results of Theorems 2 and 3 are based on the assumption of orthogonality conditions (27) and (42) correspondingly. The question naturally arise, to what extent do these assumtions limit? In order to answer this question, consider the following example.
Let Ω ⊂ R N be a bounded domain with sufficiently smooth boundary ∂Ω and denote by A 0 the operator in is the Laplace operator. Then (see, e.g., [22]) A 0 has a complete in L 2 (Ω) system of orthonormal eigenfunctions {v k (x)} and a countable set of nonnegative eigenvalues λ k (→ +∞), and λ 1 = λ 1 (Ω) > 0.
Let A stand for the operator, acting as Then it is not hard to verify, that A is a positive self-adjoint extension in L 2 (Ω) of operator A 0 . Therefore, one can apply Theorems 2 and 3 to operator A and consequently to the problem: Suppose α ∈ (0, 1) and ξ ∈ (0, T] are such that Then problem (44) has the unique solution for any function ψ ∈ L 2 (Ω), which satisfies the orthogonality condition or in other words, for any function

Inverse Problem of Determining the Heat Source Density
The inverse problems of determining the right-hand side (the heat source density) of various subdiffusion equations have been considered by a number of authors (see, e.g., survey papers [23] and the bibliography therein). However, there is no general closed theory for the abstract case of the source function F(x, t). Known results deal with the separated source term F(x, t) = h(t) f (x). The appropriate choice of the overdetermination depends on the choice whether the unknown is h(t) or f (x).
The most difficult to study is the case when the function h(t) is unknown (see the survey work [3,23] for the case of subdiffusion equations, and, for example, monographs [6,16] and [24][25][26] for the classical heat equation). In inverse problems of this type, the condition u(x 0 , t) = u 0 (t) is taken as an additional one, and the operator A does not depend on t. In all the above-mentioned works on the subdiffusion equations, only the uniqueness of the solution to the inverse problem was proved.
Uniqueness questions in the inverse problem of finding the function has been studied in, e.g., [27,28]. Paper [27] shows the uniqueness result for f (x) if h(t) > 0 in a case of a self-adjoint operator A with time-independent coefficients. Article [28], Example 3.1, shows the non-uniqueness result if h(t) changes its sign. The time-dependent operator A(t) has been considered in [28], Theorem 3.1, where uniqueness was established under some monotonic in time assumptions on h(t) and A(t). In paper [29], the author studies a general multi-dimensional case for a linear time-fractional partial differential equations with time-dependent coeffitients. The author uses separation of variable technique combined with the maximum principle to derive the uniqueness result assuming h ∈ C[0, T] and 0 ≡ h ≥ 0. In this result, we do not need any monotonic premises on h(t) and/or the coefficients of A(t), which is the new aspect (and highlight) in this area of inverse problems. Many authors have considered an Equation (45) in which h(t) ≡ 1 and f (x) is unknown (see, e.g., [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48]). Let us mention only some of these works. The case of subdif-fusion equations, the elliptic part of A of which is an ordinary differential expression, is considered in [30][31][32][33][34][35][36]. The authors of the articles [37][38][39][40][41] studied subdiffusion equations in which the elliptic part of A is either a Laplace operator or a second-order operator. The paper [42] studied the inverse problem for the subdiffusion Equation (2) with the Cauchy condition. In this article [42] and most other articles, including [37][38][39][40], the Caputo derivative is used as a fractional derivative. The recent articles [43,44] are devoted to the inverse problem for the subdiffusion equation with Riemann-Liouville derivatives. In [31,41], the fractional derivative in the subdiffusion equation is a two-parameter generalized Hilfer fractional derivative; this type of fractional derivative contains a parameter belonging to the interval [0, 1], and its extreme values correspond to the Caputo and Riemann-Liouville derivatives. Various models of applied problems leading to Hilfer fractional derivatives are investigated in [49]. Note also that the papers [31,37,40] contain a survey of papers dealing with inverse problems of determining the right-hand side of the subdiffusion equation.
In [41,46,47], non-self-adjoint differential operators (with nonlocal boundary conditions) were taken as A, and solutions to the inverse problem were found in the form of biortagonal series.
In their previous work [48], the authors of this article considered an inverse problem for simultaneously determining the order of the Riemann-Liouville fractional derivative and the source function in the subdiffusion equations. Using the classical Fourier method, the authors proved the uniqueness and existence of a solution to this inverse problem.
It should be noted that in all of the listed works, the Cauchy conditions in time are considered (an exception is work [45], where the integral condition is set with respect to the variable t). In the present paper, for the best of our knowledge, an inverse problem for subdiffusion equation with a nonlocal condition in time (see (46)) is considered for the first time.
The papers [50,51] deal with the inverse problem of determining the order of the fractional derivative in the subdiffusion equation and in the wave equation, respectively.
Let us consider the inverse problem with the additional condition in which the unknown element f ∈ H, characterizing the action of heat sources, does not depend on t and Ψ, ϕ ∈ H are given elements, α is an arbitrary given constant. Note that if τ = ξ, then the nonlocal condition in (46) coincides with the Cauchy condition u(0) = ϕ 1 (see (4)). In this case, this inverse problem was studied in [42]. In what follows we shall deal only with the case α ≥ 1, since in this case the uniqueness of the solution is relatively easy to prove.
Proof of Theorem 4. Existence. If f is known, then the unique solution of problem (46) has the form (39), and since f does not depend on t, then, thanks to formulas (22) and (14), it is easy to verify that the formal solution of problem (46) has the form (49). By virtue of additional condition (47) and completeness of the system {v k } we obtain: After simple calculations, we obtain With these Fourier coefficients we have the above formal series (48) for the unknown Let us show the convergence of series (48). If F j are the partial sums of series (48), then by virtue of the Parseval equality we may write Since ξ > 0, then E ρ (−λ k τ ρ )ξ ρ E ρ,ρ+1 (−λ k ξ ρ ) > 0. Therefore, Using the asymptotic estimate (see (8)) we obtain Since τ > 0 and α ≥ 1, then By virtue of (52), Thus, if ϕ, Ψ ∈ D(A), then from estimates of I i,j and (51) we obtain f ∈ H.
After finding the unknown function f ∈ H, the fulfilment of the conditions of Definition 2 for function u(t), defined by series (49) is proved in exactly the same way as with Corollary 1 and Theorem 2.
Uniqueness. Suppose we have two solutions: {u 1 (t), f 1 } and {u 2 (t), f 2 }. It is required to prove u(t) ≡ u 1 (t) − u 2 (t) ≡ 0 and f ≡ f 1 − f 2 = 0. Since the problem is linear, to determine u(t) and f we have the problem: where ξ-fixed point. Let u(t) be a solution to this problem and u k (t) = (u(t), v k ). Then, by virtue of Equation (53) and the self-adjointness of operator A, Thus, taking into account (54) and (55), we have the following problem Suppose that f k is known and use the nonlocal condition to obtain (see, e.g., [20], p. 174) Note, that if α ≥ 1 and t > 0 then, α − E ρ (−λ k t ρ ) > 0, E ρ (−λ k t ρ ) > 0 and E ρ,ρ+1 (−λ k t ρ ) > 0. Therefore, for all k one has f k = 0. Hence, from the completeness of the system of eigenfunctions {v k }, we finally obtain f = 0 and u(t) ≡ 0, as required. The uniqueness is proved.
In Theorem 3, condition E ρ (−λ k ξ ρ ) = α for all k, ensured the uniqueness of the solution and no orthogonality conditions were required for its existence. In this inverse problem, the same condition does not guarantee the uniqueness of the solution. The reason is that at α < 1 the bracket in (57) can vanish for some k 0 , and these numbers k 0 depend on α and the location of points ξ and τ. For the existence of an unknown function f in this case, it is necessary to require the orthogonality conditions (ϕ, v k 0 ) = 0 and (Ψ, v k 0 ) = 0. Finding out the dependence of k 0 on these parameters is the subject of a separate paper. This is the reason why only case α ≥ 1 is considered in Theorem 4.
On the other hand, in the next two theorems the condition E ρ (−λ k ξ ρ ) = α for all k ensures the unique solvability of the problems under study. If E ρ (−λ k ξ ρ ) = α for some k, then, just as above (see the proof of Theorem 3), there is no uniqueness, and to ensure the existence of the solution, it is necessary to require the fulfillment of certain orthogonality conditions.

The Inverse Problem of Determining the Boundary Function ϕ
Consider the problem (2) and assume that, together with function u(t), function ϕ in the nonlocal condition u(ξ) = αu(0) + ϕ is also unknown. To solve this inverse problem, we need an additional condition, and as such we again take the condition that was used in the previous inverse problem: If τ = ξ, then the nonlocal condition u(ξ) = αu(0) + ϕ coincides with the Cauchy condition u(0) = ϕ 1 (see (4)) and we have the inverse problem, considered in [2][3][4]. As noted at the end of Section 4, an additional condition is imposed on α, which simplifies the proof of the uniqueness of the solution.
Theorem 5. Let Ψ ∈ D(A), f ∈ C([0, T]; D(A ε )) for some ε ∈ (0, 1) and E ρ (−λ k ξ ρ ) = α for all k. Then the inverse problem (2), (58) has a unique solution {u(t), ϕ} and this solution has the form where Proof of Theorem 5. Existence. If ϕ is known, then the solution of problem (3) has the form (60) (see Theorem 3). Condition (58) implies: Let us expand function Ψ ∈ H in a Fourier series in system {v k }. Then Therefore, equality (59) is formally established. It remains to prove that ϕ ∈ H. Let Φ j be the partial sum of series (59). Then due to Parseval's equality Since |E ρ (−λ k ξ ρ ) − α| ≤ C, then by virtue of the asymptotic estimate (8), we obtain Similarly, by virtue of the estimate (11) and the definition of ω k , we have For the sum Φ 3 j one has Thus, it is shown that ϕ ∈ H (see (61)). Fulfillment of the conditions of Definition 3 for function u(t), defined by the series (60) is proved in exactly the same way as Theorem 3.
Uniqueness. Let us prove that if {u(t), ϕ} is a solution to the homogeneous problem: then u(t) ≡ 0 and ϕ = 0. Let u(t) be a solution to this problem and let u k (t) = (u(t), v k ). Then The solution to this problem has the form (this is the same problem as (23)) Since E ρ (−λ k ξ ρ ) = α for all k ≥ 1, then by virtue of the properties of the Mittag-Leffler functions (see Proposition 2) we have ϕ k = 0 for all k. This, in turn, means u k (t) ≡ 0 for all k. Therefore, due to the completeness of the system of eigenfunctions {v k }, we have ϕ = 0 and u(t) ≡ 0, as required.
Theorem 5 is completely proved.

Remark 1.
If f does not depend on t, then the statement of Theorem 5 is true for all f ∈ H, and Ψ ∈ D(A).
Using Corollary 1, this statement can be proved in the same way as above.

Well-Posedness of Problem (3)
In the case of fractional Riemann-Liouville derivatives, we consider only the forward problem for the homogeneous subdiffusion equation. The inhomogeneous equations and inverse problems considered above are studied in exactly the same way as in the case of the Caputo derivatives.
Consider the forward problem: where φ ∈ H and the number α are given. As was mentioned above (see Section 4) we will only consider the case E ρ (−λ k ξ ρ ) = α for all k ≥ 1.

Theorem 6.
Let φ ∈ H and E ρ (−λ k ξ ρ ) = α for all k ≥ 1. Then problem (65) has a unique solution and this solution has the form Moreover, there is a constant C ξ > 0 such that the following coercive type inequality holds: Proof. As in the case of problem (2), we will seek a solution to problem (65) in the form of a series: ∑ ∞ k=1 T k (t)v k , where T k (t) is a solution to the problem: Let us denote lim t→0 ∂ ρ−1 t T k (t) = b k . Then the unique solution to the differential Equation (68) with this initial condition has the form T k (t) = b k t ρ−1 E ρ,ρ (−λ k t ρ ) (see, e.g., [20], p. 173, [52]). Proposition 5 implies Then from the nonlocal conditions (69) we find the unknown numbers b k : Hence, function u(t) defined by series (66) is a formal solution to problem (65).

Now let us show that series (66) is indeed a solution.
To do this we denote by S j (t) the partial sum of series (66). Then Due to the Parseval equality we may write Estimate (10) of function E ρ,ρ (−t) implies (note, E ρ (−λ k ξ ρ ) = α) Thus A 2 u(t) ∈ C((0, T]; H) (and even more so Au(t) ∈ C((0, T]; H)).
Since  Thus, (66) is a solution to problem (65). The coercivity inequality follows from the last two estimates.
The uniqueness is proved in the same way as in the proof of Theorem 2.

Remark 2.
Having considered problem (65), we would like to draw the reader's attention to the fact that although in the case of the Riemann-Liouville derivatives function T k (t) is determined using Mittag-Leffler function with two parameters, but thanks to Proposition 5, the coefficients b k are determined in exactly the same way as in the case of the Caputo derivatives. It is this fact that allows us to apply the same reasoning as above to solve the problem for an inhomogeneous equation with Riemann-Liouville derivatives and inverse problems.

Conclusions
We have considered the following three problems: (1) the nonlocal in time forward problems (2) and (3); (2) the inverse source problem (46) with additional condition (47); (3) the inverse problem of determining the boundary function ϕ. To the best of the authors knowledge, forward problems (2) and (3) with such nonlocal conditions for differential equations with fractional derivatives have not been discussed before. In the study of forward problems, we determined the influence of parameter α on the correctness of the problems under study. It should be noted, as proved above, for some values of α there is no uniqueness of the solution, and to ensure the existence of a solution, it is necessary to require orthogonality conditions for the given functions ϕ and f . From the proofs of Theorems 4 and 5, it can be seen that the same values of the parameter α violate the uniqueness of solutions to inverse problems.
In [45] nonlocal in time (forward and space-dependent inverse source) problems for a so-called generalized time fractional diffusion equation were considered. The fractional derivative used in [45] is the Hilfer derivative, which is a generalization of the well-known Riemann-Liouville and Caputo fractional derivatives. Indeed, the generalized fractional derivative interpolates the Riemann-Liouville and Caputo fractional derivatives. The results presented in this paper can be extended by considering the Hilfer derivative.