The heat equation with multiplicative stable Lévy noise

Carl Mueller 1 Dept. of Mathematics University of Rochester Rochester, NY 14627 E-mail: cmlr@math.rochester.edu

Leonid Mytnik 2 Faculty of Industrial Engineering and Management Technion – Israel Institute of Technology Haifa 32000, Israel E-mail: leonid@ie.technion.ac.il

Aurel Stan 1 Dept. of Mathematics University of Rochester Rochester, NY 14627 E-mail: astan@math.rochester.edu

Abstract
We study the heat equation with a random potential term. The potential is a one-sided stable noise, with positive jumps, which does not depend on time. To avoid singularities, we define the equation in terms of a construction similar to the Skorokhod integral or Wick product. We give a criterion for existence based on the dimension of the space variable, and the parameter p   of the stable noise. Our arguments are different for p < 1   and p ≥ 1   .

1 Introduction

Our goal is to study the heat equation with stable Lévy noise L ˙ ( x )   which depends only on space.
Assumption 1 Throughout the paper, we assume that L ˙   has only positive jumps.
Our noise is multiplicative, in the sense that it is multiplied by u   in the equation.
∂ u ∂ t = Δ u + u ⋄ L ˙ ( x ) (1.1)
u ( 0 , x ) = u 0 ( x ) .
Here, L ˙ ( x )   is a one-sided stable noise of index p ≥ 1   defined on x ∈ R d   .
The product u ⋄ L ˙ ( x )   is related to the Skorokhod integral or Wick product.
We will give precise conditions later. Our main goal is to give conditions on p   and the dimension d   such that the solution exists and is unique. In fact, for p < 1   , the stochastic integral is defined in a different way and can include more singular terms. In this case, we use the notation
∂ u ∂ t = Δ u + u ◆ L ˙ ( x ) (1.2)
u ( 0 , x ) = u 0 ( x ) .
There are many papers concerning the heat equation with a multiplicative noise term (see [Par93] ). In most cases, however, the noise depends on time, and is even independent from one time point to another. This is a very natural assumption from the point of view of Itô integration, since the relevant integrals often yield martingales. For example, the equation
∂ u ∂ t = Δ u + f ( u ) W ˙ ( t , x ) (1.3)
u ( 0 , x ) = u 0 ( x )
where W ˙ ( t , x )   is space-time white noise, is usually formulated in terms of an integral equation
u ( t , x ) = ∫ R d G ( t , x − y ) u 0 ( y ) d y
+ ∫ 0 t ∫ R d G ( t − s , x − y ) f ( u ( s , y ) ) W ( d y d s )
where G ( t , x )   is the heat kernel, and where the final integral can be defined in much the same way as the Itô integral.
But if the noise does not depend on time, such an integral would be anticipating. Further evidence of the difficulty of noise depending only on space is given by the following equation:
∂ u ∂ t = Δ u + δ ( x ) u (1.4)
u ( 0 , x ) = u 0 ( x )
where δ ( x )   is the Dirac delta function. When u   does not depend on t   , this equation has been extensively studied by mathematical physicists. See the book [AGHKH88] by Albeverio et al., for example. This equation (without t   dependence) models the quantum mechanics of a particle under the influence of a point interaction. As usual, we can think of this equation in terms of the density of Brownian particles. The δ ( x )   term means that Brownian particles at 0 would give birth to many new particles. But these new particles would themselves be close to 0, and so they might give birth to other particles, ad infinitum. Because of this unstable behavior, when the dimension d   is 2 or higher we need to take an approximate δ ( x )   function, multiply it by a very small term, and take the limit of the approximation. The case of x ∈ R   is better behaved.
Yaozhong Hu has written two papers on time-independent noise, [Hu01] and [Hu02] . He deals with the equation
∂ u ∂ t = Δ u + u ⋄ W ˙ ( x )
u ( 0 , x ) = u 0 ( x )
where ⋄   denotes Skorokhod integration. One can also think of this product as the Wick product, see [Jan97] Chapter 3. W ˙ ( x )   is white noise in [Hu02] , but in [Hu01] it is colored Gaussian noise. By the Skorokhod integral, roughly speaking, we mean that in the relevant stochastic integrals, we should drop diagonal terms such as ( W ˙ ( x ) ) 2   . These terms correspond to repeated influence of W ˙ ( x )   for the same point x   , so they are related to the singular behavior in the equation ( 1.4 ).
The goal of this article is to study ( 1.1 ) for time-independent Lévy noise L ˙ ( x )   . The product ⋄   , like the Wick product, involves the deletion of diagonal terms in the appropriate integrals. But it is not defined in the same way as the Wick product for Brownian functionals, at least with the usual definitions.
Once again, the reason for using such a product is to avoid singular terms such as those coming from ( 1.4 ). The product ◆   used in ( 1.2 ) is similar, but we delete fewer terms.
Our idea is that, once we understand the equation without the singular terms, we can try gingerly adding back the singular terms. But the understanding of the simpler situation must come first.
In the time dependent situation, there have been several papers involving Lévy noise, such as Kallianpur, Xiong, Hardy, and Ramasubramanian [KXHR94] .
In most cases, the noise has been additive. That is, it appears without any multiple of u   or f ( u )   .
In the case p < 1   , and for noise depending on both time and space, [Mue98] dealt with the equation
∂ u ∂ t = Δ u + u γ L ˙ ( t , x ) (1.5)
u ( 0 , x ) = u 0 ( x ) .
If x ∈ D   for a smooth and bounded domain D ⊂ R d   , and with some conditions on u 0   , short-time existence was obtained in the case d < ( 1 − p ) α γ p − ( 1 − p ) .   For p ≥ 1   , Mytnik [Myt02] found that if
p γ < 2 d + 1 ,
1 < p < min ( 2 , 2 d + 1 )
then there exists a weak solution to ( 1.5 ).
Both of these papers were motivated by the superprocess with stable branching, see [Daw93] . We will not describe this superprocess, but it is believed that in some weak sense, it is associated with the equation ( 1.5 ) with γ = 1 / p   . This motivation led us to study equations with one-sided noise. We leave possible generalizations to the reader.
Our results are also divided into the two cases p < 1   and p ≥ 1   . Our main tools are multiple stochastic integrals, and this is the reason that we restrict ourselves to noise terms of the form f ( u ) ⋄ L ˙ ( x )   or f ( u ) ◆ L ˙ ( x )   with f ( u ) = u   . Hu restricted himself to similar equations in the Gaussian case, for the same reason. Of course, multiple stochastic integrals have been used many times to study equations such as ( 1.3 ). For example, see Nualart and Zakai [NZ89] and Nualart and Rozovskii [NR97] . But in these articles, the noise always depends on time.

2 Theorems

In this section we list our main results. They depend on some definitions which will be explained later. We assume that u 0 ( x )   is a bounded function on D   . Since our equation is linear, we may assume without loss of generality that ∥ u 0 ∥ ∞ ≤ 1   .
Here are our main theorems.
Theorem 1 Suppose that 0 < p < 1   , and let u 0 ( ⋅ )   be a bounded function on D   . Assume that either d ≤ 4   or p < 1 2 + 1 d − 2 , d ≥ 5 .   Then ( 1.2 ) has a unique solution u ( t , x )   .
We are not sure that the conditions on p , d   in the above theorem are optimal.
Theorem 2 Suppose that p ≥ 1   and u 0 ( ⋅ )   is a bounded function on D   . Assume that p < 1 + 2 d .   Then equation ( 1.1 ) has a unique solution u ( t , x )   . Moreover, for all t > 0   and x ∈ D   , u ( t , x ) ∈ L l o c q ( Ω )   , for all q ∈ [ 1 , 1 + ( 2 / d ) )   .
Here L l o c q ( Ω )   denotes the space of all random variables f : Ω → R   , such that, for all ε > 0   , there exists an event A ⊂ Ω   , of probability less than ε   , such that E [ | f 1 A c | q ] < ∞   .
Remarks:
  • 1. In Theorem  2 , in the Brownian case of p = 2   , the critical dimension is d = 2   , which agrees with [Hu02] .
  • 2. For p < 1   , if we use the product ⋄   instead of ◆   , the proof of Theorem  1 is trivial. The reader can check that by removing the large atoms of the measure L ˙ ( x )   , we can take q = 1   and use the proof of Theorem  2 .

3 Basic Definitions

We will index sequences of real numbers, such as y i   , by subscripts. A sequence of vectors in R n   will be indexed by superscripts. Thus, x ( i )   is the i t h   vector in the sequence, and x k ( i )   is the k t h   component of the i t h   vector.
If it is not stated otherwise, C   in this paper will denote a constant which value my change from line to line.
Throughout the paper, G ( t , x )   denotes the fundamental solution of the heat equation for x ∈ R d   ,
∂ u ∂ t = Δ u
u ( 0 , x ) = δ ( x ) .
Explicitly, G ( t , x ) = 1 ( 4 π t ) d / 2 exp ( − x 2 4 t ) .   Given a domain D ⊂ R d   , we let G ( t , x , y )   be the fundamental solution of the heat equation in D   , with Dirichlet boundary conditions. That is, G ( t , x , y )   satisfies
∂ u ∂ t = Δ u
u ( 0 , x ) = δ y ( x )
u ( t , x ) = 0 for x ∈ ∂ D .
We note the elementary inequality 0 ≤ G ( t , x , y ) ≤ G ( t , x − y ) .   In other words, when G   has two arguments, it is the heat kernel on R d   , and when it has three arguments, it is the heat kernel on D   .
Finally, we write T n = T n ( t ) : = { ( s 1 , . . . , s n ) : 0 < s 1 < . . . < s n < t } .  

