International Journal of Partial Differential Equations

Volume 2013 (2013), Article ID 297857, 7 pages

http://dx.doi.org/10.1155/2013/297857

## Probabilistic Representations for the Solution of Higher Order Differential Equations

^{1}CIRM Centro Internazionale per la Ricerca Matematica, Fondazione Bruno Kessler, 38123 Trento, Italy^{2}Dipartimento di Matematica, Università di Trento, Via Sommarive 14, 38123 Trento, Italy

Received 16 April 2013; Revised 12 June 2013; Accepted 13 June 2013

Academic Editor: David Holcman

Copyright © 2013 S. Mazzucchi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A probabilistic representation for the solution of the partial differential equation , is constructed in terms of the expectation with respect to the measure associated to a complex-valued stochastic process.

#### 1. Introduction

The connection between the solution of parabolic equations associated to second-order elliptic operators and the theory of stochastic processes is a largely studied topic [1]. The main instance is the * Feynman-Kac formula*, providing a representation of the solution of the heat equation with potential (the continuous functions vanishing at infinity):
in terms of an integral with respect to the measure of the Wiener process, the mathematical model of the Brownian motion [2]:

If the Laplacian in (1) is replaced with a higher order differential operator, that is, if we consider a Cauchy problem of the form with , , then a formula analogous to (2), giving the solution of (3) in terms of the expectation with respect to the measure associated to a Markov process, is lacking. In fact, such a formula cannot be proved for semigroups whose generator does not satisfy the maximum principle, as in the case of with . In fact in the case where the Cauchy problem (3) is not well posed on the space of continuous bounded functions [3]. In other words it is not possible to find a stochastic process which plays for the parabolic equation (3) the same role that the Wiener process plays for the heat equation.

We would like to point out that the problem of the probabilistic representation of the solution of the Cauchy problem (3) presents some similarities with the problem of the mathematical definition of Feynman path integrals (see [4–7] for a discussion of this topic). Indeed in both cases it is not possible to implement an integration theory of Lebesgue type in terms of a bounded variation measure on a space of continuous paths [8]. An analogous of the Feynman-Kac formula for the parabolic equation (3), namely, an equation of the form (where should be some “measure” on a space of “paths” ), can be obtained only under some restrictions on and and by giving up a traditional integration theory in the Lebesgue sense with respect to a bounded variation measure on a space of (real) continuous paths.

In the mathematical literature two main approaches have been proposed. The first one [9, 10] realizes formula (4) in terms of the expectation with respect to a * signed* measure on a space of paths on the interval .

It is worthwhile to mention that an analogous of the arc-sine law [10, 11], of the central limit theorem [12], and of Ito formula and Ito stochastic calculus [10, 13] have been proved for the (finite additive) signed measure . Moreover, a Feynman-Kac formula has been proved [9–11], for the representation of the solution of the Cauchy problem (3) in the case where is a bounded piecewise continuous function and for an initial datum , by realizing formula (4) as limit of finite dimensional cylindrical approximations [14]. We also mention the work by Levin and Lyons [15] on rough paths, conjecturing that the signed measure could exist on the quotient space of equivalence classes of paths corresponding to different parametrization of the same path.

A different approach is based on the construction of a stochastic process on a space of * complex paths*. In this case the integration is performed with respect to a well-defined positive probability measure on a complex space. One of the first results was given by Funaki [16], who constructed a complex stochastic process by composing two independent Brownian motions and in the following way:
and proving that, for a suitable class of analytic initial datum , the solution of the Cauchy problem
is given by the expectation

In fact this result can be generalized to partial differential equations of order , by multiple iterations of suitable processes [16–18]. These results are also related to Bochner subordination [19].

There are also similarities between the Funaki’s process and the “iterated Brownian motion” [20], but the latter is not connected to the probabilistic representation of the solution of a partial differential equation with regular coefficients. In fact the processes constructed by iterating copies of independent BMs (or other process) are associated to higher order PDE of particular form, where the initial datum plays a particular role and enters also in the differential equation [21].

Complex-valued processes, connected to PDE of the form (3), have been also proposed by other authors by means of different techniques. In [22–24] Madrecki et al. consider the fourth degree heat-type equation (6) and construct a probabilistic representation for its solution in terms of a stable probabilistic Borel measure on the space of continuous mappings on with values in the set of complex-valued sequences, endowed with the product topology. In this setting a Feynman-Kac type formula is proved, for the fourth order heat equation with linear potential.

