Fractional Cauchy problem on random snowflakes

We consider time-changed Brownian motions on random Koch (pre-fractal and fractal) domains where the time change is given by the inverse to a subordinator. In particular, we study the fractional Cauchy problem with Robin condition on the pre-fractal boundary obtaining asymptotic results for the corresponding fractional diffusions with Robin, Neumann and Dirichlet boundary conditions on the fractal domain.


Introduction
Many physical and biological phenomena take place across irregular and wild structures in which boundaries are "large"while bulk is "small". In this framework, domains with fractal boundaries provide a suitable setting to model phenomena in which the surface effects are enhanced like, for example, pulmonary system, root infiltration, tree foliage, etc..
In this paper, we consider random Koch domains which are domains whose boundary are constructed by mixtures of Koch curves with random scales. These domains are obtained as limit of domains with Lipschitz boundary whereas for the limit object, the fractal given by the random Koch domain, the boundary has Hausdorff dimension between 1 and 2.
Our attention will be focused on fractional Cauchy problems on the random Koch domains with boundary conditions. Literature on fractional Cauchy problems is extensive both from the probability and the analysis point of view. Here, our aim is not providing a large list of references. We mention here only few works investigating basic and fondamental aspects: [1], [3], [14], [17], [18], [20], [22], [28], [32].
The non-local time-operator we deal with is very general and covers a huge class of nonlocal (convolution type) operators. Such operators have been recently considered in the papers [13; 31]. From the probabilistic point of view, we consider time-changed Brownian motions where the time change is given by an inverse to a subordinator characterized by a symbol which is a Bernstein function. Thus, with this time-fractional operator at hand, we study the fractional Cauchy problem with Robin condition on the pre-fractal boundary and we obtain asymptotic results for the corresponding fractional diffusions with Robin, Neumann and Dirichlet boundary conditions on the fractal domain.
The asymptotic problem we deal with can be illustrated, in the simple case, by the following parabolic Dirichlet-Robin problem on the interval (0, a), a > 0. More precisely, we consider the heat equation ∂ t u n = ∂ xx u n , t > 0, x ∈ (0, a) u n (t, 0) = 0, t > 0 (1.1) with Robin boundary condition ∂ x u n (t, a) + c n u n (t, a) = 0, t > 0 (1.2) where c n > 0, n ∈ N. The solution can be written as follows u n (t, x) = k≥1 e −tλ (n) k φ (n) Our aim here is to point out the asymptotic behaviour of the solution u n as n → ∞. We obtain three different limit problems. If c n → 0, then z Now we wonder if a similar asymptotic behaviour holds for the analogue time-fractional problem. A simple example is given by the problem with c n > 0, n ∈ N where ∂ β t u is the Caputo fractional derivative of u (see formula (4.5) below). The solution can be written as follows , w ≥ 0 is the Mittag-Leffler function and the system {φ k : k ∈ N} has been introduced before. By simple arguments, we get that the solution (1.3) uniformly converges to a function u which turns out to be analogously related to the boundary problems above (Neumann, Dirichlet, Robin) with the Caputo time-fractional derivative ∂ β t in place of the ordinary derivative ∂ t . This is due to the fact that we have explicit representation of the system {φ Following the same spirit, in the present paper, we move on to general domains like the random snowflakes we have introduced before and we address the same asymptotic problem with a general time-fractional operator. In this case we do not have the same informations about the associated system and the compact representation of the solution. We overcome this difficulty by using the theory of Dirichlet forms and Markov processes. An essential tool will be given by the convergence of forms associated with time-changed processes.
We remark that the peculiarity in studying the asymptotic behaviour of these approximating problems is that one has to deal with an increasing sequence of Lipschitzian domains which converges in the limit to the domain whose boundary is a fractal.
The plan of the paper is the following: in Section 2 we introduce the random Koch domains; in Section 3 we recall the definition of Dirichlet forms with associated base processes; in Section 4 we introduce time-fractional equations and time changes; in the last section we prove our main results. More precisely, in Theorem 5.1 we solve the asymptotic problem for the the timechanged processes and in Theorem 5.2 we point out some peculiar aspects arising by passing from the ordinary to the fractional Cauchy problem.
The fractal K (ξ) associated with the random environment sequence ξ is therefore defined by Let Ω (ξ|n) be the planar domain obtained from a regular polygon by replacing each side with a pre-fractal curve K (ξ) n and Ω (ξ) be the planar domain obtained by replacing each side with the corresponding fractal curve K (ξ) . We introduce the random planar domains Ω (ξ|n) and Ω (ξ) by considering the random curves K (ξ) n and K (ξ) . Examples of (pre-fractal) random Koch domains are given in figures 1 (outward curves), 2 (inward curves), 3 (inward curves) by choosing as regular polygon the square.
Since ξ i law = ξ 1 , ∀ i, we have that the Hausdorff dimension d (ξ) of the curve K (ξ) can be obtained by considering the strong law of large numbers and the fact that , n → ∞.
Moreover the measure µ (ξ) in (2.3) has the property that there exist two positive constants C 1 , C 2 , such that, where B(P, r) denotes the Euclidean ball with center in P and radius 0 < r ≤ 1 (see [2]). According to Jonsson and Wallin (see [25]), we say that K (ξ) is a d-set with respect to the Hausdorff measure H d , with d = d (ξ) . The sequence is obtained from the realization of ξ|n and therefore, from the realization of the random variable (ξ|n) with mean value given by Thus, for α = E ξ 1 ∈ (2, 4) we find the mean value E[σ (ξ|n) ] = α n /4 n . The realization ξ|n can be regarded as the vector a|n = (a 1 , . . . , a n ) which is a n-dimensional vector with N different values of I, that is a|n ∈ I n . We introduce the multinomial distribution ..,N . Thus, for the realization of the vector ξ|n we have that or equivalently We notice that

Dirichlet forms and Base processes
Let E be a locally compact, separable metric space and The point ∂ is the cemetery point for X and a function f on E can be extended to E ∂ by setting f (∂) = 0. The associated semigroup is uniquely defined by be the Dirichlet form associated with (the non-positive definite, self-adjoint operator) A. Then X is equivalent to an m-symmetric Hunt process whose Dirichlet form (E, D(E)) is on L 2 (E) (see the books [15; 23]). Without restrictions we assume that the form is regular ([23, page 143]).
We say that X is the base process. Our aim is to consider time changes of the base process X. Such random times will be introduced in the next section.

Time fractional equations and Time changes
We first introduce the subordinator H = {H t , t ≥ 0} for which where Φ is the symbol of H. The symbol Φ may be associated also to the inverse L of H, that is L = {L t , t ≥ 0} defined as L t = inf{s ≥ 0 : H s > t}, t ≥ 0.
We assume that H 0 = 0, L 0 = 0. By definition, we also have that P 0 (H t < s) = P 0 (L s > t), s, t > 0. (4.1) The symbol Φ we consider hereafter is a Bernstein function with representation and Π is the so called tail of the Lévy measure. Both random times H, L are non-decreasing. We do not consider step-processes with Π((0, ∞)) < ∞ and therefore we focus only on strictly increasing subordinators with infinite measures. Thus, the inverse process L turns out to be a continuous process. For details, see the books [5; 30].
We now introduce the fractional operators and the fractional equations governing the timechanged process X L = {X • L t , t ≥ 0}, that is the base process X = {X t , t ≥ 0} with the time change L characterized by the symbol Φ.
Let M > 0 and w ≥ 0. Let M w be the set of (piecewise) continuous function on [0, ∞) of exponential order w such that |u(t)| ≤ M e wt . Denote by u the Laplace transform of u. Then, we define the operator where Φ is given in (4.2). Since u is exponentially bounded, the integral u is absolutely convergent for λ > w. By Lerch's theorem the inverse Laplace transforms u and D Φ t u are uniquely defined. Notice that Simple arguments say that D Φ t can be written as a convolution involving the ordinary derivative and the inverse transform of (4.3) iff u ∈ M w ∩ C([0, ∞), R + ) and u ∈ M w , that is, We notice that when Φ(λ) = λ (that is, the ordinary derivative) we have that a.s. H t = t and L t = t. We also notice that for Φ(λ) = λ β , the symbol of a stable subordinator, the operator D Φ t becomes the Caputo fractional derivative with u (s) = du/ds. For Φ(λ) = (λ + η) β − η β , with η ≥ 0 and β ∈ (0, 1), the operator D Φ t becomes the Caputo tempered fractional derivative For explicit representation of the operator D Φ t see also the recent works [13; 31].
Let X be the process with generator (A, D(A)) introduced above. In the present work we consider the time fractional equation The probabilistic representation of the solution to (4.6) is written in terms of the time-changed process X L , that is We notice that (4.7) is not a semigroup, indeed the random time L is not Markovian and therefore, the composition X L is not a Markov process. The fractional Cauchy problem has been investigated by many authors by considering Caputo derivative and only recently, by taking into account more general operators. The following theorem has been obtained in [11] for Feller processes (not necessarily Feller diffusions, see [11]) and we mention here such a result for the reader's convenience. (1) ϕ : t → u(t, ·) is such that ϕ ∈ C([0, ∞), R + ) and ϕ ∈ M 0 , (2) ϑ : x → u(·, x) is such that ϑ, Aϑ ∈ D(A), In [13] the author proves existence and uniqueness of strong solutions to general time fractional equations with initial datum f ∈ D(A). In [16] the authors establish existence and uniqueness for weak solutions and initial datum f ∈ L 2 . The result in Theorem 4.1 has been proved in a general setting, that is by considering a generator of a Feller process as in [13] but following a very different approach. We notice that the condition on the initial datum f must be better specified for the compact representation of the solution, this is the case investigated in [17] for instance (the domain has no boundary) or the case investigated in [14] (with Dirichlet condition on the boundary).
In the next section we will study continuous base processes time changed by continuous random times, thus we do not stress the fact that the previous result holds for Feller process (right-continuous with no discontinuity other than jumps).

Main results
We consider the prefractal RKD Ω (ξ|n) defined in Section 2 and we construct the set Ω (ξ|n) \ B where B ⊂ Ω (ξ|1) is a ball.
Then, we consider Brownian diffusions on the Random Koch Domain Ω (ξ|n) \ B. Let X n = {X n t , t ≥ 0} with X n 0 = x ∈ Ω (ξ|n) \ B be a sequence of planar Brownian motions for a given ξ ∈ Ξ. Let (A n , D(A n )) be the generator of X n , in particular A n = ∆ and where c n ≥ 0, n(x) denote the inward normal vector at x ∈ ∂Ω (ξ|n) and σ (ξ|n) is defined in (2.6). It is well-known that there is one to one correspondence between the infinitesimal generator of X n and the closed symmetric form (E n , D(E n ) (see [23,Theorem 1.3

.1]).
We recall that a form (E n , D(E n ) can be defined in the whole of L 2 (F, m) by setting E n (u, u) = +∞ ∀ u ∈ L 2 (F, m) \ D(E n ). Similarly a forms E, E can be defined in the whole of L 2 (F, m) by setting E(u, u) = +∞ ∀ u ∈ L 2 (F, m) \ D(E).
For the convenience of the readers we recall the definition of convergence of forms introduced by Mosco in [27], denoted by M -convergence. In our framework, we consider the pre-fractal form E n (·, ·) on L 2 (Ω (ξ) \ B) by defining We now introduce the time-changed process X L,n = X n • L and we study the asymptotic behaviour of X L,n depending on the asymptotics for c n . The process X L,n can be considered in order to study the corresponding time-fractional Cauchy problem on Ω (ξ|n) \ B   [21]. Random time change theorem). Suppose that X n , X are in D and L n , L are in D 0 . If (X n , L n ) converges to (X, L) in distribution as n → ∞, then X n • L n converges to X • L in distribution as n → ∞.
Proof. The proof follows from part b) of Theorem 1.1 and part a) of Lemma 2.3 in [21]. Lemma 2.3 gives convergence for strictly increasing time changes. Since H is strictly increasing, we use part c) of Theorem 1.1 and find results for L which is non-decreasing and continuous. Then, part b) holds for the random time changes L n . Theorem 5.1. As n → ∞, In particular, as c n → c ≥ 0, i) if c = 0, then X L is reflected on ∂Ω (ξ) , that is the process driven by ii) if c ∈ (0, ∞), then X L is (elastic) partially reflected on ∂Ω (ξ) , that is the process driven by ii) if c = ∞, then X L is killed on ∂Ω (ξ) , that is the process driven by Remark 5.1. We point out that the condition on the boundary ∂Ω (ξ) must be meant in the dual of certain Besov spaces (for details, see [8], [26] and the references therein).
Proof. Fix ξ ∈ Ξ. First we prove the M-convergence in L 2 (Ω (ξ) \ B) of the Dirichlet forms E n .
The case of finite limit has been addressed in Theorem 5.2 in [9]: in particular, it has been proved that if c n → c ≥ 0, then the sequence of forms E n (·, ·) M -converges in the space The last form E c , for c ∈ (0, ∞), is associated with the semigroup ( [6; 15]) where the multiplicative functional M t is associated to the Revuz measure given by the perturbation of the form E c . Thus, (5.3) is the solution to ∂ t u = ∆ R u, u 0 = f ∈ D(∆ R ). For c = 0, the form E c is associated with solution to ∂ t u = ∆ N u, u 0 = f ∈ D(A). Now we prove that if c n → ∞, the sequence of forms E n M-converges on L 2 (Ω (ξ) ) to the form First we prove condition (a) of Definition 1. Up to passing to a subsequence, which we still denote by v n , we can suppose that v n | Ω (ξ|n) \B ∈ H 1 (Ω (ξ|n) \ B), (5.4) and, for every n, ||v n || H 1 (Ω (ξ|n) \B) c * , (5.5) with c * independent of n. First we extend v n by Jones extension operator (Theorem 1 in [24]) and after we restrict it to the domain Ω (ξ) \ B : more precisely, we extend v n to a function v * We point out that the constant C J independent of n (see Theorem 3.4 in [9]) that is the norm of extension operator is independent of the (increasing) number of sides.
Then, there exists v * such that the sequence v * n weakly converges to v * in H 1 (Ω (ξ) \ B) : for the uniqueness of the limit in the weak topology, we obtain that v * = u and, in particular, u ∈ H 1 (Ω (ξ) \ B). Since the sequence v * n weakly converges to u in H 1 (Ω (ξ) \ B), we have that From the compact embedding of and by using Trace theorems (see [25] and [10]) we obtain that when n → ∞ (see Theorem 2.1 in [9]). We stress the fact that the value of σ (ξ|n) play a crucial role in the previous limit. Now, if c n → ∞, for any k > 0 there exists n 1 such that, for all n > n 1 , c n ≥ k. Then c n σ (ξ|n) when n → ∞. Dividing for k and letting k → ∞ we obtain that ∂Ω (ξ) |u| 2 dµ (ξ) = 0 (5.11) and so u = 0 on ∂Ω (ξ) . By combining (5.7), (5.9), (5.11) we have proved condition (a) of Definition 1.
In order to prove condition (b) of Definition 1, we can assume that u ∈ H 1 0 (Ω (ξ) \ B) without loss of generality: then, the choice of v n = u suffices to achieve the result. So we have proved the M-convergence of the forms E n (·, ·) on L 2 (Ω (ξ) \ B) to the form E ∞ when c n → ∞.
From the M-convergence of the forms E n (·, ·) on L 2 (Ω (ξ) \B), by using the results in the recent paper [11], we obtain the convergence of the time changed processes.
More precisely, from the M -convergence of the forms we have the strong convergence of semigroups. From Theorem 17.25 (Trotter, Sova, Kurtz, Mackevičius) in [19] we have that strong convergence of semigroups (Feller semigroups) is equivalent to weak convergence of measures if X n 0 → X 0 in distribution. Then we obtain that X n d → X in D. From Proposition 5.1, we have that ∀ ξ ∈ Ξ, X n • L =: X L,n → X L := X • L on Ω (ξ) in distribution as n → ∞ in D.
From the pointwise convergence, we get that ξ−a.s.
in distribution as n → ∞ in D.
Let us consider now the process X L on E. We point out some peculiar aspects of X L and the corresponding lifetimes.
Theorem 5.2. Let us consider the Cauchy problems and We have that, ∀ x ∈ E: Proof. The solution to (5.12) has the following probabilistic representation where M t = 1 (t<ζ) is the multiplicative functional written in terms of the lifetime ζ of the process X on E. Then, we consider the part process X of * X where * X 0 = x ∈ E. It is well known that M t characterizes uniquely the associated semigroup ( [6]), that is the solution w. We also have that is the solution to the elliptic problem on E −Aw = f.
From Theorem 4.1 we have that the time-changed process X L can be considered in order to solve the problem (5.13), that is where ζ L is the lifetime of X L . As before we introduce the which is the solution to the elliptic problem associated with the fractional Cauchy problem (5.13). We are able to obtain the key relation between w and u by taking into consideration the following plain calculations. First we recall (4.7) where P s f (x) here is given by w(s, x) with w(s, x) → f (x) as s → 0. Moreover (see [12]), ∞ 0 e −λt P 0 (L t ∈ ds)dt = Φ(λ) λ e −sΦ(λ) ds. (5.14) We have that u(x) = lim That is and this gives a connection between solutions of elliptic problems introduced above in the proof. Since Φ is a Bernstein function with Φ(0) = 0 we get the result.
The characterization given in the previous result admits a probabilistic interpretation in terms of mean lifetime of the base and time-changed processes. The problems (5.12) and (5.13) with f = 1 E are associated with w(x) = E x [ζ] and u(x) = E x [ζ L ] as described in the previous proof and the mean lifetime says how much the time change L modifies the base process X. By following the definition given in [12] and the relation between w and u we say that X is delayed or rushed on E by L. An example is given by the tempered fractional derivative ([4; 29]) associated with the symbol Φ(λ) = (λ + η) β − η β with η > 0 and β ∈ (0, 1). We get that that is, if βη β−1 < 1 then the process X is rushed by L, whereas if βη β−1 > 1 then the process X is delayed by L.
The previous discussion on either delayed or rushed processes holds according to specific regularity conditions on the boundary ∂E. We must have that sup E w(x) < ∞ which is the characterization of trap domains (written here for X with generator A) given in [7] for the Brownian motion. By applying the result in [7] it follows that the following proposition holds true.