4 The Noise

Roughly speaking, our noise L ˙ ( x )   will be the analogue of d L ( t ) / d t   where L ( t )   is a one-sided stable process of index p ∈ ( 0 , 2 )   , and t ∈ [ 0 , ∞ )   . Of course, this derivative only exists in the generalized sense. For example, we could use the theory of Schwartz distributions.
We will start with a review of the one-sided stable Lévy processes, and refer the reader to Bertoin's book [Ber96] for details. There is a basic difference between the cases p < 1   and p ≥ 1   .
For p < 1   , we can construct L ( t )   via a Poisson process, as follows.
Consider the following measure on ( t , y ) ∈ [ 0 , ∞ ) 2   given by ν ( d y d t ) = c p y − ( p + 1 ) d y d t ,   where c p   is a constant depending on p   . This is the same constant that occurs in the Lévy measure of the one-sided stable processes (see [Ber96] ). Let Y ( ⋅ )   be a Poisson random measure on [ 0 , ∞ ) 2   with intensity ν   . Then, the p   -stable process L ( t )   is defined as L t = ∫ [ 0 , t ] ∫ 0 ∞ y Y ( d y d s ) .   Note that if we let ( t i , y i )   be the locations of the points (or atoms) of the Poisson measure Y   , then L t = ∑ t i ≤ t y i .   In a similar way, we define the random measure L ( A ) : A ⊂ [ 0 , ∞ )   as L ( A ) = ∑ t i ∈ A y i   and observe that L t = ∫ 0 t L ( d s ) .   For p ≥ 1   , we must introduce an approximation and compensation procedure, as follows. Let Y n ( d y d t ) = Y ( d y d t ) 1 ( y ≥ 1 / n ) .   Note that if we define ν n ( d y d t ) : = ν ( d y d t ) 1 ( y ≥ 1 / n )   then Y n   becomes a Poisson measure with intensity ν n   . Note that ν n ( [ 0 , t ] × [ 0 , ∞ ) ) = ν ( [ 0 , t ] × [ 1 / n , ∞ ) ) = c p p n p t .   Next, let
L t ( n ) = ∫ 0 t ∫ 0 ∞ y Y n ( d y d s ) − t ∫ 0 ∞ y d ν n ( y )
= ∫ 0 t ∫ 0 ∞ y Y n ( d y d s ) − c p p − 1 n p − 1 t .
There is a problem with p = 1   . We leave the details to the reader, but roughly speaking, we would truncate the measure by removing large jumps.
Note that L t ( n )   is a martingale with respect to its own filtration. The p   -stable process L ( t )   is the limit of the processes L t ( n )   as n → ∞   . For details of this limit, we refer the reader to [Ber96] . We could define the random measure L ( d t )   in a similar way.
The notation L ˙ ( t )   stands for the the density of the measure L   , even though this is only defined in the generalized sense. This is similar to the convention of writing B ˙ ( t )   for the derivative of Brownian motion.
Our construction of L ( A ) : A ⊂ R d   is similar. First consider the case p < 1   . Fix a bounded open region D ⊂ R d   ; we will give some conditions on D   later. For ( x , y ) ∈ B ⊂ D × [ 0 , ∞ )   , let Y ( B )   be the Poisson measure with intensity d x ν ( d y )   . Let ( x ( i ) , y i )   denote the locations of the points (atoms) of the Poisson measure Y   . For A ⊂ D   , we define L ( A ) = ∑ x ( i ) ∈ A y i .   For p ≥ 1   , we must approximate. Given n   , let Y , ν   be defined as above, and let
Y n ( d x d y ) = 1 ( y ≥ 1 / n ) Y ( d x d y )
d x ν ( d y ) = c p y − ( p + 1 ) d x d y
d x ν n ( d y ) = 1 ( y ≥ 1 / n ) d x ν ( d y ) .
As above, for A ⊂ D   , we let
L n ( A ) = ∑ x ( i ) ∈ A y i 1 ( y i ≥ 1 / n ) − ∫ 0 ∞ ∫ A y d x ν n ( d y )
= ∑ x ( i ) ∈ A y i 1 ( y i ≥ 1 / n ) − m ( A ) c p p − 1 n p − 1
where m ( A )   denotes the Lebesgue measure of A   . Again, we leave the case p = 1   to the reader. To get L ( d x )   , we again take a limit as n → ∞   . Since this procedure is not very different than for the one-sided stable processes with p ∈ [ 1 , 2 )   , we refer the reader to [Ber96] .
For the future, we label the atoms of the measure L   by ( x ( i ) , y i ) : x ( i ) ∈ R d   .
Then, x ( i )   is the location of the atom, and y i   is the mass of the measure L ( { x ( i ) } )   .
Fix K > 0   . Since D   is a bounded domain, there are only a finite number of atoms ( x ( i ) , y i )   with x ( i ) ∈ D   and y i > K   . Let L K   denote L   with the preceding atoms deleted. That is, L K = L − ∑ y i > K y i δ x ( i ) .  
Assumption 2 In the succeeding sections, we will tacitly assume that L ˙   is actually L ˙ K   . Furthermore, let A K   denote the event that L = L K   , that is, there are no atoms larger than K   .
Note that lim K → ∞ P ( A K ) = lim K → ∞ P ( L = L K ) = 1   From now on, we will replace L   by L K   , and let K → ∞   at the end.

5 Multiple Stochastic Integrals

In the one-dimensional case, such integrals have been considered earlier, for example in [RW86] . We suppose that D ⊂ R d   is an open, bounded domain with smooth boundary. As described earlier, for p ∈ ( 0 , 2 )   , let L ( A )   be the one-sided stable random measure of index p   , defined for Borel sets A ⊂ D   .
Assume that p ≥ 1   . Now we describe our multiple integrals with respect to the noise L ˙   , recalling that L ˙ ( x ) d x   is just another notation for L ( d x )   .
Assume that we have a symmetric function f n ( x ( 1 ) , . . . , x ( n ) )   defined for x ( 1 ) , . . . , x ( n ) ∈ D   . Define D ~ n ⊂ D n   as follows.
D ~ n = { ( x ( 1 ) , x ( 2 ) , . . . , x ( n ) ) ∈ D n : x 1 ( 1 ) < x 1 ( 2 ) < . . . < x 1 ( n ) } .   We regard x 1   as the time variable, so that the following integral is defined in the Itô sense.
I n ( f n ) = n ! ∫ D ~ n f n ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) .   Sufficient conditions for the existence of the integral I n ( f n )   are given in Lemma  3 .
Since f n   is symmetric, I n ( f n )   is invariant under permutations of the indices i = 1 , . . . , n   of the x ( i )   . For f = ( f 0 , f 1 , . . . )   , where f n   is a symmetric function of the variables x ( 1 ) , . . . , x ( n )   , let I ( f ) = ∑ n = 0 ∞ I n ( f n ) .   Here, f 0   is interpreted as a constant, and I 0 ( f 0 ) = f 0   .
For future use, we define the symmetrization of a function as follows.
s y m ( f ) ( x ( 1 ) , . . . , x ( n ) ) = 1 n ! ∑ σ f ( x ( σ ( 1 ) ) , . . . , x ( σ ( n ) ) )   where the sum is taken over all permutations σ   of the indices 1 , . . . , n   . If we wish to symmetrize over only some of the variables, we write those in a subscript.
s y m x ( 1 ) , . . . , x ( n ) ( f ) ( x ( 1 ) , . . . , x ( n ) , z ( 1 ) , . . . , z ( m ) )
= 1 n ! ∑ σ f ( x ( σ ( 1 ) ) , . . . , x ( σ ( n ) ) , z ( 1 ) , . . . , z ( m ) )
Note that the integral I ( f )   has the following property. Because of the ordering we have imposed on the variables x ( i )   , it follows that I ( f )   depends linearly on the masses y i   of the atoms of the measure L   . That is, I ( f )   can have terms like y 1 y 2   , but no quadratic terms like y 1 2   .
We note in passing that our multiple stochastic integrals are related to the Wick product for the compensated Poisson process Y   ; see Lytvynov [Lyt03] , for example.
For the case p < 1   , we propose a different integral J ( f )   , which allows for more singularities. Let X n ( D ) ⊂ D n   be the set of coordinates ( x ( 1 ) , . . . , x ( n ) )   with no adjacent indices equal. In other words, x ( i ) ≠ x ( i + 1 )   for i = 1 , . . . , n − 1   . Let J n ( f n ) = ∫ X n ( D ) f n ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) .   Since L ( d x ( 1 ) ) . . . L ( d x ( n ) )   is a measure, the integral J n ( f n )   can be defined for each ω   , provided f n   is integrable with respect to this product measure.
We define J ( f )   in a similar way:
J ( f ) = ∑ n = 0 ∞ J n ( f n ) .  

6 A Framework for the Equation