Another probabilistic approach is presented by Sainty in [25], where a representation for the solution of is given in terms of the expectation with respect to a particular complex-valued process , called “Brownian motion of order ”. It is worthwhile also to mention a completely different approach proposed by Léandre [26, 27], which has some analogies with the mathematical realization of Feynman path integrals by means of white noise calculus [28]. Indeed Léandre has recently constructed a “probabilistic representation” of the solution of the Cauchy problem (3) in terms of an infinite dimensional distribution on the Connes space [27, 29, 30]. Another interesting approach related to the theory of pseudoprocesses introduced by Daletsky and Fomin [31] has been recently proposed by Smorodina and Faddeev [32].

We eventually mention an interesting probabilistic approach to the equation described in [33].

The present paper presents the construction of an alternative complex-valued stochastic process generalizing Funaki's result [16] and a corresponding probabilistic representation for the solution of the Cauchy problem: with and for satisfying suitable analyticity assumptions.

The paper is organized as follows. Section 2 presents the construction of a complex random variable and the representation of the solution of (8) with in terms of the expectation with respect to the probability measure associated with . Section 3 presents the proof of a Feynman-Kac type formula for the solution of (8) in the cases where is linear in the space variables and presents an explicit time dependence.

#### 2. A Complex-Valued Random Variable Associated to the 4-Order Heat-Type Equation

In the present section we construct a probabilistic representation for the solution of (8) in the case where , namely, An equation of this form in the case where , as mentioned in the introduction, has been studied by several authors by means of different techniques [9, 10, 16, 25]. In this section we show that the results in [9, 10, 16] can be seen as particular cases of a general theory presented in [34–37], connecting the solution of parabolic problems with the solution of related hyperbolic problems.

Given a Banach space and a strongly continuous group of operators on with generator , it is possible to construct the holomorphic semigroup with generator in terms of a Gaussian expectation of the group : More generally, given a polynomial in with complex coefficients, whose leading term has the form , with , then generates a holomorphic semigroup on . Its action on a vector belonging to the domain of is given by with (see [36]). Equation (11) can be regarded as a particular example of the “method of transmutation” described in [38], which allows to relate different classes of problems by means of a suitable transformation technique.

Let us consider (11) in the case where is the Hilbert space , , and . One gets the following representation: giving the semigroup in terms of a Gaussian expectation, with respect to the time variable , of the Schrödinger group . By applying (12) to an initial datum one can write the following chain of equalities: Under the assumptions that can be extended to an entire function on the complex plane , denoted again with , such that for any one has that is a bounded function on , a rotation of the integration path in the complex -plane gives Analogously, if is a real positive constant, a change of variable in the latter integral gives Now the integral on the right hand side of (15) makes sense for a class of functions larger than those considered so far, as well as for complex values of the parameter . Let us consider the Cauchy problem with . The next theorems give a characterization of some classes of functions such that the integral is absolutely convergent and provides a representation for the solution of the Cauchy problem (16). The following theorem is a generalization of Funaki's result [16] to the case where is an arbitrary complex parameter.

Theorem 1. *Let be a real valued function satisfying the following properties: *(1)* can be extended to an entire function on the complex plane , denoted again with ;*(2)*for any , one has that , , and , with , are bounded functions on . **
Then the solution of (16) is given by (17). *

*Proof. *By the assumption that is a bounded functions on for any , one can easily verify that the integral
is absolutely convergent. By estimates
and an integration by parts formula, one has that is given by
which is equal to .

*Remark 2. *The functions satisfying the assumptions (1) and (2) of Theorem 1 include, for instance, polynomials of arbitrary degree, as well as the functions that are Fourier transform of measures, that is, of the form
where is a complex bounded variation measure such that : (i), , being the total variation measure of , (ii)for any , the functions and are bounded on .

*Remark 3. *The Funaki formula (7) for the solution of (16) in the case where can be written in the following form:
and can be obtained as a special case of (17).

The following theorem describes a particular class of functions such that the integral (17) is absolutely convergent and provides a representation for the solution of the Cauchy problem (16).

Theorem 4. * Let be a real valued function of the form
**
where is a complex bounded variation measure on satisfying the following bound:
** being the total variation measure of .**Then the solution of the Cauchy problem (16) is given by (17). *

*Proof. *Under the stated assumption on , the integral (17) assumes the following form:
by Fubini theorem, which holds because of condition (24); the latter is equal to
which is the solution of the Cauchy problem (16).

