Differentiation of Solutions of Caputo Boundary Value Problems with Respect to Boundary Data

: Under suitable continuity and uniqueness conditions, solutions of an α order Caputo fractional boundary value problem are differentiated with respect to boundary values and boundary points. This extends well-known results for nth order boundary value problems. The approach used applies a standard algorithm to achieve the result and makes heavy use of recent results for differentiation of solutions of Caputo fractional intial value problems with respect to initial conditions and continuous dependence for Caputo fractional boundary value problems.


Introduction
Let n ∈ N with α ∈ (n − 1, n) and a < t 0 < b in R. Our concern is characterizing partial derivatives with respect to the boundary data for solutions to the Caputo fractional boundary value problem D α * t 0 x(t) = f (t, x(t), x ′ (t), . . ., x (n−1) (t)), a < t 0 < t < b, satisfying conjugate boundary conditions where D α * t 0 x is the Caputo fractional derivative of order α of the function x(t) and a < t 0 ≤ t 1 < t 2 < . . .< t n < b and x i ∈ R for 1 ≤ i ≤ n.These partial derivatives solve the associated Caputo fractional variational equation.∂ f ∂x j (t, x(t), x ′ (t), . . ., x (n−1) (t))z (j) . ( In this paper, we impose suitable continuity and uniqueness hypotheses so that given a solution of (1), (2), one may take the derivative with respect to the boundary data.This derivative solves the variational Equation (3) with interesting boundary data where all but one of the boundary values are zero.Colloquially, we refer to this as studying the smoothness of conditions.
The history of initial and boundary data smoothness dates back to Peano and his work on the smoothness of initial conditions as cited by Hartman [1].Subsequently, Peterson [2], Spencer [3], and Sukup [4] were among the first to shift to studying the smoothness of boundary conditions.In the following decade, these results were then extended by Henderson to right-focal boundary conditions [5,6] and conjugate-type boundary conditions [7].Over the next several decades, results were introduced for nonlocal boundary conditions [8][9][10], difference equations [11][12][13], dynamic equations on time scales [14][15][16], and researchers incorporated parameters into the nonlinearity [15,17].
With this work, we broaden the scope even further by analyzing smoothness of solutions to Caputo fractional boundary value problems.Research into fractional differential equations has seen an explosion of articles in the past decade that seek to generalize results for integer order differential equations to fractional order.To name a few, we cite [18][19][20][21][22][23][24][25][26][27].In fact, there also seem to be a limitless number of different ways to define a fractional derivative.However, two definitions have become the source of focus amongst a broad range of researchers in the field; namely the Riemann-Liouville and Caputo fractional derivatives.Brief definitions may be found in Section 2. For expository material on fractional differential equations, we refer the reader to [28][29][30][31].
The theorems and proof in this article are novel as no other research to date has attempted to extend boundary data smoothness to fractional differential equations.The reason is that the results found in this article rely heavily upon two recent results for Caputo fractional differential equations.The first establishes differentiation of solutions of Caputo initial value problems with respect to the initial data [22], and the second establishes the continuous dependence on boundary conditions for Caputo boundary value problems [32].
The idea behind the proof of our main result is to first assume a unique solution to a Caputo boundary value problem.Then, we define a difference quotient with respect to the boundary point or boundary value of interest.We view this difference quotient in terms of an initial value problem and apply Theorem 3.2 from [22].This yields that the difference quotient solves the variational equation.Finally, we take a limit by applying the continuous dependence result, Theorem 4.2, from [32] which yields the desired result.
The remainder of the paper is organized as follows.In Section 2, one will find brief definitions of fractional integrals and derivatives.Section 3 is where we establish our sufficient hypotheses.For Section 4, we present important recent results in continuous dependence and smoothness of initial conditions.Following this, we have Section 5 that contains the main result and its proof.Finally, we conclude with a summary of project and thoughts on future research avenues.

Fractional Derivatives
Let α > 0. The Riemann-Liouville fractional integral of a function x of order α, denoted I α t 0 x, is defined as provided the right-hand side exists.Moreover, let n ∈ N denote a positive integer and assume n − 1 < α ≤ n.The Riemann-Liouville fractional derivative of order α of the function x, denoted D α t 0 x, is defined as provided the right-hand side exists.If a function x is such that exists, then the Caputo fractional derivative of order α of x is defined by Remark 1.A sufficient condition to guarantee the existence of the Caputo fractional derivative is the absolute continuity of the (n − 1)st derivative of x(t).See Theorem 3.1 in [28] and discussion thereafter.

Preliminaries
Throughout this work, we make use of the following assumptions which are required to apply the continuous dependence and differentiation results from [22,32] (1) (3) solutions to initial value problems for (1) are unique on (a, b); Next, we present two more hypotheses which establish a uniqueness condition for (1) and (3), respectively.
(4) Given points a < t 0 ≤ t 1 < t 2 < . . .< t n < b, if y and z are solutions of (1) such that for 1 Next, we present two crucial results that make this work possible.
We seek a boundary value problem result as an analog of the initial value problem result from Eloe et al [22].
Theorem 1. Assume that hypotheses (1)-(3) hold.Let y(t) := y(t; t 0 , y 0 , . . ., y n−1 ) be the unique solution of the initial value problem (1) satisfying with maximal interval of existence [t 0 , ω).Choose [c, d] ⊂ [t 0 , ω).Then, (a) for each 0 ≤ j ≤ n − 1, γ j (t) := ∂y(t)/∂y j exists and is the solution of the variational Equation (3) along y(t) on [c, d] and hence, [t 0 , ω) satisfying the initial conditions (b) if, in addition, f has a continuous first derivative with respect to t and f (t 0 , y 0 , y 1 , . . ., y n−1 ) = 0, then β(t) := ∂y(t)/∂t 0 exists and is the solution of the variational Equation (3) along y(t) on [c, d] and hence, [t 0 , ω) satisfying the initial conditions We also use recent continuous dependence on boundary conditions results for Caputo fractional differential equations [32].The first one is if the left-most boundary condition is to the right of the starting point of the Caputo fractional derivative; namely t 0 < t 1 , and the second is if they are equal; namely t 0 = t 1 .Note that the second result has an additional condition to establish continuous dependence to the left of t 0 .

Main Results
In this section, we present our boundary value problem analog.First, we state and prove the result when t 0 < t 1 .
Proof.We will only prove part (a) as the proof of part (b) is similar.Part (c) is an immediate consequence from parts (a) and (b) when coupled with hypothesis (5).Let 1 ≤ j ≤ n, and consider ∂x(t)/∂x j .In the interests of conserving space and lessening the tedious notation, we denote x(t; t 1 , . . ., t n , x 1 , . . ., x j , . . ., x n ) by x(t; x j ) as x j is the boundary value of interest.
Let δ > 0 be as in Theorem 2, 0 < |h| < δ be given, and define the difference quotient with respect to x j by Note that for every h ̸ = 0, Also, for every Now that we have established the boundary conditions for z jh (t), we show that z jh (t) solves the variational equation.To that end, for 1 ≤ i ≤ n − 1, let By Theorem 2, for 1 Using the notation of Theorem 1 for solutions of initial value problems for (1), viewing x(t) as the solution of an initial value problem, and denoting the solution x(t) by y(t; t j , x j , v 1 , . . ., v n−1 ), we have Then, by utilizing telescoping sums, we have By Theorem 1 and the Mean Value Theorem, we obtain where, for 0 is the solution of the variational Equation (3) along y(•) satisfying is between v i and v i + ϵ i .Thus, to show lim h→0 z jh (x) exists, it suffices to show, for each 1 ≤ i ≤ n − 1, lim h→0 ϵ i /h exists.Now, from the construction of z jh (t), we have Hence, for 1 ≤ k ≤ n with k ̸ = j, we have a system of n − 1 linear equations with n − 1 unknowns: γ n−1 (t k , y(t; t j , x j , v 1 , . . ., v n−1 + εn−1 )).
In the system of equations above, we notice that y(•) is not always the same.Therefore, we consider the coefficient matrix M based on y(t) We claim det(M) ̸ = 0. Suppose to the contrary that det(M) = 0.Then, there exist (c) for each 2 ≤ j ≤ n, if f has a continuous first derivative with respect to t and f t j , x(t j ), x ′ (t j ), . . ., x (n−1) (t j ) = 0, then w j (t) := ∂x(t)/∂t j exists and is the solution of the variational Equation (3) along x(t) on [c, d] and hence, [t 0 , ω) satisfying the boundary conditions (d) Under the conditions of (b) and (c), for each 1 ≤ j ≤ n, w j (t) = −x ′ (t j )z j (t).

Conclusions
In this paper, we showed that under suitable continuity and uniqueness conditions that a solution Caputo fractional conjugate boundary value problem may be differentiated with respect to the boundary points and the boundary values.The resulting function solves the Caputo fractional version of the variational equation.This work only recently became possible as its proof relies extensively upon the differentiation of a solution to a Caputo fractional initial value problem [22] and the continuous dependence of solutions to Caputo fractional boundary value problems with respect to boundary data [32].
The results contained herein are novel and have not been explored or considered previously.We believe this result is foundational in smoothness of solutions for Caputo fractional boundary value problems and the proof sets a template for how to proceed for several future research avenues such as Caputo fractional differential equations with varying types of boundary conditions including parameter dependence, Caputo fractional difference equations, and Caputo fractional dynamic equations.
When looking at the wide breadth and depth of research conducted for smoothness of solutions for integer order differential, difference, and dynamic equations, it is clear there is a lot of work to be done in this area.
Other future work could entail loosening the hypothesis that the nonlinearity be continuously differentiable.This was posited for future study in [22].Another avenue would be finding sufficient conditions to guarantee the uniqueness condition in hypothesis (3) in certain contexts.