The goal of this section is to set up a framework for discussing the equation, so that existence and uniqueness can be rigorously discussed. In the white noise case, there are some related concepts in [NR97] and [MR98] .
We define
ℐ n ≡ { I n ( f n ) : f n is a symmetric function on D n } . (6.1)
Definition 1 Let g ( ⋅ , ⋅ )   be a measurable function on D × Ω   , such that g ( x , ⋅ ) ∈ ℐ n   for almost every x ∈ D   . That is, there exists a function f n   on D n + 1   such that g ( x , ⋅ ) = n ! ∫ D ~ n f n ( x ( 1 ) , . . . , x ( n ) , x ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) , a . e . x ,   and f n ( ⋅ , x )   is a symmetric function on D n   for almost every x   .
Then we define
∫ D g ( x ) ⋄ L ( d x )
≡ ( n + 1 ) ! ∫ D ~ n + 1 s y m x ( 1 ) , . . . , x ( n + 1 ) f n ( x ( 1 ) , . . . , x ( n ) , x ( n + 1 ) )
L ( d x ( 1 ) ) . . . L ( d x ( n + 1 ) ) .
Set ℐ = { ∑ n = 0 ∞ g n : g n ∈ ℐ n a n d t h e s u m i s c o n v e r g i n g i n p r o b a b i l i t y } .  
Definition 2 Let g ( ⋅ , ⋅ )   be a measurable function on D × Ω   , such that g ( x , ⋅ ) ∈ ℐ   for almost every x ∈ D   . That is, there exists a sequence { g n } n ≥ 0   such that g n ( x , ⋅ ) ∈ ℐ n , n ≥ 0 , a . e . x ∈ D ,   and g ( x , ⋅ ) = ∑ n = 0 ∞ g n ( x , ⋅ ) , a . e . x ∈ D .   Then we define ∫ D g ( x ) ⋄ L ( d x ) ≡ ∑ n = 0 ∞ ∫ D g n ( x ) ⋄ L ( d x ) ,   where the sum is converging in probability.
Next definition is analogous to the definition of Hu, from [Hu01] .
Definition 3 A measurable function u : R + × D × Ω ↦ R   is called a solution to the equation (1.1) if u ( t , x ) ∈ ℐ   for every x ∈ R d   and t ∈ ( 0 , T ]   and the following equation is satisfied:
u ( t , x ) = ∫ D G ( t , x , y ) u 0 ( y ) d y (6.2)
+ ∫ D ( ∫ 0 t G ( t − s , x , y ) u ( s , y ) d s ) ⋄ L ( d y ) ,
∀ x ∈ R d , t ∈ ( 0 , T ] .
Our strategy is to expand u ( t , x )   in a recursively defined series:
u ( t , x ) = ∑ n = 0 ∞ u n ( t , x ) ,   where
u 0 ( t , x ) = ∫ D G ( t , x , y ) u 0 ( y ) d y (6.3)
u n + 1 ( t , x ) = ∫ 0 t ∫ D G ( t − s , x , y ) u n ( s , y ) ⋄ L ( d y ) d s .
Note that u 0   has two meanings. If u 0 ( x )   is a function of one variable, then it is the initial value for our SPDE ( 1.1 ). If u 0 ( t , x )   depends on two variables, then it is the first term of the expansion for u ( t , x )   .
To deal with the case p < 1   , we define the operation ◆   . We define
J n ≡ { J n ( f n ) : f n is a function on D n } . (6.4)
Definition 4 Let g ( ⋅ , ⋅ )   be a measurable function on D × Ω   , such that g ( x , ⋅ ) ∈ J n   for almost every x ∈ D   . That is, there exists a function f n   on D n + 1   such that g ( x , ⋅ ) = ∫ X n ( D ) f n ( x ( 1 ) , . . . , x ( n ) , x ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) , a . e . x ,   and f n   is a function on D n + 1   .
Then we define
∫ D g ( x ) ◆ L ( d x )
≡ ∫ X n + 1 ( D ) f n ( x ( 1 ) , . . . , x ( n + 1 ) ) L ( d x ( 1 ) ) . . . L ( d x ( n + 1 ) ) .
The other definitions for p < 1   are completely analogous to the case p ≥ 1   , with ◆   instead of ⋄   . We leave this part to the reader.

7 Proof of the existence part of Theorem  1 ( p < 1   )

7.1 Calculations for the heat equation

As mentioned at the end of the Section  4 , we will assume that A K   occurs, so that L ˙ = L ˙ K   . We suppose that K = 2 N   for some integer N   . We will obtain conclusions that hold almost surely for each integer K = 2 N   . Since P ( A K ) → 1   as K → ∞   , our assertion will almost surely hold for L   .
We will use the notation of Mueller [Mue98] . We call the atom at x ( i )   with mass y i   a particle of type n   if 2 − ( n + 1 ) < y i ≤ 2 − n .   Since there are no atoms of mass greater than K = 2 N   , we need only consider particles of type n ≥ − N   .
We will use the expansion defined in ( 6.3 ), with ◆   instead of ⋄   . Then u 0 ( t , x )   is the solution of the heat equation with initial function u 0 ( x )   .

7.2 Distance between points

We start with an elementary lemma about Poisson probabilities.
Lemma 1 Let X   be a Poisson random variable with parameter λ   . Then for any nonnegative integer n 0   we have P ( X ≥ n 0 ) ≤ λ n 0 n 0 ! .   As a consequence P ( X ≥ n 0 ) ≤ λ n 0 n 0 n 0 e n 0 .  
Proof.
P ( X ≥ n 0 ) = ∑ n = n 0 ∞ P ( X = n )
= ∑ n = n 0 ∞ λ n n ! e − λ
= λ n 0 n 0 ! ∑ n = n 0 ∞ λ n − n 0 ( n − n 0 ) ! ⋅ n 0 ! ( n − n 0 ) ! n ! e − λ
But for all n ≥ n 0   , we have n 0 ! ( n − n 0 ) ! n ! = 1 ( n n 0 ) ≤ 1 .   Thus,
P ( X ≥ n 0 ) ≤ λ n 0 n 0 ! ∑ n = n 0 ∞ λ n − n 0 ( n − n 0 ) ! e − λ
= λ n 0 n 0 ! ∑ k = 0 ∞ λ k k ! e − λ
= λ n 0 n 0 ! .
Since
e n 0 = 1 + n 0 1 ! + n 0 2 2 ! + . . . + n 0 n 0 n 0 ! + . . .
≥ n 0 n 0 n 0 !
we have P ( X ≥ n 0 ) ≤ λ n 0 n 0 ! ≤ λ n 0 e n 0 n 0 n 0 .   Suppose that D ⊂ B   , where B   is a cube in R d   with side r ≥ 1   . Fix a , ℓ > 0   . Consider the set of closed cubes of side length ℓ   whose centers lie in the scaled lattice a Z d   . Let C ℓ ( B , a )   be the subset of such cubes which intersect B   .
Suppose that x , y ∈ D   is a pair of points whose distance | x − y | ≤ a   .
We claim that there exists a cube in C 3 a ( B , a )   which contains both x   and y   . Indeed, since the cubes C a ( B , a )   cover B   , it follows that x ∈ S a   for some cube S a ∈ C a ( B , a )   . Let S 3 a   be the cube with the same center as S a   , but with side length 3 a   . It is easy to see that every point within distance a   of S a   lies in S 3 a   . Therefore, S 3 a   is a cube from the set C 3 a ( B , a )   which contains both x   and y   .
Furthermore, let R ⊂ R d   and let N ( R , n )   be the number of particles of type n   which lie in R   . Let | R |   be the Lebesgue measure of R   . Note that N ( R , n )   is a Poisson random variable with parameter
λ ( R , n ) = | R | ∫ 2 − ( n + 1 ) 2 − n ν ( d y ) (7.1)
= | R | ∫ 2 − ( n + 1 ) 2 − n c p y p + 1 d y
= C p | R | 2 n p .
where C p = c p ( 2 p − 1 ) / p   .
For a , M > 0   , let A ( n , m , M , S )   be the event that within S ∩ B   , there is at least one particle of type n   and at least M   particles of type m   , where S   is a cube in C 3 a ( B , a )   . Let S ~ = S ∩ B   and note that | S ~ | ≤ ( a ∧ r ) d   where r   is the side of B   . We will estimate P ( A ( n , m , M , S ) )   for M = 2 k   , k ≥ 0   . There are two cases.
Case 1: n ≠ m   .
According to Lemma  1 ,
P ( A ( n , m , M , S ) ) = P ( N ( S ~ , n ) ≥ 1 ) P ( N ( S ~ , m ) ≥ M )
≤ C e λ ( S ~ , n ) λ ( S ~ , m ) M e M M M
= C e | S ~ | M + 1 2 n p + m M p e M M M
≤ C e ⋅ ( a ∧ r ) d ( M + 1 ) 2 n p + m M p e M M M .
Case 2: n = m   .
P ( A ( n , m , M , S ) ) = P ( N ( S ~ , m ) ≥ M + 1 )
≤ C λ ( S ~ , n ) M + 1 e M M M
= C | S ~ | M + 1 2 n p + m M p e M M M
≤ C ( a ∧ r ) d ( M + 1 ) 2 n p + m M p e M M M .
Next, for some p 0 > 0   , let a n , m , M = e − 1 / d δ 4 / ( M d ) 2 − p n / ( d M ) 2 − p m / d M 1 / d 2 − δ ( n + m + k ) / ( d M ) .   Note that a n , m , M d M = δ 4 [ 2 n p + m M p e M M M ] − 1 2 − δ ( n + m + k )   and therefore, for S ∈ C 3 a n , m , M ( B , a )   ,
( a n , m , 2 k ∧ r ) − d P ( A ( n , m , M , S ) )
≤ ( a n , m , 2 k ∧ r ) − d ( a n , m , 2 k ∧ r ) d ( M + 1 ) 2 n p + m M p e M M M
= ( a n , m , 2 k ∧ r ) d M 2 n p + m M p e M M M
≤ δ 4 2 − δ ( n + m + k ) .
Let A ( B )   be the event that there exist n , m ≥ − N   and M = 2 k   with k ≥ 0   , such that there exists a particle of type n   in B   , and within distance a n , m , M   of this particle there are M   particles of type m   in B   . Note that the number of cubes in C 3 a ( B , a )   is the same as the number of cubes in C a ( B , a )   . This number is less than or equal to a constant times ( a ∧ r ) − d   . Therefore, we have the following. We write ∑ S   for the sum over S ∈ C 3 a n , m , 2 k ( B , a n , m , 2 k )   .
As before, S ~ = S ∩ B   .
P ( A ( B ) )
≤ ∑ k = 0 ∞ ∑ n , m = − N ∞ ∑ S P ( A ( n , m , M , S ) )
= ∑ k = 0 ∞ ∑ n , m = − N ∞ ( a n , m , 2 k ∧ r ) d ∑ S ( a n , m , 2 k ∧ r ) − d P ( A ( n , m , M , S ) )
≤ C ∑ k = 0 ∞ ∑ n , m = − N ∞ ( a n , m , 2 k ∧ r ) d ∑ S δ 4 2 − δ ( n + m + k )
= C ∑ k = 0 ∞ ∑ n , m = − N ∞ ( a n , m , 2 k ∧ r ) d δ 4 2 − δ ( n + m + k ) # S
= C ∑ k = 0 ∞ ∑ n , m = − N ∞ δ 4 2 − δ ( n + m + k )
= C δ 4 ∑ k = 0 ∞ 2 − δ k ( ∑ n = − N ∞ 2 − δ n ) 2
= C δ 4 1 1 − 2 − δ ( 2 δ N 1 − 2 − δ ) 2
≤ C δ
where C   depends on N   , but not on δ   , as long as δ < 1   . Here, # S   denotes the number of cubes S   that we are considering.
Thus, we have shown:
Lemma 2 There exists a constant C > 0   such that the following holds. Fix a ball B ⊂ R d   of radius at least 1, let δ ∈ ( 0 , 1 )   , and let a n , m , M = e − 1 / d δ 4 / ( M d ) 2 − p n / ( d M ) 2 − p m / d M 1 / d 2 − δ ( n + m + k ) / ( d M ) .   Let A   be the event that for some n , m , M   , there exists a particle of type n   in B   with M   particles of type m   in B   within distance a n , m , M   . Then P ( A ) ≤ C δ .  