Equation (17) can be written in terms of the expectation with respect to the measure associated to a complex random variable : where has the following distribution: being a Borel subset of the complex plane and being its characteristic function. Clearly the measure is concentrated on two rays of the complex plane and .

One can easily verify that random variable has the following properties: (i), (ii), , (iii), (iv), (v).

Moreover formula (27) can be written in Funaki's notation (see (6) and (7)). Indeed let us consider two independent Brownian motions and and define the process as Then (17) can then be written in the following form:

*Remark 5. *Analogous results can be obtained also in the case where is replaced with higher powers of the Laplacian, namely, . It is sufficient to iterate -times formula (10). One obtains a formula with multiple Gaussian integrations, similar to the one proposed, for instance, in [17]. As in the Funaki approach, the probability measure of the complex random variable can also be obtained by composing three independent Brownian motions in a suitable way. In fact any even power of the Laplacian can be handled by means of the general formula (11), but a probabilistic interpretation in terms of the composition of several independent Brownian motions is not always possible. For instance, in the case where one considers , (11) gives the following result:
with
We will not further develop these formulae here, but we will only focus on the case of .

#### 3. Feynman-Kac Type Formulae

The process defined in (29) and appearing in (30) provides a probabilistic representation for the solution of (16) (under suitable analyticity assumptions on the initial datum ). On the other hand it has not independent increments, so it does not naturally give rise to generalizations of formula (30) to the case where (16) contains also a potential , that is, to a Feynman-Kac formula of the form Indeed, by applying formally the Trotter product formula, one gets where , , is a family of independent identically distributed random variables, distributed as (see (28)). In the latter line one would be tempted to interpret the random variables as the independent increments of a complex-valued stochastic process , different from , that is, to interpret the limit (34) as the cylindrical approximations of an integral of the following form: In fact such a process cannot exist, as its construction would be possible providing the weak convergence of the sequence of complex random variables as . By using the scaling properties of the random variable , this is equivalent to the weak convergence of the sequence . As are independent identically distributed complex random variables with finite covariance, then the sequence has a Gaussian limit. Consequently, the sequence cannot converge weakly, as erroneously stated in [25], and formula (35) cannot be interpreted as an integral with respect to the measure associated to a process but just as the limit (34). Indeed the present section is devoted to the proof, for a suitable class of continuous functions and initial datum , of a Feynman-Kac formula representing the solution of the Cauchy problem (8) as the limit of a sequence of finite dimensional approximations.

The implementation of formula (35) presents some technical problems, which do not appear in the proof of the classical Feynman-Kac formula (for the heat equation with potential). The first one is the definition of the integrals involved. In fact, since the random variables are complex valued, the real valued function as well as the initial datum has to admit an analytic extension to an entire function on the complex plane. We cannot require that it is bounded on ; otherwise we could consider only the trivial case. Consequently we will integrate unbounded function, and in principle the convergence of the integrals has to be checked. In fact, for a large class of potentials, the integrals are not absolutely convergent and have to be defined in a suitable way.

The second problem concerns the proof that the integral (35) represents the solution of the Cauchy problem (8). Even if the second line of (35) recalls Trotter’s product, this formula cannot be directly applied since it does not hold for general and continuous real valued function.

The problem of the proof of a Feynman-Kac type formula for equations of the form (8) has been analyzed in [22], where the case of a linear potential is handled, and in [39], but a detailed proof for a sufficiently large class of potentials is still lacking. We generalized these results to the case where is linear in the space variable and presents an explicit time dependence.

Theorem 6. *Let be a continuous function and let us denote . Let be of the form , where is a complex bounded variation measure on satisfying the following condition:
**
Then the solution of the Cauchy problem:
**
is given by
*

*Proof. *Under the stated assumptions, the finite dimensional integrals appearing in formula (38) assume the following form:
where , , is a family of independent identically distributed random variables, distributed as . Now by applying the Fubini theorem, which holds because of condition (36), the latter is equal to
By dominated convergence theorem, the limit as of the last line is equal to
which is the solution of Cauchy problem (37).

#### 4. Conclusions

In this paper we have proposed the construction of a particular probabilistic representation for the solution of the equation in terms of a Feynman-Kac type formula. The class of potentials which can be handled by requiring that the probabilistic integrals are defined in Lebesgue sense, that is, as absolutely convergent integrals, is rather restricted because of the complex nature of the process. A generalization of these results to more general potentials requires the implementation of an integration technique, in infinite dimensions, of a different type, by relaxing the absolute convergence of the integrals, as in the cases handled, for instance, in [6] concerning the functional integral representation for the solution of Schrödinger equations. This problem will be handled in a forthcoming paper.

