2000 Mathematics Subject Classification. Primary 60H20; Secondary 60H05, 45D05. The second author is supported in part by FONDECYT Grant #1010675 .
<ph f="cmbx">Regularity of solutions to stochastic Volterra equations with infinite delay</ph>

Anna Karczewska

Carlos Lizama

Department of Mathematics, University of Zielona Gora, ul. Szafrana 4a, 65-246 Zielona Gora, Poland. E-mail address : A.Karczewska@im.uz.zgora.pl Universidad de Santiago de Chile, Departamento de Matematica, Facultad de Ciencias, Casilla 307-Correo 2, Santiago, Chile. E-mail address : clizama@lauca.usach.cl

1 Introduction

We consider the following integro-differential stochastic equation with infinite delay
X ( t , θ ) = t b ( t s ) [ Δ X ( s , θ ) + W Γ s ( s , θ ) ] d s , t 0 , θ T d , (1.1)
where b L 1 ( R + )   , Δ   is the Laplace operator and T d   is the d   -dimensional torus.
In ( 1.1 ), W Γ   is a spatially homogeneous Wiener process with the space covariance Γ   and taking values in the space of tempered distributions S ( T d ) .   Such equation arises, in the deterministic case, in the study of heat flow in materials of fading memory type (see [3, [14).
In this paper we address the following question: under what conditions on the covariance Γ   the process X   takes values in a Sobolev space H α ( T d )   , particularly in L 2 ( T d )   ?
We study a particular case of weak solutions under the basis of an explicit representation of the solution to ( 1.1 ) (cf. Definition  3.4 ).
Observe that equation ( 1.1 ) can be viewed (see [16) as the limiting equation for the stochastic Volterra equation
X ( t , θ ) = 0 t b ( t s ) [ Δ X ( s , θ ) + W Γ s ( s , θ ) ] d s , t 0 , θ T d . (1.2)
If b   is sufficiently regular, we get, by differentiating ( 1.2 ) with respect to t   ,
X t ( t , θ ) = b ( 0 ) [ Δ X ( t , θ ) + W Γ t ( t , θ ) ] + 0 t b ( t s ) [ Δ X ( s , θ ) + W Γ s ( s , θ ) ] d s , (1.3)
where t 0   and θ T d .   Taking in ( 1.3 ), b ( t ) 1   we obtain
{ X t ( t , θ ) = Δ X ( t , θ ) + W Γ t ( t , θ ) , t > 0 , θ T d , X ( 0 , θ ) = 0 , θ T d . (1.4)
Similarly, taking b ( t ) t   and differentiating ( 1.2 ) twice with respect to t   we obtain
{ 2 X t 2 ( t , θ ) = Δ X ( t , θ ) + W Γ t ( t , θ ) , t > 0 , θ T d , X ( 0 , θ ) = 0 , θ T d , X t ( 0 , θ ) = 0 , θ T d . (1.5)
It has been shown in [8,Theorem5.1(see also [7,Theorem1) that equations ( 1.4 ) and ( 1.5 ) on the d   -dimensional torus T d   have an H α + 1 ( T d )   -valued solutions if and only if the Fourier coefficients ( γ n )   of the space covariance Γ   of the process W Γ   satisfy
n Z d γ n ( 1 + | n | 2 ) α < . (1.6)
Observe that for both, stochastic heat ( 1.4 ) and wave ( 1.5 ) equations, the conditions are exactly the same, despite of the different nature of the equations. On the other hand, the obtained characterization form a natural framework in which nonlinear heat and wave equations can be studied.
In the present paper, we will prove that condition ( 1.6 ) even characterizes H α + 1 ( T d )   -valued solutions for the stochastic Volterra equation ( 1.1 ), provided certain conditions on the kernel b   are satisfied. This is a strong contrast with the deterministic case, where regularity of ( 1.1 ) is dependent on the kernel b   .
We use, instead of resolvent families, a direct approach to the equation ( 1.1 ) finding an explicit expression for the solution in terms of the kernel b   . This approach reduces the considered problem to questions in harmonic analysis and lead us with a complete answer.

2 Preliminaries

Let ( Ω , , ( t ) t 0 , P )   be a complete filtered probability space. By T d   we denote the d   -dimensional torus which can be identified with the product ( π , π ] d   . Let D ( T d )   and D ( T d )   denote, respectively, the space of test functions on T d   and the space of distributions. By ξ , φ   we denote the value of a distribution ξ   on a test function. We assume that W Γ   is a D ( T d )   -valued spatially homogeneous Wiener process (abbr. s.h. Wiener process) with covariance Γ   which is a positive definite distribution. For more information about s.h. Wiener process we refer to [9.
Any arbitrary s.h. Wiener process W Γ   is uniquely determined by its covariance Γ   according to the formula
E W Γ ( t , θ ) , φ W Γ ( τ , θ ) , ψ = m i n ( t , τ ) Γ , φ ψ ( s ) , (2.1)
where φ , ψ D ( T d )   and ψ ( s ) ( η ) = ψ ( η )   , for η T d   . Because W Γ   is spatially homogeneous process, the distribution Γ = Γ ( θ η )   for θ , η T d   .
The space covariance Γ   , like distribution in D ( T d )   , may be uniquely expanded (see, e.g. [18or [5) into its Fourier series (with parameter w = 1   because the period is 2 π   )
Γ ( θ ) = n Z d e i ( n , θ ) γ n , θ T d , (2.2)
convergent in D ( T d )   . In ( 2.2 ), ( n , θ ) = i = 1 d n i θ i   and Z d   denotes the product of integers.
The coefficients γ n   , in the Fourier series ( 2.2 ), satisfy:
  • (1) γ n = γ n   for n Z d   ,
  • (2) the sequence ( γ n )   is slowly increasing, that is
    n Z d γ n 1 + | n | r < + , for some r > 0 . (2.3)
Let us introduce, by induction, the following set of indexes. Denote Z s 1 : = N   , the set of natural numbers and define Z s d + 1 : = ( Z s 1 × Z d ) { ( 0 , n ) : n Z s d }   . Let us notice that Z d = Z s d ( Z s d ) { 0 }   . For instance, for d = 2   , Z s 2 = N × Z { ( 0 , n ) : n Z }   .
Now, the s.h. Wiener process W Γ   corresponding to the covariance Γ   given by ( 2.2 ) may be represented in the form
W Γ ( t , θ ) = γ 0 β 0 ( t ) + n Z s d 2 γ n [ cos ( n , θ ) β n 1 ( t ) + sin ( n , θ ) β n 2 ( t ) ] , (2.4)
t 0 and θ T d .
In ( 2.4 ) β 0 , β n 1 , β n 2 , n Z s d   , are independent real Brownian motions and γ 0 , γ n   are coefficients of the series ( 2.2 ). The series ( 2.4 ) is convergent in the sense of D ( T d )   .
Because any periodic distribution with positive period is a tempered distribution (see, e.g [5), we may restrict our considerations to the space S ( T d )   of tempered distributions. By S ( T d )   we denote the space of infinitely differentiable rapidly decreasing functions on the torus T d   .
Let us denote by H α = H α ( T d )   , α R   , the real Sobolev space of order α   on the torus T d   . The norms in such spaces may be expressed in terms of the Fourier coefficients (see [1) | | ξ | | H α = ( n Z d ( 1 + | n | 2 ) α | ξ n | 2 ) 1 2 = ( | ξ 0 | 2 + 2 n Z s d ( 1 + | n | 2 ) α ( ( ξ n 1 ) 2 + ( ξ n 2 ) 2 ) ) 1 2 ,   where ξ n = ξ n 1 + i ξ n 2 , ξ = ξ ¯ n , n Z d   .
There is another possibility to define the Sobolev spaces (see, e.g. [17). We say that a distribution ξ S ( T d )   belongs to H α , α R   , if its Fourier transform ξ ^   is a measurable function and T d ( 1 + | λ | 2 ) α | ξ ^ ( | λ | ) | 2 d λ < + .  

3 Main results

If b L l o c 1 ( R + )   and μ C   , we shall denote by r ( t , μ )   the unique solution in L l o c 1 ( R + )   to the linear Volterra equation
r ( t , μ ) = b ( t ) + μ 0 t b ( t s ) r ( s , μ ) d s , t 0 . (3.1)
In many cases the function r ( t , μ )   may be found explicitly. For instance:
b ( t ) 1 , r ( t , μ ) = e μ t
b ( t ) = t , r ( t , μ ) = sinh μ t μ μ 0
b ( t ) = e t , r ( t , μ ) = e ( 1 + μ ) t
b ( t ) = t e t , r ( t , μ ) = e t sinh μ t μ .
Remark 3.1. Let us denote by f ^ ( k ) , k Z   , the k   -th Fourier coefficient of f   :
f ^ ( k ) = 1 2 π 0 2 π e i k t f ( t ) d t .   Given b L 1 ( R + )   , we find that, for F ( t ) : = t b ( t s ) f ( s ) d s   , we have
F ^ ( k ) = b ~ ( i k ) f ^ ( k ) , k Z , (3.2)
where b ~ ( λ ) : = 0 e λ t b ( t ) d t   denotes the Laplace transform of b   .
We state the following hypothesis:
Definition 3.2. We say that a kernel b L 1 ( R + )   is admissible for equation ( 1.1 ) if lim | n | | n | 2 0 [ r ( s , | n | 2 ) ] 2 d s = : C b   exists.
Example 3.3. i) In case b ( t ) = e t   we obtain | n | 2 0 [ r ( s , | n | 2 ) ] 2 d s = | n | 2 2 ( 1 | n | 2 )   and hence C b = 1 2 .   ii) In case b ( t ) = t e t   we obtain | n | 2 0 [ r ( s , | n | 2 ) ] 2 d s = | n | 2 4 + 4 | n | 2   and hence C b = 1 4 .  
Denote by W n ( t , θ ) : = cos ( n , θ ) β n 1 ( t ) + sin ( n , θ ) β n 2 ( t ) , n Z s d   , that is the n   -th element in the expansion ( 2.4 ).
Definition 3.4. By a solution X ( t , θ )   to the stochastic Volterra equation ( 1.1 ) we will understand the process of the form
X ( t , θ ) = γ 0 β 0 ( t ) + n Z s d 2 γ n t r ( t s , | n | 2 ) d W n ( s , θ ) , (3.3)
where the function r   is as above, t 0   and θ T d   .
The process X   given by ( 3.3 ) is a particular form of the weak solution to the equation ( 1.1 ) (cf. [9) and takes values in the space S ( T d )   .
The following is our main result.
Theorem 3.5. Assume b L 1 ( R + )   is admissible for ( 1.1 ). Then, the equation ( 1.1 ) has an H α + 1 ( T d )   -valued solution if and only if the Fourier coefficients ( γ n )   of the covariance Γ   satisfy
n Z d γ n ( 1 + | n | 2 ) α < . (3.4)
Remark 3.6. Particularly, the eq. ( 1.1 ) has L 2 ( T d )   -valued solution if and only if n Z d γ n 1 + | n | 2 < .  
Proof. We shall use the representation ( 2.4 ) for the Wiener process W Γ ( t , θ )   with respect to the basis: 1 , cos ( n , θ ) , sin ( n , θ )   , where n Z s d   and θ T d   . Equation ( 1.1 ) may be solved coordinatewise as follows.
Assume that
X ( t , θ ) = n Z s d [ cos ( n , θ ) X n 1 ( t ) + sin ( n , θ ) X n 2 ( t ) ] + X 0 ( t ) . (3.5)
Introducing ( 3.5 ) into ( 1.1 ), we obtain
cos ( n , θ ) X n 1 ( t ) + sin ( n , θ ) X n 2 ( t ) = | n | 2 t b ( t s ) [ cos ( n , θ ) X n 1 ( s ) + sin ( n , θ ) X n 2 ( s ) ] d s
+ 2 γ n t b ( t s ) [ cos ( n , θ ) β n 1 ( s ) + sin ( n , θ ) β n 2 ( s ) ] d s ,
or, equivalently
[ cos ( n , θ ) , sin ( n , θ ) ] [ X n 1 ( t ) X n 2 ( t ) ] = | n | 2 t b ( t s ) [ cos ( n , θ ) , sin ( n , θ ) ] [ X n 1 ( s ) X n 2 ( s ) ] d s
+ 2 γ n t b ( t s ) [ cos ( n , θ ) , sin ( n , θ ) ] [ d β n 1 ( s ) d β n 2 ( s ) ] .
Denoting X n ( t ) : = [ X n 1 ( t ) X n 2 ( t ) ] and β n ( t ) : = [ β n 1 ( t ) β n 2 ( t ) ]   we arrive at the equation
X n ( t ) = | n | 2 t b ( t s ) X n ( s ) d s + 2 γ n t b ( t s ) d β n ( s ) . (3.6)
Taking Fourier transform in t   , and making use of ( 3.2 ), we get the following solution to the equation ( 3.6 ):
X n ( t ) = t r ( t s , | n | 2 ) 2 γ n d β n ( s ) = 0 r ( s , | n | 2 ) 2 γ n d β n ( t s ) .   Hence, we deduce the following explicit formula for the solution to the equation ( 1.1 ):
X ( t , θ ) = γ 0 β 0 ( t ) + n Z s d 2 γ n [ cos ( n , θ ) 0 r ( s , | n | 2 ) d β n 1 ( t s ) + sin ( n , θ ) 0 r ( s , | n | 2 ) d β n 2 ( t s ) ] . (3.7)
Since the series defining the process X   converges in S ( T d ) , P   -almost surely, it follows from the definition of the space H α   that X ( t ) H α + 1 , P   -almost surely if and only if
n Z d ( 1 + | n | 2 ) α + 1 γ n [ ( t r ( t s , | n | 2 ) d β n 1 ( s ) ) 2 + ( t r ( t s , | n | 2 ) d β n 2 ( s ) ) 2 ] < . (3.8)
Because the stochastic integrals in ( 3.8 ) are independent Gaussian random variables, we obtain that ( 3.8 ) holds P   -almost surely if and only if
n Z d ( 1 + | n | 2 ) α + 1 γ n E [ ( t r ( t s , | n | 2 ) d β n 1 ( s ) ) 2 + ( t r ( t s , | n | 2 ) d β n 2 ( s ) ) 2 ] < . (3.9)
Or equivalently, using properties of stochastic integrals, if and only if
n Z d ( 1 + | n | 2 ) α + 1 γ n 0 [ r ( s , | n | 2 ) ] 2 d s < . (3.10)
Since b   is admissible for the equation ( 1.1 ), we conclude that ( 3.10 ) holds if and only if n Z d ( 1 + | n | 2 ) α + 1 γ n | n | 2 < ,   and the proof is complete.
Concerning uniqueness, we have the following result.
Proposition 3.7. Assume b L 1 ( R + )   is admissible for ( 1.1 ) and the following conditions hold:
(i) n Z d γ n ( 1 + | n | 2 ) α < ,   (ii) { 1 b ^ ( i k ) } k Z C \ { | n | 2 : n Z d }   .
Then, ( 1.1 ) has a unique H α + 1 ( T d )   -valued solution.
Proof. Let X ( t , θ )   be solution of X ( t , θ ) = t b ( t s ) Δ X ( s , θ ) d s .   Taking Fourier transform in θ   and denoting by X n ( t )   the n   -th Fourier coefficient of X ( t )   , we obtain X n ( t ) = | n | 2 t b ( t s ) X n ( s ) d s   for all n Z d .   Taking now Fourier transform in t   , we get ( 1 + | n | 2 b ^ ( i k ) ) X ~ n ( k ) = 0   for all k Z   and all n Z d .   According to (ii) we obtain that X ~ n ( k ) = 0   for all k Z   and all n Z d .   Hence, the assertion follows by uniqueness of the Fourier transform.
The following corollaries are an immediate consequence of Theorem  3.5 . The arguments are the same as in [7. We give here the proof for the sake of completeness.
Corollary 3.8. Suppose b L 1 ( R + )   is admissible for ( 1.1 ) and assume Γ L 2 ( T d ) .   Then the integro-differential stochastic equation ( 1.1 ) has a solution with values in L 2 ( T d )   for d = 1 , 2 , 3 .  
Proof. We have to check equation ( 3.4 ) with α = 1   . Note, that if Γ L 2 ( T d )   then Γ ^ = ( γ n ) l 2 ( Z d ) .   Consequently n Z d γ n 1 + | n | 2 ( n Z d γ n 2 ) 1 / 2 ( n Z d 1 ( 1 + | n | 2 ) 2 ) 1 / 2 .   But n Z d γ n 2 <   and
n Z d 1 ( 1 + | n | 2 ) q < if and only if 2 q > d . (3.11)
Hence, the result follows.
Corollary 3.9. Suppose b L 1 ( R + )   is admissible for ( 1.1 ) and assume Γ ^ l p ( Z d )   for 1 < p 2 .   Then the integro-differential stochastic equation ( 1.1 ) has a solution with values in L 2 ( T d )   for all d < 2 p p 1 .  
Proof. Note that n Z d γ n 1 + | n | 2 ( n Z d γ n p ) 1 / p ( n Z d 1 ( 1 + | n | 2 ) q ) 1 / q ,   where 1 p + 1 q = 1 .   Hence the result follows from ( 3.11 ) with q = p p 1 .   For α = 1   , the condition ( 3.4 ) can be written as follows.
Theorem 3.10. Let b   be admissible for the equation ( 1.1 ). Assume that the covariance Γ   is not only a positive definite distribution but is also a non-negative measure. Then the equation ( 1.1 ) has L 2 ( T d )   -valued solution if and only if
( Γ , G d ) < + , (3.12)
where
G d ( x ) = n Z d 0 + 1 ( 4 π t ) d e t e | x + 2 π n | 2 t d t , x T d . (3.13)
The proof of Theorem  3.10 is the same that for Theorem 2, part 2) in [7, so we omit it. For more details concerning the function G d   we refer to [7and [10.
Additionally, from properties of function G d   defined by ( 3.13 ) and the condition ( 3.12 ) we obtain the following result (see [8,Theorem4).
Corollary 3.11. Assume that Γ   is a non-negative measure. Then equation ( 1.1 ) has function valued solutions:
i) for all Γ   if d = 1   ; ii) for exactly those Γ   for which | θ | 1 ln | θ | Γ ( d θ ) < +   if d = 2   ; iii) for exactly those Γ   for which | θ | 1 1 | θ | d 2 Γ ( d θ ) < +   if d 3   .
Let us see that the formula ( 3.7 ) of the solution provides the next regularity result. In order to prove hölderianity of X   with respect to t   , we need assumptions very similar to those in [2.
Hypothesis (H) Assume that there exist δ ( 0 , 1 )   and C δ > 0   such that, for all τ ( , t )   we have:
(i) s t r 2 ( t τ , | n | 2 ) d τ C δ | n | 2 ( δ 1 ) | t s | δ   ; (ii) s [ r ( t τ , | n | 2 ) r ( s τ , | n | 2 ) ] 2 d τ C δ | n | 2 ( δ 1 ) | t s | δ   .
Proposition 3.12. Assume that n Z d γ n 1 + | n | 2 < +   . Under Hypothesis (H), the trajectories of the solution X   to the equation ( 1.1 ) are almost surely η   -Hölder continuous with respect to t   , for every η ( 0 , δ / 2 )   .
Proof. From the expansion ( 3.7 ) and properties of stochastic integral, we have
E | | X ( t , θ ) X ( s , θ ) | | L 2 2 = E | | γ 0 ( β 0 ( t ) β 0 ( s ) )
+ n Z s d 2 γ n [ cos ( n , θ ) ( t r ( t τ , | n | 2 ) d β n 1 ( τ ) s r ( s τ , | n | 2 ) d β n 1 ( τ ) )
+ sin ( n , θ ) ( t r ( t τ , | n | 2 ) d β n 2 ( τ ) s r ( s τ , | n | 2 ) d β n 2 ( τ ) ) ] | | L 2 2
= ( 2 π ) d ( γ 0 | t s |
+ n Z s d 2 γ n [ s [ r ( t τ , | n | 2 ) r ( s τ , | n | 2 ) ] 2 d τ + s t r 2 ( t τ , | n | 2 ) d τ ] ) .
According to assumptions (i) and (ii) of the Hypothesis (H), we have E | | X ( t , θ ) X ( s , θ ) | | L 2 2 C δ n Z s d 2 γ n | n | 2 ( δ 1 ) | t s | δ .   Because X   is a Gaussian process, then for any m N   , there exists a constant C m > 0   that E | | X ( t , θ ) X ( s , θ ) | | L 2 2 m C m [ C δ n Z s d 2 γ n | n | 2 ( δ 1 ) ] m | t s | m δ .   Taking m   such that m δ > 1   and using the Kolmogorov test, we see that the solution X ( t , θ )   is η   -Hölder continuous, with respect to t   , for η = δ / 2 1 / ( 2 m )   .
Example 3.13. Let us consider the particular case b ( t ) = e t , t 0   .
Then, by previous considerations, r ( t , | n | 2 ) = e ( 1 | n | 2 ) t   . One can check that in this case the Hypothesis (H) is fulfilled.
References

  1. R. Adams, Sobolev Spaces, Academic Press, New York, 1975.
  2. Ph. Clément, G. Da Prato, Some results on stochastic convolutions arising in Volterra equations perturbed by noise, Rend. Math. Acc. Lincei. s. 9, 7, (1996) 147–153.
  3. Ph. Clément, G. Da Prato, Existence and regularity results for an integral equation with infinite delay in a Banach space, Integral Equations and Operator Theory, 11, (1988), 480-500.
  4. R. Dalang, N. Frangos, The stochastic wave equation in two spatial dimensions, The Annals of Probability 26 (1998), 187-212.
  5. C. Gasquet, P. Witomski, Fourier Analysis and Applications, Springer-Verlag, New York Berlin, 1999.
  6. A. Karczewska, Regularity and continuity of solutions to stochastic evolution equations, G. Da Prato and L. Tubaro (eds.), Stochastic partial differential equations and applications, Marcel Dekker, Lect. Notes Pure Appl. Math. 227, 309-323 (2002).
  7. A. Karczewska, J. Zabczyk, A note on stochastic wave equations, G. Lumer and L. Weis (eds.), Evolution equations and their applications in physical and life sciences, Marcel Dekker, Lect. Notes Pure Appl. Math. 215, 501-511 (2001).
  8. A. Karczewska, J. Zabczyk, Stochastic PDE's with function-valued solutions. Ph. Clément (ed.) et al., Infinite dimensional stochastic analysis. Royal Netherlands Academy of Arts and Sciences. Verh. Afd. Natuurkd., 1. Reeks, K. Ned. Akad. Wet. 52, 197-216 (2000).
  9. A. Karczewska, J. Zabczyk, Regularity of solutions to stochastic Volterra equations, Rend. Math. Acc. Lincei. s. 9, 11, (2001) 141–154.
  10. N.S. Landkof, Foundations of modern potential theory, Springer-Verlag, Berlin, 1975.
  11. O. Lévêque, Hyperbolic stochastic partial differential equations driven by boundary noises, Ph.D. Thesis, Ecole Polytechnique Fédérale de Lausanne, 2001.
  12. R.C. Mac Camy, Stability theorems for a class of functional differential equations, SIAM J. Appl. Math., 30 (1976), 557-576.
  13. C. Mueller, Long time existence for the wave equation with a noise term, The Annals of Probability 25 (1) (1997), 133-151.
  14. J.W. Nunziato, On heat conduction in materials with memory, Quart. Appl. Math. 29, (1971), 187-204.
  15. G. Da Prato, A. Lunardi, Periodic solutions for linear integrodiferential equations with infinite delay in Banach spaces. Differential Equations in Banach spaces, Springer-Verlag, Lecture Notes in Mathematics 1223, 49-60 (1985).
  16. J. Prüss. Evolutionary Integral Equations and Applications. Monographs Math., vol. 87, Birkhäuser Verlag, Basel, 1993.
  17. M. Reed, B. Simon, Methods of modern mathematical physics, Vol. II, Academic Press, New York, 1975.
  18. L. Schwartz, Méthodes mathématiques pour les sciences physiques, Hermann, Paris,1965.

Department of Mathematics, University of Zielona Gora, ul. Szafrana 4a, 65-246 Zielona Gora, Poland. E-mail address : A.Karczewska@im.uz.zgora.pl Universidad de Santiago de Chile, Departamento de Matematica, Facultad de Ciencias, Casilla 307-Correo 2, Santiago, Chile. E-mail address : clizama@lauca.usach.cl