7.3 Estimation of the solution for the heat equation.

Now we will use Lemma  2 to bound the solution u ( t , x )   . Assume that event A   does not occur, where this event was defined in Lemma  2 .
Taking x   as fixed, define A 1   to be the event that for some m , M   , there exist M   particles of type m   within distance a 0 , m , M   of x   . The same arguments as before easily lead to the estimate P ( A 1 ) ≤ C δ   where C   is the same constant as in Lemma  2 .
For the rest of the section, fix n   , and consider u n ( t , x )   . Here n   is no longer a subscript for a n , m , M   . We will assume that neither A 1   nor A   occur. Also, we write
x ( n + 1 ) = x
M = 2 k .
Define ℳ ( n ) = { ( i 1 , . . . , i n ) ∈ N n : ∀ 1 ≤ k ≤ n − 1 , i k ≠ i k + 1 } .   Let v ( t , ⋅ )   be the solution of the heat equation on D   with initial function u 0   , and recall that G ( t , x , y ) ≤ G ( t , x − y )   . Setting i n + 1 = n + 1   and m n + 1 = 0   , we get
| u n ( t , x ) | 1 ( A c ∩ A 1 c ) (7.2)
≤ ∑ ( i 1 , . . . , i n ) ∈ ℳ ( n ) ∫ T n
n ∏ j = 1 [ G ( s j + 1 − s j , x ( i j + 1 ) − x ( i j ) ) y i j ] | v ( s 1 , x ( i 1 ) ) | d s
≤ ∥ u 0 ∥ ∞ ∑ k 1 , . . . , k n = 0 ∞ ∑ m 1 , . . . , m n = − N ∞ ∫ T n
n ∏ j = 1 [ G ( s j + 1 − s j , a m j + 1 , m j , 2 k j ) 2 k j + 1 2 − m j ] d s .
We assumed that ∥ u 0 ∥ ∞ ≤ 1   . It is known that the solution of the heat equation is bounded when the initial function u 0 ( x )   is bounded, and ∥ v ( t , ⋅ ) ∥ ∞ ≤ ∥ u 0 ∥ ∞   . When summing over all ( i 1 , . . . , i n ) ∈ ℳ n   we sum first over all possible types m 1 , . . . , m n   of particles x ( i 1 ) , . . . , x ( i n )   . For two fixed consecutive particles x ( i j )   and x ( i j + 1 )   , we let k j   be the smallest nonnegative integer k   such that a m j + 1 , m j , 2 k ≤ | x ( i j + 1 ) − x ( i j ) | < a m j + 1 , m j , 2 k + 1   where m j + 1   and m j   are the types of the particles x ( i j + 1 )   and x ( i j )   . Observe that such a k j   always exists. Because we are on the complement of both sets A   and A 1   , there is no particle of type m j   at a distance less than a m j + 1 , m j , 0   from x ( i j + 1 )   . Thus, | x ( i j + 1 ) − x ( i j ) | ≥ a m j + 1 , m j , 2 0 .   On the other hand, for fixed m j + 1   and m j   , we have lim k → ∞ a m j + 1 , m j , 2 k = ∞ .   Thus for large enough k   , we have | x ( i j + 1 ) − x ( i j ) | ≥ a m j + 1 , m j , 2 k .   Therefore, G ( s j + 1 − s j , x ( i j + 1 ) − x ( i j ) ) ≤ G ( s j + 1 − s j , a m j + 1 , m j , 2 k ) .   Since x ( i j )   and x ( i j + 1 )   are within distance a m j + 1 , m j , 2 k j + 1   for a fixed k j   and a fixed x ( i j + 1 )   , there could be at most 2 k j + 1 − 1 < 2 k j + 1   particles x ( i j )   of type m j   within this distance from x ( i j + 1 )   . This is the reason why the factor 2 k j + 1   appears in the previous inequality. Since x ( i j )   is of type m j   , we have y i j ≤ 2 − m j   and this is why 2 − m j   appears in ( 7.2 ).
Estimating the term G   in ( 7.2 ), with M = 2 k   , we get
G ( r , a m j + 1 , m j , 2 k )
= C r − d / 2 exp ( − a m j , m j + 1 , M 2 4 r )
= C r − d / 2 exp ( − ( 4 r ) − 1 e − 2 / d δ 8 / ( M d ) 2 − 2 p m j + 1 / ( d M ) 2 − 2 p m j / d
⋅ M 2 / d 2 − 2 δ ( m j + 1 + m j + k ) / ( M d ) )
≤ C r − d / 2 exp ( − ( 4 r ) − 1 e − 2 / d δ 8 / d
2 2 [ k ( M − δ ) − m j + 1 ( p + δ ) − m j ( p M + δ ) ] / ( d M ) ) .
Therefore, G ( r , a m j + 1 , m j , 2 k ) ≤ C r − d / 2 exp ( − C 0 r − 1 2 2 [ k ( 1 − δ ) − ( m j + m j + 1 / M ) ( p + δ ) ] / d )   where C 0   depends on δ   . Thus,
u n ( t , x ) ≤ ∫ T n ∑ k 1 , . . . , k n = 0 ∞ ∑ m 1 , . . . , m n = − N ∞
n ∏ j = 1 [ C ( s j + 1 − s j ) − d / 2 2 k j − m j
⋅ exp ( − C 0 ( s j + 1 − s j ) − 1 2 2 [ k j ( 1 − δ ) − ( m j + m j + 1 / M ) ( p + δ ) ] / d ) ] d s .
Rearranging the terms and using 2 − m 1 / ( M + 1 ) ≤ 1   and since m n + 1 = 0   , we get
u n ( t , x ) ≤ ∫ T n ∑ k 1 , . . . , k n = 0 ∞ ∑ m 1 , . . . , m n = − N ∞ (7.3)
n ∏ j = 1 [ C ( s j + 1 − s j ) − d / 2 2 k j − ( M / ( M + 1 ) ) ( m j + m j + 1 / M )
⋅ exp ( − C 0 ( s j + 1 − s j ) − 1 2 2 [ k j ( 1 − δ ) − ( m j + m j + 1 / M ) ( p + δ ) ] / d ) ] d s .
We wish to replace the above sums over m j   by integrals. This would involve replacing each variable m j   by a continuous variable w j   taking values in [ m j , m j + 1 ]   . This would decrease the exponential, but would increase the first power of 2 by a multiplicative factor of 2 n   . We would have
| u n ( t , x ) | ≤ ∫ T n ∑ k 1 , . . . , k n = 0 ∞ ∫ − N ∞ . . . ∫ − N ∞ (7.4)
n ∏ j = 1 [ C ( s j + 1 − s j ) − d / 2 2 k j − ( M / ( M + 1 ) ) ( w j + w j + 1 / M )
⋅ exp ( − C 0 ( s j + 1 − s j ) − 1 2 2 [ k j ( 1 − δ ) − ( w j + w j + 1 / M ) ( p + δ ) ] / d ) ] d w d s
where, as with d s   , we let d w = d w 1 . . . d w n   . Next, we wish to make the change of variables v j = w j + w j + 1 M , j = 1 , . . . , n .   Now w n + 1   is fixed, and the reader can easily check that the Jacobian determinant for this change of variables is 1. Since M = 2 k ≥ 1   and w j ≥ − N   , we have that v j ≥ − 2 N   . After making this change of variables, we can put the integral signs over v   inside the product as follows.
| u n ( t , x ) | ≤ ∫ T n ∑ k 1 , . . . , k n = 0 ∞ ∫ − 2 N ∞ . . . ∫ − 2 N ∞ (7.5)
n ∏ j = 1 [ C ( s j + 1 − s j ) − d / 2 2 k j − ( M / ( M + 1 ) ) v j
⋅ exp ( − C 0 ( s j + 1 − s j ) − 1 2 2 [ k j ( 1 − δ ) − v j ( p + δ ) ] / d ) ] d v d s
= ∫ T n n ∏ j = 1 [ ∑ k j = 0 ∞ ∫ − 2 N ∞ C ( s j + 1 − s j ) − d / 2 2 k j − ( M / ( M + 1 ) ) v j
⋅ exp ( − C 0 ( s j + 1 − s j ) − 1 2 2 [ k j ( 1 − δ ) − v j ( p + δ ) ] / d ) d v j ] d s
where once more, d v = d v 1 . . . d v n   . Let T j , k j   denote the term inside the final sum in ( 7.5 ). Dropping the subscripts on T , s , k   and setting r = s j + 1 − s j   , we can write
T = ∫ − 2 N ∞ C r − d / 2 2 k − ( M / ( M + 1 ) ) v (7.6)
⋅ exp ( − C 0 r − 1 2 2 [ k ( 1 − δ ) − v ( p + δ ) ] / d ) d v .
We will split the integral into two ranges for the variable v   .
For some values of v   , the exponential in ( 7.6 ) is close to 1. This will happen if the terms following C 0   in ( 7.6 ) are less than or equal to 1. We must solve r − 1 2 2 [ k ( 1 − δ ) − v ( p + δ ) ] / d ≤ 1 .   Simplifying the above expression, we get
2 − 2 v ( p + δ ) / d ≤ r 2 − 2 k ( 1 − δ ) / d (7.7)
or
2 − M M + 1 v ≤ r M d 2 ( M + 1 ) ( p + δ ) 2 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) k . (7.8)
Let Λ = Λ ( k , d , p )   be the value of v   for which equality is attained in ( 7.8 ).
We will split T   into two integrals, T = T ( 1 ) + T ( 2 )   , where
T ( 1 ) = ∫ Λ ∞ I ( v ) d v
T ( 2 ) = ∫ − 2 N Λ I ( v ) d v
and I ( v )   is the integrand in ( 7.6 ).
Case 1: An estimate for T ( 1 )   .
By ( 7.8 ), and using the fact that 1 2 ≤ M M + 1 ≤ 1   we have
T ( 1 ) = ∫ Λ ∞ C r − d / 2 2 k − ( M / ( M + 1 ) ) v exp ( − C 0 r − 1 2 2 [ k ( 1 − δ ) − v ( p + δ ) ] / d ) d v
≤ ∫ Λ ∞ C r − d / 2 2 k − ( M / ( M + 1 ) ) v d v
≤ C r − d / 2 2 k − ( M / ( M + 1 ) ) Λ
= C r − d 2 + M d 2 ( M + 1 ) ( p + δ ) 2 k ( 1 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) )
≤ C r − d 2 + d 4 ( p + δ ) 2 k ( 1 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) ) .
since 0 < r < t   . Here, C   depends on t   .
If p < 1   and if δ   is sufficiently small, then (remembering that M = 2 k   ), ∑ k = 0 ∞ 2 k ( 1 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) ) = ∑ k = 0 ∞ 2 k ( 1 − 2 k ( 1 − δ ) ( 2 k + 1 ) ( p + δ ) ) < ∞   so that
∑ k j = 0 ∞ T j , k j ( 1 ) ≤ C r − d 2 + d 4 ( p + δ ) . (7.9)
Case 2: An estimate for T ( 2 )   .
Now we turn to the case v ≤ Λ   . Here, we assume that Λ ≥ − 2 N   or there will be no values of v   for which − 2 N ≤ v ≤ Λ   . We will again drop the subscripts and set r = s j + 1 − s j   . Since equality is attained in ( 7.7 ) when v = Λ   , we have
T ( 2 )
= ∫ − 2 N Λ C r − d / 2 2 k − ( M / ( M + 1 ) ) v exp ( − C 0 r − 1 2 2 [ k ( 1 − δ ) − v ( p + δ ) ] / d ) d v
= ∫ − 2 N Λ C r − d / 2 2 k − ( M / ( M + 1 ) ) v exp ( − C 0 2 2 ( Λ − v ) ( p + δ ) / d ) d v .
Changing variables to z = Λ − v   , and since M / ( M + 1 ) < 1   , we see that
T ( 2 ) = ∫ 0 Λ + 2 N C r − d / 2 2 k − ( M / ( M + 1 ) ) ( Λ − z ) exp ( − C 0 2 2 z ( p + δ ) / d ) d z
≤ C r − d / 2 2 k − ( M / ( M + 1 ) ) Λ ∫ 0 ∞ 2 z exp ( − C 0 2 2 z ( p + δ ) / d ) d z
≤ C r − d / 2 2 k − ( M / ( M + 1 ) ) Λ
≤ C r − d 2 + d 4 ( p + δ ) 2 k ( 1 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) ) .
Once again, if p < 1   and δ   is small, then ∑ k = 0 ∞ 2 k ( 1 − M ( 1 − δ ) ( M + 1 ) ( p + δ ) ) = ∑ k = 0 ∞ 2 k ( 1 − 2 k ( 1 − δ ) ( 2 k + 1 ) ( p + δ ) ) < ∞   so that
T ( 2 ) ≤ C r − d 2 + d 4 ( p + δ ) . (7.10)
Combining ( 7.9 ) and ( 7.10 ), if p < 1   and δ   is sufficiently small, then T = T ( 1 ) + T ( 2 ) ≤ C r − d 2 + d 4 ( p + δ )   and u n ( t , x ) ≤ C ∫ T n n ∏ i = 1 ( s i + 1 − s i ) − d 2 + d 4 ( p + δ ) d s .   Now, a calculation of Hu ([Hu01] , proof of Theorem 4.1, equation (4.10)) gives, for β > 0   ,
∫ T n n ∏ k = 1 ( s k + 1 − s k ) − β d s = t n ( 1 − β ) Γ n ( 1 − β ) Γ ( 1 + n ( 1 − β ) ) . (7.11)
We give a quick derivation of Hu's estimate in the appendix, since he merely cites the result.
In ( 7.11 ), substitute β = d 2 − d 4 ( p + δ ) .   If β < 1   , then Γ ( 1 − β )   is defined. With C n   denoting a term of order C n   for some C > 0   , we get u n ( t , x ) ≤ t n ( 1 − β ) Γ n ( 1 − β ) Γ ( 1 + n ( 1 − β ) ) ≤ C n t n ( 1 − β ) n − n ( 1 − β )   since only the denominator has order greater than C n   . Thus, if β < 1   , and if A   and A 1   do not occur, u ( t , x ) = ∑ n = 0 ∞ u n ( t , x ) < ∞ .   Solving for β < 1   , we obtain the requirement that d 2 − d 4 ( p + δ ) < 1 .   For small enough values of δ   , this is satisfied if either d ≤ 4   or p < 1 2 + 1 d − 2 , d ≥ 5 .   This completes the proof of Theorem  1 , for L K   . As mentioned earlier, since P ( A K ) = P ( L = L K ) → 1   as K → ∞   , it follows that Theorem  1 holds almost surely.