#### Acknowledgments

Interesting discussions with Sergio Albeverio, Giuseppe Da Prato, Paolo Dai Pra, Franco Flandoli, Enrico Priola, and Luciano Tubaro are gratefully acknowledged, as well as the financial support of Fondazione Bruno Kessler, Trento, Italy.

#### References

- E. B. Dynkin,
*Theory of Markov Processes*, Dover, Mineola, NY, USA, 2006. - I. Karatzas and S. E. Shreve,
*Brownian Motion and Stochastic Calculus*, Springer, New York, NY, USA, 1991. - E. Sinestrari, “Accretive differential operators,”
*Bollettino della Unione Matematica Italiana B*, vol. 13, no. 1, pp. 19–31, 1976. View at Google Scholar - S. Albeverio, “Wiener and feynman path integrals and their applications,”
*Proceedings of Symposia in Applied Mathematics*, vol. 52, pp. 163–194, 1997. View at Google Scholar - S. Albeverio, R. Hoegh-Krohn, and S. Mazzucchi,
*Mathematical Theory of Feynman Path Integrals—An Introduction*, vol. 523 of*Lecture Notes in Mathematics*, Springer, Berlin, Germany, 2nd edition, 2008. - S. Mazzucchi,
*Mathematical Feynman Path Integrals and Applications*, World Scientific, Singapore, 2009. - G. W. Johnson and M. L. Lapidus,
*The Feynman Integral and Feynman’s Operational Calculus*, Oxford University Press, New York, NY, USA, 2000. - R. H. Cameron, “A family of integrals serving to connect the wiener and feynman integrals,”
*Journal of Mathematical Physics*, vol. 39, pp. 126–140, 1960. View at Google Scholar - V. Y. Krylov, “Some properties of the distribution corresponding to the equation $\partial u/\partial t={(-1)}^{q+1}{\partial}^{2q}u/\partial {x}^{2q}$,”
*Doklady Akademii Nauk SSSR*, vol. 132, pp. 1254–1257, 1960 (Russian), translated as*Soviet Mathematics Doklady*, vol. 1, pp. 760–763. View at Google Scholar - K. Hochberg, “A signed measure on path space related to wiener measure,”
*Annals of Probability*, vol. 6, no. 3, pp. 433–458, 1978. View at Publisher · View at Google Scholar - K. Hochberg and E. Orsingher, “The arc-sine law and its analogs for processes governed by signed and complex measures,”
*Stochastic Processes and Their Applications*, vol. 52, no. 2, pp. 273–292, 1994. View at Google Scholar · View at Scopus - K. Hochberg, “Central limit theorem for signed distributions,”
*Proceedings of the American Mathematical Society*, vol. 79, no. 2, pp. 298–302, 1980. View at Google Scholar - K. Nishioka, “Stochastic calculus for a class of evolution equations,”
*Japanese Journal of Mathematics*, vol. 11, no. 1, pp. 59–102, 1985. View at Google Scholar - L. Beghin, K. J. Hochberg, and E. Orsingher, “Conditional maximal distributions of processes related to higher-order heat-type equations,”
*Stochastic Processes and Their Applications*, vol. 85, no. 2, pp. 209–223, 2000. View at Google Scholar · View at Scopus - D. Levin and T. Lyons, “A signed measure on rough paths associated to a PDE of high order: results and conjectures,”
*Revista Matematica Iberoamericana*, vol. 25, no. 3, pp. 971–994, 2009. View at Google Scholar · View at Scopus - T. Funaki, “Probabilistic construction of the solution of some higher order parabolic differential equation,”
*Proceedings of the Japan Academy A*, vol. 55, no. 5, pp. 176–179, 1979. View at Publisher · View at Google Scholar - K. Hochberg and E. Orsingher, “Composition of stochastic processes governed by higher-order parabolic and hyperbolic equations,”
*Journal of Theoretical Probability*, vol. 9, no. 2, pp. 511–532, 1996. View at Google Scholar · View at Scopus - E. Orsingher and X. Zhao, “Iterated processes and their applications to higher order differential equations,”
*Acta Mathematica Sinica*, vol. 15, no. 2, pp. 173–180, 1999. View at Google Scholar · View at Scopus - S. Bochner,
*Harmonic Analysis and the Theory of Probability*, University of California Press, Berkeley, Calif, USA, 1955. - K. Burdzy, “Some path properties of iterated brownian motion,” in
*Seminar on Stochastic Processes*, vol. 33 of*Progress in Probability*, pp. 67–87, Birkhäuser, Boston, Mass, USA, 1993. View at Publisher · View at Google Scholar - H. Allouba, “Brownian-time processes: the PDE connection II and the corresponding feynman-kac formula,”
*Transactions of the American Mathematical Society*, vol. 354, no. 11, pp. 4627–4637, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. Madrecki and M. Rybaczuk, “New feynman-kac type formula,”
*Reports on Mathematical Physics*, vol. 32, no. 3, pp. 301–327, 1993. View at Google Scholar · View at Scopus - K. Burdzy and A. Madrecki, “An asymptotically 4-stable process,” in
*Proceedings of the the Conference in Honor of Jean-Pierre Kahane*, pp. 97–117, Journal of Fourier Analysis and Applications, Orsay, France, June 1993. - K. Burdzy and A. Madrecki, “Ito formula for an asymptotically 4-stable process,”
*Annals of Applied Probability*, vol. 6, no. 1, pp. 200–217, 1996. View at Google Scholar · View at Scopus - P. Sainty, “Construction of a complex-valued fractional brownian motion of order N,”
*Journal of Mathematical Physics*, vol. 33, no. 9, pp. 3126–3149, 1992. View at Google Scholar · View at Scopus - R. Léandre, “Stochastic analysis without probability: study of some basic tools,”
*Journal of Pseudo-Differential Operators and Applications*, vol. 1, no. 4, pp. 389–400, 2010. View at Google Scholar - R. Léandre, “A path-integral approach to the Cameron-Martin-Maruyama-Girsanov formula associated to a bilaplacian,”
*Journal of Function Spaces and Applications*, vol. 2012, Article 458738, 9 pages, 2012. View at Publisher · View at Google Scholar - T. Hida, H. H. Kuo, J. Potthoff, and L. Streit,
*White Noise*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1995. - R. Léandre, “Path integrals in noncommutative geometry,” in
*Encyclopedia of Mathematical Physics*, J. P. Francoise, G. Naber, and S. T. Tsou, Eds., pp. 8–12, Elsevier, Oxford, UK, 2006. View at Google Scholar - R. Léandre, “Theory of distribution in the sense of connes-hida and feynman path integral on a manifold,”
*Infinite Dimensional Analysis, Quantum Probability and Related Topics*, vol. 6, no. 4, pp. 505–517, 2003. View at Publisher · View at Google Scholar · View at Scopus - Yu. L. Daletsky and S. V. Fomin,
*Measures and Differential Equations in Infinitedimensional Space*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1991. - N. V. Smorodina and M. M. Faddeev, “Probabilistic representation of solutions for a class of evolution equations,”
*Journal of Mathematical Sciences*, vol. 176, no. 2, pp. 239–254, 2011. View at Publisher · View at Google Scholar · View at Scopus - L. L. Helms, “Biharmonic functions and brownian motion,”
*Journal of Applied Probability*, vol. 4, pp. 130–136, 1967. View at Google Scholar - L. R. Bragg and J. W. Dettman, “Related problems in partial differential equations,”
*Bulletin of the American Mathematical Society*, vol. 74, pp. 375–378, 1968. View at Google Scholar - R. J. Griego and R. Hersh, “Random evolutions, Markov chains, and systems of partial differential equations,”
*Proceedings of the National Academy of Sciences of the USA*, vol. 62, pp. 305–308, 1969. View at Google Scholar - R. Hersh, “Direct solution of general one-dimensional linear parabolic equation via an abstract plancherel formula,”
*Proceedings of the National Academy of Sciences of the USA*, vol. 63, pp. 648–654, 1969. View at Google Scholar - R. Hersh, “Explicit solution of a class of higher-order abstract cauchy problems,”
*Journal of Differential Equations*, vol. 8, no. 3, pp. 570–579, 1970. View at Google Scholar · View at Scopus - R. Hersh, “The method of transmutations,” in
*Partial Differential Equations and Related Topics*, vol. 446 of*Lecture Notes in Mathematics*, pp. 264–282, Springer, Berlin, Germany, 1975. View at Publisher · View at Google Scholar - B. Gaveau and P. Sainty, “A path integral formula for certain fourth-order elliptic operators,”
*Letters in Mathematical Physics*, vol. 15, no. 4, pp. 345–350, 1988. View at Publisher · View at Google Scholar · View at Scopus