8 Proof of the existence part of Theorem  2 ( p ≥ 1   )

For p ≥ 1   , we also wish to replace L   by L K   . But in this case our stochastic integrals will no longer be martingales unless we delete the compensator for L − L K   . We define this compensator as follows.
ℒ ( d t d x ) = d t d x ∫ K ∞ c p y − ( p + 1 ) d y = c p ′ K − p d t d x .   Then, on the event A K   , we can replace ( 1.1 ) by
∂ u ∂ t = Δ u + c p ′ K − p u + u ⋄ L ˙ K ( x )
u ( 0 , x ) = u 0 ( x ) .
If we let w ( t , x ) = exp ( − c p ′ K − p t ) u ( t , x )   then w ( t , x )   solves
∂ u ∂ t = Δ u + u ⋄ L ˙ K ( x ) (8.1)
u ( 0 , x ) = u 0 ( x ) .
Therefore, any conclusions which hold almost surely for w ( t , x )   , for all K   , must hold for solutions to ( 1.1 ).
For ease of notation, we write u ( t , x )   instead of w ( t , x )   and L   instead of L K   . With this notation, we can prove the following lemma.
Lemma 3 For any q ∈ ( p , 2 ]   there exists a positive constant C q   such that for any measurable symmetric function f : D n → R   we have:
E [ | n ! ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ] (8.2)
≤ C q n ( n ! ) q − 1 ∫ D n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n ) .
Note that the above lemma gives sufficient conditions for the existence of the multiple stochastic integral I n ( f n )   .
Proof. We would like to use the Burkholder-Davis-Gundy inequality, and for this reason we need to expand on our formulation of the multiple stochastic integral given in Section ( 5 ), so that we have a martingale. We start with the σ   -algebra. Recall that we have chosen the first coordinate of x   as our time variable. In particular, we can define a filtration. Let m ( D ) , M ( D )   be the minimum and maximum values of x 1   , the first coordinate of x   , for x ∈ D   . For z ∈ [ m ( D ) , M ( D ) ]   , we define the σ   -algebra ℋ z   as follows.
ℋ z = σ { L ( A ) : A ⊂ { x ∈ D : x 1 ≤ z } } .   Next, let R ( z , D ) = { x ∈ D : x 1 ≤ z } .   Let h ( x )   be a predictable function with respect to the ℋ z   filtration. Let Y z = ∫ R ( z , D ) h ( x ) L ( d x ) .   Then ( Y z , ℋ z )   is a martingale, for appropriate conditions on h   , and its quadratic variation is 〈 Y 〉 z = ∑ x 1 ( i ) ≤ z h 2 ( x ( i ) ) y i 2 .   For ease of notation, we drop the subscript z   when z = M ( D )   , so that
Y = Y M ( D ) ,
〈 Y 〉 = 〈 Y 〉 M ( D ) .
The Burkholder-Davis-Gundy inequality now implies that for q > 1   , there is a constant C q   depending only on q   , such that E [ | Y | q ] ≤ C q E [ 〈 Y 〉 q / 2 ] .   Recall that for r ≤ 1   and a sequence of nonnegative numbers a i   , the following elementary inequality holds.
( ∑ i = 1 ∞ a i ) r ≤ ∑ i = 1 ∞ a i r .   Setting r = q / 2   , and using our formula for 〈 Y 〉   , we get
E [ | Y | q ] ≤ C E [ ( ∑ i = 1 ∞ h 2 ( x ( i ) ) y i 2 ) q / 2 ] (8.3)
≤ C E [ ∑ i = 1 ∞ | h 2 ( x ( i ) ) | q / 2 y i q ]
= C ∫ D E | h ( x ) | q d x
provided q > p   . The constant C   depends on p   , q   , and K   . Applying ( 8.3 ) repeatedly, we get
E [ | ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ] (8.4)
≤ C n ∫ D ~ n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n ) .
Note, that the repeated application of inequality ( 8.3 ) was possible due to the order x 1 ( 1 ) < x 1 ( 2 ) < . . . < x 1 ( n )   that we imposed on D ~ n   . To see this, let us remember that the function h   involved in ( 8.3 ) needed to be predictable.
Thus,
E [ | ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ]
= E [ | ∫ D h ( x ( n ) ) L ( d x ( n ) ) | q ] ,
where h ( x ( n ) ) = ∫ R ~ n − 1 ( x 1 ( n ) ) f ( x ( 1 ) , . . . , x ( n − 1 ) , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n − 1 ) )   and
R ~ n − 1 ( y ) = { ( x ( 1 ) , . . . , x ( n − 1 ) ) ∈ D ~ n − 1 : x 1 ( n − 1 ) < y } (8.5)
where y ∈ R   . The function h ( x ( n ) )   is predictable since x 1 ( 1 ) < x 1 ( n )   , . . .   , x 1 ( n − 1 ) < x 1 ( n )   on R ~ ( x 1 ( n ) , D )   . So, we can apply the inequality ( 8.3 ) and obtain
E [ | ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ]
≤ C ∫ D E [ | h ( x ( n ) ) | q ] d x ( n ) ,
and so on. Using ( 8.3 ), and the fact that f   is symmetric, we get:
E [ | n ! ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ]
= ( n ! ) q E [ | ∫ D ~ n f ( x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) | q ]
≤ ( n ! ) q C n ∫ D ~ n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n )
= ( n ! ) q C n 1 n ! ∫ D n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n )
= C n ( n ! ) q − 1 ∫ D n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n ) .
Before proving Theorem  2 , we present two simple lemmas.
Lemma 4 Let n   be a positive integer and q   a real number, such that q ≥ 1   . Let f : D n → R   be a function in L q ( D n , d x )   , where d x   is Lebesgue measure. If f ~   is the symmetrization of f   , then f ~ ∈ L q ( D n , d x )   and
∥ f ~ ∥ q ≤ ∥ f ∥ q (8.6)
where ∥ ⋅ ∥ q   denotes the norm of the space L q ( D n , d x )   .
Proof. For any permutation σ   of the set { 1 , 2 , . . . , n }   , we can consider the function f σ : D n → R   , defined by f σ ( x ( 1 ) , . . . , x ( n ) ) = f ( x ( σ ( 1 ) ) , . . . , x ( σ ( n ) ) ) .   Making the change of variable y ( 1 ) = x ( σ ( 1 ) )   , . . .   , y ( n ) = x ( σ ( n ) )   , we can see that
∫ D n | f σ ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n )
= ∫ D n | f ( x ( 1 ) , . . . , x ( n ) ) | q d x ( 1 ) . . . d x ( n ) .
Thus for all permutations σ   , we have
∥ f σ ∥ q = ∥ f ∥ q .
Since the symmetrization f ~   of f   is defined as f ~ = 1 n ! ∑ σ f σ ,   we can apply the triangle inequality to get
∥ f σ ∥ q = ∥ 1 n ! ∑ σ f σ ∥ q
≤ 1 n ! ∑ σ ∥ f σ ∥ q
= 1 n ! ∑ σ ∥ f ∥ q
= ∥ f ∥ q .
In the proof of Lemma  1 we used the inequality e n ≥ n n / n !   , for all positive integers n   . We would like to prove this inequality for all positive real numbers x   and therefore we have to replace n !   by Γ ( x + 1 )   . Thus we will prove the following lemma.
Lemma 5 For any positive real number x   , we have
Γ ( x + 1 ) ≥ x x e − x . (8.7)
Proof. If x > 0   , then for all t ≥ x   , t x ≥ x x   . Thus we have:
Γ ( x + 1 ) = ∫ 0 ∞ t x e − t d t
≥ ∫ x ∞ t x e − t d t
≥ ∫ x ∞ x x e − t d t
= x x ∫ x ∞ e − t d t
= x x e − x .
We can prove now Theorem  2 . According to inequalities ( 8.2 ) and ( 8.6 ), we have:
E [ | u n ( t , x ) | q ] (8.8)
≤ C n ( n ! ) q − 1 ∫ D n [ s y m x ( 1 ) , . . . , x ( n ) | u 0 ( x ( 1 ) ) |
⋅ ∫ T n n ∏ k = 1 G ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s ] q d x
≤ a C n ( n ! ) q − 1 ∫ D n [ ∫ T n n ∏ k = 1 G ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s ] q d x ,
where a = ∥ u 0 ∥ ∞   and T n = { ( s 1 , . . . , s n ) ∈ R n : 0 < s 1 < s 2 < . . . < s n < t } .   We may assume that ∥ u 0 ∥ ∞ = 1   . Let V ( T n ) = t n / n !   be the volume of T n   .
Using Jensen's inequality, we obtain
[ ∫ T n n ∏ k = 1 G ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s ] q (8.9)
= V ( T n ) q [ V ( T n ) − 1 ∫ T n n ∏ k = 1 G ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s ] q
≤ V ( T n ) q V ( T n ) − 1 ∫ T n n ∏ k = 1 G q ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s
= ( t n n ! ) q − 1 ∫ T n n ∏ k = 1 G q ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d s .
Recall that
∫ R d G q ( r , x ) d x = C r − ( q − 1 ) d / 2 ∫ R d r − d / 2 exp ( − q x 2 4 r ) d x
= C r − ( q − 1 ) d / 2 . (8.10)
Combining ( 8.8 ) and ( 8.9 ), we see that ( n ! ) q − 1   is canceled out by ( 1 / n ! ) q − 1   , and by using Fubini's theorem and ( 8.10 ) repeatedly (i.e. integrated first over x ( 1 )   , then x ( 2 )   , and so on), we obtain:
E [ | u n ( t , x ) | q ] (8.11)
≤ ( C t q − 1 ) n ∫ T n [ ∫ D n n ∏ k = 1 G q ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d x ] d s
≤ ( C t q − 1 ) n ∫ T n [ ∫ R n n ∏ k = 1 G q ( s k + 1 − s k , x ( k + 1 ) − x ( k ) ) d x ] d s
= ( C t q − 1 ) n ∫ T n n ∏ k = 1 ( s k + 1 − s k ) − ( q − 1 ) d / 2 d s .
Now, an estimate of Hu ( 7.11 ) gives ∫ T n n ∏ k = 1 ( s k + 1 − s k ) − ( q − 1 ) d / 2 d s = t n ( 1 − ( d / 2 ) ( q − 1 ) ) Γ n ( 1 − d 2 ( q − 1 ) ) Γ ( 1 + n ( 1 − d 2 ( q − 1 ) ) ) .   If d < 2 / ( q − 1 )   , then Γ ( 1 − ( d / 2 ) ( q − 1 ) )   is defined. For d < 2 / ( q − 1 )   , let α = 1 − ( d / 2 ) ( q − 1 ) > 0   . Then from ( 8.11 ) we conclude that
E [ | u n ( t , x ) | q ] ≤ ρ n Γ ( 1 + n α ) , (8.12)
for some constant ρ   depending on t   , q   , p   , and K   . Using the well-known inequality ( b + c ) q ≤ 2 q − 1 ( b q + c q )   , for all q ≥ 1   and b   and c   non-negative numbers, an easy induction argument shows that for a n ≥ 0   , and c = 2 q − 1   , ( ∑ n = 0 ∞ a n ) q ≤ ∑ n = 0 ∞ c n + 1 a n q .   Applying this inequality to u ( t , x )   , we find
E [ | u ( t , x ) | q ] = E [ ( ∑ n = 0 ∞ | u n ( t , x ) | ) q ] (8.13)
≤ E [ ∑ n = 0 ∞ c n + 1 | u n ( t , x ) | q ]
= ∑ n = 0 ∞ c n + 1 E [ | u n ( t , x ) | q ]
≤ c ∑ n = 0 ∞ 1 Γ ( 1 + n α ) ( c ρ ) n .
Using now ( 8.7 ) we obtain:
E [ | u ( t , x ) | q ] ≤ c ∑ n = 0 ∞ e n α ( n α ) n α ( c ρ ) n . (8.14)
Thus, we can see from ( 8.14 ) that if d < 2 / ( q − 1 )   (or equivalently q < 1 + ( 2 / d )   ) and p < q   , then E [ | u ( t , x ) | q ] < ∞   . (The condition p < q   was used in the Lemma  3 .) For such a q   to exist we need that p < 1 + 2 / d   .
Let Π   be the random Poisson set { ( x ( i ) , y i ) } i ≥ 1   . For any positive integer k   we define the event:
A k = { ω ∈ Ω : ∀ ( x , y ) ∈ Π ( ω ) , y ≤ k } .   That means, A k   is the event that all atoms have masses less than or equal to k   .
Thus if p < 1 + ( 2 / d )   , then for all q ∈ ( p , 1 + ( 2 / d ) )   and for all k ≥ 1   , E [ | u ( t , x ) | 1 A k q ] < ∞ .   Since, the expectation is the integration with respect to a probability (finite) measure, it follows, that for any q ′   and q   , such that 1 ≤ q ′ ≤ p < q < 1 + ( 2 / d )   , we have:
( E [ | u ( t , x ) 1 A k | q ′ ] ) 1 / q ′ ≤ ( E [ | u ( t , x ) 1 A k | q ] ) 1 / q < ∞ .   Thus the solution u ( t , x )   exists, for all t ≥ 0   and all x ∈ D   , and belongs to the space L l o c q ( Ω )   , for all 1 ≤ q < 1 + ( 2 / d )   .
This completes the proof of Theorem  2 , if u ( t , x )   means w ( t , x )   and L   means L K   . As mentioned at the beginning of Section  8 , this proof carries over to u ( t , x )   and L   , since P ( A K ) = P ( L = L K ) → 1   as K → ∞   .

9 Uniqueness

9.1 The Main Lemma

We consider both cases p < 1   and p ≥ 1   . Suppose u ( t , x )   is a solution to ( 1.1 ) in case p ≥ 1   or ( 1.2 ) in case p < 1   . Then we have
u ( t , x ) = ∑ n = 0 ∞ g n ( t , x ) (9.1)
where g n ∈ J n   in case p < 1   or g n ∈ ℐ n   in case p ≥ 1   .
Lemma 6 Let u ( t , x )   be as in ( 9.1 ). If the g n   are uniquely determined from u   , then uniqueness holds for ( 1.1 ) or ( 1.2 ).
Proof. Both cases have the same proof, so we will only deal with the case p ≥ 1   . Suppose there are two solutions u i ( t , x )   for i = 1 , 2   , where u i ( t , x ) = ∑ n = 0 ∞ g n ( i ) ( t , x )   and g n ( i ) ( t , x ) ∈ ℐ n   . Let Y ( t , x ) = u 1 ( t , x ) − u 2 ( t , x ) = ∑ n = 0 ∞ g n ( t , x )   where g n ( t , x ) = g n ( 1 ) ( t , x ) − g n ( 2 ) ( t , x )   . Then Y ( t , x )   is also a solution to ( 1.1 ) with initial condition Y ( 0 , x ) = 0   , so by ( 6.2 ), we have
∑ n = 0 ∞ g n ( t , x ) = Y ( t , x )
= ∫ D ( ∫ 0 t G ( t − s , x , y ) Y ( s , y ) d s ) ⋄ L ( d y )
= ∑ n = 0 ∞ ∫ D ( ∫ 0 t G ( t − s , x , y ) g n ( s , y ) d s ) ⋄ L ( d y )
= ∑ n = 1 ∞ g ~ n ( t , x )
∀ x ∈ R d , t ∈ ( 0 , T ] .
where
g ~ n ( t , x ) = ∫ D ( ∫ 0 t G ( t − s , x , y ) g n − 1 ( s , y ) d s ) ⋄ L ( d y ) (9.2)
∈ ℐ n .
Since we are assuming that Y ( t , x )   uniquely determines both g n ( t , x )   and g ~ n ( t , x )   , it follows that g 0 ( t , x ) = 0   and
g n ( t , x ) = g ~ n ( t , x ) (9.3)
for any positive value of n   . Thus, ( 9.2 ) implies that g ~ 1 ( t , x ) = 0   , and hence by ( 9.3 ) g 1 ( t , x ) = 0   . An induction argument now yields g n ( t , x ) = 0   for all values of n   .

9.2 Proof of Uniqueness for p < 1  

For simplicity in this proof we will assume from the beginning that L ˙   is a stable noise (not L ˙ K   ) without truncation of large atoms.
We begin by ordering the atoms of L ˙   according to size, so that y 1 > y 2 > ⋯   .
For z ∈ ( 1 , 2 )   , let R i ( z ) L ˙   be the transformation of the noise L ˙   which replaces y i   with z y i   .
Lemma 7 Fix z ∈ ( 1 , 2 ) , i ≥ 1   . Let M ˙ = R i ( z ) L ˙   Also, let P M , P L   be the probability measures on the canonical probability space of atomic measures, generated by M ˙   and L ˙   , respectively. Then P M < < P L .  
Proof. For any k ≥ 1 , n ≥ 0   define an event
A k , n ≡ { there are exactly n atoms of L ˙
with masses greater than or equal to 1 / k }
Note that for any k ≥ 1  
∞ ⋃ n = 0 A k , n = Ω
A k , l ⋂ A k , m = ∅ , ∀ m ≠ l .
It is obvious that
lim k → ∞ P L ( i ⋃ l = 0 A k , l ) = 0 . (9.4)
Fix δ > 0   arbitrary small. By ( 9.4 ) we can fix k   sufficiently large such that
P L ( i ⋃ l = 0 A k , l ) ≤ δ . (9.5)
Define
A ¯ k , i ≡ ( i ⋃ l = 0 A k , l ) c = ∞ ⋃ l = i + 1 A k , l . (9.6)
Now will show that
P M < < P L on A ¯ k , i . (9.7)
Let us decompose L ˙   as follows:
L ˙ = L ˙ 1 , k + L ˙ 2 , k   where L ˙ 1   includes atoms of size y ≥ 1 / k   and L ˙ 2   includes atoms of size y < 1 / k   . Note that L ˙ 1 , k , L ˙ 2 , k   are independent, and
R i ( z ) L ˙ = R i ( z ) L ˙ 1 , k + L ˙ 2 , k on A ¯ k , i .
Hence, if we define M ˙ 1 = R i ( z ) L ˙ 1 , k   then, in order to get ( 9.7 ), it suffices to show that
P M 1 < < P L 1 , k on A k , l , ∀ l ≥ i + 1 (9.8)
where P M 1 , P L 1 , k   are the measures induced by M ˙ 1 , L ˙ 1 , k   respectively. Fix an arbitrary l ≥ i + 1   . On A k , l   , M ˙ 1   and L ˙ 1 , k   have l   atoms at the same locations. Therefore, to get the desired absolute continuity we just need to show the absolute continuity of the laws of the corresponding atom sizes. To be more precise, from now on we assume that A k , l   occurs. Let Y ^ 1 L , … , Y ^ l L   (resp. Y ^ 1 M , … , Y ^ l M   ) be the unordered sizes of the atoms of L 1 , k   (resp.
M 1   ). ( Y ^ 1 L , … , Y ^ l L )   and ( Y ^ 1 M , … , Y ^ l M )   are continuous l   -dimensional random variables whose laws are supported on ( 1 / k , ∞ ) l   . Moreover, the probability density function of ( Y ^ 1 L , … , Y ^ l L )   is positive at any point of ( 1 / k , ∞ ) l   . This immediately implies that conditioned on A k , l   the probability law of ( Y ^ 1 M , … , Y ^ l M )   is absolutely continuous with respect to the law of ( Y ^ 1 L , … , Y ^ l L )   , and since l ≥ i + 1   was arbitrary, ( 9.8 ) and ( 9.7 ) follow.
The end of the proof of the lemma is trivial. Let B   be any event such that P L ( B ) = 0 .   Then
P M ( B ) = ∑ l = 0 i P M ( B ⋂ A k , l ) + ∑ l = i + 1 ∞ P M ( B ⋂ A k , l )
= ∑ l = 0 i P M ( B ⋂ A k , l )
≤ P M ( i ⋃ l = 0 A k , l ) , (9.9)
where the second equality follows from ( 9.7 ). Now it is easy to check the following equality of events:
{ there are no more than i atoms of L ˙ with masses y ≥ 1 / k }
= { there are no more than i atoms of M ˙ with masses y ≥ 1 / k }
Therefore, we get P M ( i ⋃ l = 0 A k , l ) = P L ( i ⋃ l = 0 A k , l ) .   Thus, ( 9.5 ) and ( 9.9 ) imply that P M ( B ) ≤ δ .   Since δ > 0   was arbitrary, we get P M ( B ) = 0   and the result follows.
By Lemma  6 the following lemma gives the required uniqueness for the case p < 1   .
Lemma 8 Let p < 1   . In ( 9.1 ), the functions g n   are uniquely determined.
Proof We argue by contradiction. Suppose Lemma  8 were false. Then by subtracting the two series, we form a series G ≡ ∑ n = 0 ∞ g n = 0   where g n ∈ J n   , and at least one of the g n   is not identically 0. In fact, we may assume that g n = ∑ ( i 1 , . . . , i n ) ∈ J n y i 1 . . . y i n h n ( x i 1 , . . . , x i n )   with
∑ n = 0 ∞ ∑ ( i 1 , . . . , i n ) ∈ J n y i 1 . . . y i n | h n ( x i 1 , . . . , x i n ) | < ∞ (9.10)
with probability 1. Actually, each solution is a difference of two series for which the functions h n   are nonnegative.
Recall that { y 1 > y 2 > … }   are the ordered sizes of the atoms of our noise L ˙   . Fix an arbitrary n ≥ 1   . Our goal is to show that, almost surely, y 1 . . . y n h n ( x 1 , . . . , x n ) = 0   from which it follows that h n = 0   almost everywhere. By Lemma  7 , for z i ∈ ( 1 , 2 )   , the probability induced by R i ( z i ) L ˙   is absolutely continuous with respect to the probability induced by L ˙   . Furthermore, R i ( z i )   extends to an operator on J   , as follows. We set R i ( z i ) ∑ n = 0 ∞ g n ≡ ∑ n = 0 ∞ R i ( z i ) g n ,   where we define R i ( z i ) g n   to be a multiple stochastic integral with L ˙   replaced by R i ( z i ) L ˙   . Therefore, if G = ∑ n = 0 ∞ g n = 0   almost surely, then R i ( z i ) ∑ n = 0 ∞ g n ≡ ∑ n = 0 ∞ R i ( z i ) g n = 0   almost surely as well. Furthermore, if z i ∈ ( 1 , 2 )   , then R 1 ( z 1 ) . . . R m ( z m ) G   is an analytic function in z 1 , . . . , z m   . Let Λ n , m G   be the terms in the expansion of this analytic function which contain z 1 . . . z n   but no other z i   . Note that each Λ n , m G   equals 0 almost surely, and that lim m → ∞ Λ n , m G = y 1 . . . y n h n ( x 1 , . . . , x n ) = 0   almost surely. In fact, the convergence to 0 follows from ( 9.10 ), since the operator Λ n , m   removes more and more terms as m   increases. This proves Lemma  8 .

9.3 Proof of Uniqueness for p ≥ 1  

By Lemma  6 the following lemma gives the required uniqueness for the case p ≥ 1   .
Lemma 9 Let p ≥ 1   . In ( 9.1 ), the functions g n   are uniquely determined.
Proof. Fix ( t , x )   and suppose that u 1 ( t , x ) , u 2 ( t , x )   are two solutions.
Subtracting, we obtain
Y ( t , x ) = u 1 ( t , x ) − u 2 ( t , x )
= ∑ n = 0 ∞ g n ( 1 ) ( t , x ) − ∑ n = 0 ∞ g n ( 2 ) ( t , x )
= : ∑ n = 0 ∞ g ~ n ( t , x )
where g n ( i ) ( t , x )   corresponds to the solution u i ( t , x )   . Also, each g ~ n   can be expressed as a stochastic integral g ~ n ( t , x ) = ∫ D ~ n f n ( t , x ; x ( 1 ) , . . . , x ( n ) ) L ( d x ( 1 ) ) . . . L ( d x ( n ) ) .   Adding these terms, and using the notation from the beginning of Section  8 , we find that Y ( t , x ) = f 0 ( t , x ) + ∫ R ( M ( D ) , D ) f ~ 1 ( t , x ; y ) L ( d y )   where
f ~ 1 ( t , x ; y ) = f 1 ( t , x ; y )
+ ∑ n = 2 ∞ ∫ R ~ n − 1 ( y ) f n ( t , x ; x ( 1 ) , . . . , x ( n − 1 ) , y ) L ( d x ( 1 ) ) . . . L ( d x ( n − 1 ) )
and R ~ n − 1 ( y )   was defined in ( 8.5 ).
Since Y ( t , x ) = 0   and z → ∫ R ( z , D ) f ~ 1 ( t , x ; y ) L ( d y )   is a local ℋ z   martingale, it follows that f 0 ( t , x ) = 0   and that the integrand f ~ 1 ( t , x ; y ) = 0   for almost every value of y   . But f ~ 1 ( t , x ; y ) = 0   is itself a sum of a deterministic function and a stochastic integral. Therefore we may use the same argument to show that f 1 ( t , x ; y ) = 0   for almost every value of y   . Continuing, we can use induction to show that for each value of n   , f n ( t , x ; y 1 , . . . , y n ) = 0   for almost every value of y 1 , . . . , y n   .

A A Derivation of Hu's Estimate

In the proof of Theorem 4.1, equation (4.10) of [Hu02] , Hu gives the following estimate (in his equation, α = 1 − d / 4   , and we use β = d / 4   ).
ℋ n ( t ) : = ∫ T n ( t ) n ∏ k = 1 ( s k + 1 − s k ) − β d s .   Then
ℋ n ( t ) = t n ( 1 − β ) Γ ( 1 − β ) n Γ ( 1 + n ( 1 − β ) ) . (A.1)
(Hu uses β = γ / 4   . Also, Hu has 2 instead of 1 in the above formula. His answer is equivalent if we ignore multiplicative factors of n.) Making the change of variables s = α r   , d s = α d r   , we find
ℋ n ( α t ) = ∫ T n ( α t ) n ∏ k = 1 ( s k + 1 − s k ) − β d s
= α − n β α n ℋ n ( t )
= α n ( 1 − β ) ℋ n ( t ) .
Then we get the recurrence
ℋ n + 1 ( t ) = ∫ 0 t ( t − s ) − β ℋ n ( s ) d s
= ∫ 0 t ( t − s ) − β s n ( 1 − β ) ℋ n ( 1 ) d s
= ℋ n ( 1 ) ∫ 0 t ( t − s ) − β s n ( 1 − β ) d s
= t − n ( 1 − β ) ℋ n ( t ) ∫ 0 t ( t − s ) − β s n ( 1 − β ) d s .
Let
a = β
b = n ( 1 − β ) .
Recall that the Euler beta function is defined as follows, for a , b > 0   .
B ( a , b ) : = ∫ 0 1 ( 1 − x ) a − 1 x b − 1 d x .   It is a fundamental identity that B ( a , b ) = Γ ( a ) Γ ( b ) Γ ( a + b ) .   Changing variables, we get ∫ 0 t ( t − s ) − a s b d s = t 1 − a + b Γ ( 1 − a ) Γ ( 1 + b ) Γ ( 2 − a + b ) .   Therefore,
ℋ n + 1 ( t ) = t − n ( 1 − β ) ℋ n ( t ) t 1 − β + n ( 1 − β ) Γ ( 1 − β ) Γ ( 1 + n ( 1 − β ) ) Γ ( 2 − β + n ( 1 − β ) )
= t ( 1 − β ) ℋ n ( t ) Γ ( 1 − β ) Γ ( 1 + n ( 1 − β ) ) Γ ( 1 + ( n + 1 ) ( 1 − β ) )
The product telescopes. Since ℋ 1 ( t ) = t 1 − β   , we get ( A.1 ).
References

  1. S. Albeverio, F. Gesztesy, R. Høegh-Krohn, and H. Holden. Solvable models in quantum mechanics. Texts and Monographs in Physics. Springer-Verlag, New York, 1988.
  2. J. Bertoin. Lévy processes, volume 121 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 1996.
  3. D.A. Dawson. Measure-valued Markov processes. In École d'été de probabilités de Saint-Flour, XXI-1991, P. L. Hennequin, editor, number 1180 in Lecture Notes in Mathematics, pages 1–260, Berlin, Heidelberg, New York, 1993. Springer-Verlag.
  4. Y. Hu. Heat equations with fractional white noise potentials. Appl. Math. Optim., 43(3):221–243, 2001.
  5. Y. Hu. Chaos expansion of heat equations with white noise potentials. Potential Anal., 16(1):45–66, 2002.
  6. S. Janson. Gaussian Hilbert spaces, volume 129 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 1997.
  7. G. Kallianpur, J. Xiong, G. Hardy, and S. Ramasubramanian. The existence and uniqueness of solutions of nuclear space-valued stochastic differential equations driven by Poisson random measures. Stochastics Stochastics Rep., 50(1-2):85–122, 1994.
  8. E. Lytvynov. Orthogonal decompositions for Lévy processes with an application to the gamma, Pascal, and Meixner processes. Infin. Dimens. Anal. Quantum Probab. Relat. Top., 6(1):73–102, 2003.
  9. R. Mikulevicius and B. Rozovskii. Linear parabolic stochastic PDEs and Wiener chaos. SIAM J. Math. Anal., 29(2):452–480 (electronic), 1998.
  10. C. Mueller. The heat equation with Lévy noise. Stochastic Process. Appl., 74(1):67–82, 1998.
  11. L. Mytnik. Stochastic partial differential equation driven by stable noise. Probab. Theory Related Fields, 123(2):157–201, 2002.
  12. D. Nualart and B. Rozovskii. Weighted stochastic Sobolev spaces and bilinear SPDEs driven by space-time white noise. J. Funct. Anal., 149(1):200–225, 1997.
  13. D. Nualart and M. Zakai. Generalized Brownian functionals and the solution to a stochastic partial differential equation. J. Funct. Anal., 84:279–296, 1989.
  14. E. Pardoux. Stochastic partial differential equations, a review. Bull. Sc. Math., 117:29–47, 1993.
  15. J. Rosiński and W.A. Woyczyński. On Itô stochastic integration with respect to p   -stable motion inner clock, integrability of sample paths, double and multiple integrals. Ann. Probab., 14:271–286, 1986.

1   Supported by an NSF grant.

2   Research supported in part by the Israel Science Foundation (grant No. 116/01-10.0).
Key words and phrases heat equation, white noise, stochastic partial differential equations.
AMS 1991 subject classifications Primary, 60H15; Secondary, 35R60, 35K05.