Stochastic evolution equations with singular drift and gradient noise via curvature and commutation conditions
Jonas M. T\"olle

TL;DR
This paper establishes the existence and uniqueness of solutions for a class of nonlinear stochastic evolution equations with singular drifts and gradient noise on a torus, using advanced mathematical techniques under specific geometric and commutation conditions.
Contribution
It introduces a novel approach to handle singular drifts and gradient noise in stochastic evolution equations through curvature and commutation conditions, expanding the theoretical understanding of such equations.
Findings
Proved well-posedness of the stochastic evolution equations under certain conditions.
Developed a framework using resolvent, Dirichlet forms, and an approximative Itô formula.
Extended the theory to include singular p-Laplace and total variation flows.
Abstract
We prove existence and uniqueness of solutions to a nonlinear stochastic evolution equation on the -dimensional torus with singular -Laplace-type or total variation flow-type drift with general sublinear doubling nonlinearities and Gaussian gradient Stratonovich noise with divergence-free coefficients. Assuming a weak defective commutator bound and a curvature-dimension condition, the well-posedness result is obtained in a stochastic variational inequality setup by using resolvent and Dirichlet form methods and an approximative It\^{o}-formula.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Stochastic evolution equations with singular drift and gradient noise
via curvature and commutation conditions
Jonas M. Tölle
Universität Augsburg
Institut für Mathematik
86135 Augsburg
Germany
(Date: March 18, 2024)
Abstract.
We prove existence and uniqueness of solutions to a nonlinear stochastic evolution equation on the -dimensional torus with singular -Laplace-type or total variation flow-type drift with general sublinear doubling nonlinearities and Gaussian gradient Stratonovich noise with divergence-free coefficients. Assuming a weak defective commutator bound and a curvature-dimension condition, the well-posedness result is obtained in a stochastic variational inequality setup by using resolvent and Dirichlet form methods and an approximative Itô-formula.
Key words and phrases:
nonlinear Stratonovich stochastic partial differential equation; stochastic variational inequality; singular stochastic -Laplace evolution equation; multiplicative gradient Stratonovich noise; defective commutator bound; Bakry-Émery curvature-dimension condition.
2010 Mathematics Subject Classification:
Primary: 35K55, 35K92, 60H15; Secondary: 49J40, 58J65.
The author would like to thank Max von Renesse, Hiroshi Kawabi and Lisa Beck for stimulating discussions on the topic of this work. The author is indebted to the anonymous referees for their suggestions of several improvements. Continuous support by the M.O.P.S. program is gratefully acknowledged.
1. Introduction
We shall study the nonlinear stochastic partial differential equation (nonlinear SPDE) with singular drift and gradient-type multiplicative Stratonovich noise
[TABLE]
In order to eliminate boundary or curvature effects from the underlying space, we shall consider the SPDE on , . Here, , are -coefficient fields (where denotes the transpose of ) and is an -valued standard Wiener process, . In this work, we shall prove well-posedness for initial data or, more generally, for . We assume a linear growth condition on , where we also may consider the multi-valued case . Furthermore, a geometric curvature-dimension condition on the Riemannian metric and a (defective) commutator estimate for the first order operator are assumed.
The drift operator is a distorted -Laplace-type operator, when , , which reduces to the linear diffusion operator for and includes multi-valued examples as the total variation flow operator for , , which is also called -Laplace. The equation is perturbed by independent Brownian motions , , driven by family of divergence-free -vector fields , acting in gradient direction, which may be degenerate, so that the “deterministic” PDE case is covered for .
Equations of the type (1.1), have previously been studied in [6] (for the case ), in [17] (for the case of a domain with symmetries and Neumann boundary conditions), in [42] (for Dirichlet boundary conditions), and in [7] (by an approach with weak solutions in the sense of distributions). As the coefficient field of the noise term is given by an unbounded operator in space, the Itô-analogue of (1.1) is ill-posed in general. The equation is discussed in [11, Section 7] as an example for an approach via a transformation by a group of random multipliers. Note that our approach does not rely on a transformation of (1.1) to a random PDE — rather than that, we obtain the unique solutions in terms of (stochastic) variational inequalities by a multi-step approximation procedure and a perturbation argument which relies on a weak defective commutator bound formulated in terms of Dirichlet forms, see [25, 39] for this notion. The conditions are discussed in Section 2 below. In particular, we need a curvature bound related to heat kernel estimates and lower bounds for Ricci curvature in order to derive a priori estimates for the singular equation. The higher order a priori estimates, which we need for the approximation procedure in the proof of the main result, are then derived from the defective commutator bounds. In [52], some previous results and the rough idea, which this work is based on, have been proposed by the author. In the future, equations involving more general Dirichlet operators, e.g. of nonlocal type, could be considered.
As the drift term in the SPDE (1.1) is singular, it is a known issue that solutions to the SPDE do not satisfy an Itô-equation in general (on these lines, see e.g. [28]), when, for instance, there is no Sobolev embedding for the energy of the drift available (as is e.g. in [37] for ). See [41] for an approach for nonlinearities that satisfy an superlinearity condition at infinity and see [54] for the case that the drift is perturbed with a first-order transport term. In our situation, the Itô-Stratonovich correction and the special linear structure of the noise adds some regularity to the equation, however, even for initial data in and for , to the best of our knowledge, the (limit) solutions to (1.1) can merely be characterized in terms of stochastic variational inequalities (SVI), see e.g. [10]. Among others, SVI-solutions to stochastic evolution equations have also been discussed in [45, 12, 27, 26, 29]. As seen in [30], the SVI-approach is quite robust under perturbations of the convex-subpotential-type drift with respect to the Mosco-topology. In this work, we prove existence and uniqueness for initial data in , see our main result Theorem 4.1 in Section 4 below.
Let us point out that we generalize the previously known results on solutions to (1.1) in the following aspects. As we assume periodic boundary conditions, we are able to dispense with the requirement of the boundary to be assumed to be and coefficient fields to be assumed to be perpendicular to the boundary and . We do merely require that the driving vector fields , are , divergence-free and pointwise linearly independent — we do not need to assume commutation111We remark that there is another way to dispense with commutation via a transformation to a flow of diffeomorphisms, see e.g. [15]. here, as is done e.g. in the classical linear SPDE case in [18, 19]. The introduction of the deformed diffusivity is new and may possibly be applied to more general compact Riemannian manifolds than the torus in the future. Apart from the geometric deformation, we also include more general nonlinearities (in the spirit of [28, Section 7]), extending the results of [17, 30] from homogeneous nonlinearities to ones satisfying a doubling condition, see Example 2.2 below. In fact, the nonlinearities might become multi-valued so that equation (1.1) and its Itô-Stratonovich corrected form (2.1) below are actually given by the more general subpotential equation (2.2) below, which involves the relaxed convex potentials. In the sequel, all of our analysis is generally referring to equation (2.2) below rather than to the formal equation (1.1).
Possible future topics for equations of this type include ergodicity and uniqueness of invariant measures for Markov semigroups (compare also [29, 38, 8, 22, 23] for additive noise), stability under perturbations of the drift or noise (cf. [30, 16]) and regularity (see e.g. [14, 13]) or an approach via entropy solutions [20] or renormalized solutions [48]. We also refer to [21, 53, 24] for porous media type stochastic equations with gradient noise.
Notation
Denote by \mathbb{T}^{d}:=\mathbb{R}^{d}\big{/}\mathbb{Z}^{d} the standard flat torus of dimension , equipped with the -dimensional Lebesgue measure on its Borel -algebra . On , the Euclidean norm and inner product are denoted by and respectively. For , we denote the linear operator (on some vector space) simply by . Denote , . Let denote the topological dual of . Note that the embedding is dense and compact. Let be the Dirichlet form of the Laplace-Beltrami operator on , that is, and
[TABLE]
see [25, 39]. Let be the resolvent of , i.e., for , , . For convenience, we shall also introduce the alternative resolvent , where and . Clearly, for every . The resolvent , when considered both as a map from to or as a map from to , is a contraction. Denote the Yosida-approximation of * *by , . is a negative definite bounded linear operator with operator norm bounded by . Furthermore, for , , let be inner products for with the property that , are mutually equivalent for . For , define also approximate forms \mathcal{E}^{(\beta)}(u,v):=\beta(u-\text{\beta}G_{\text{\beta}}u,v)_{H}, , see e.g. [39, Chapter I, p. 20]. Set also , , . Denote by the *heat semigroup *associated to . We shall denote the trace of a square matrix by .
Organization of the paper
In Section 2, which is subdivided into several subsections, we shall discuss the different terms of equation (1.1) and the main assumptions on those. Here, also the main hypotheses of this work are stated and discussed. In Section 3, we shall introduce and discuss the concept of so-called stochastic variational inequality (SVI) solutions to our equation. Section 4, which is subdivided into several subsections, contains the statement and the proof of our main result, including the a priori estimates, the passage to the limit of the approximating equations, the existence and the uniqueness result. In Appendix A, we recall the conditions from [28] which are needed to guarantee the existence of approximating solutions to equation (1.1). The concluding Appendix B contains the postponed proof of Proposition 2.11.
2. Main hypotheses
Let be a Wiener process on , modeled on a filtered probability space that satisfies the usual conditions.
Consider the following Itô-SPDE in ,
[TABLE]
Equation (2.1) is the formal analogue to equation (1.1), after passing over to an (merely formal) Itô-Stratonovich correction, see e.g. [34, 33] and compare with the previous works [7, 17, 6, 52].
The precise assumptions on and will be stated below. The operator is a realization of , which equals the Itô-Stratonovich correction of the noise term , whenever for all .
2.1. The nonlinearity
To improve precision once more, instead of (2.1), we shall in fact study the following Itô-SPDE (stochastic evolution equation / inclusion222The set membership symbol is to be understood as an equality whenever the r.h.s. is single-valued.)
[TABLE]
where denotes the subdifferential333The subdifferential of a convex, lower semi-continuous, proper map is defined by , , whenever for every , see e.g. [9]. of , which, in turn, is defined to be the lower semi-continuous (l.s.c.) envelope of the map
[TABLE]
where is some convex potential of , that is, for every and all . The l.s.c. envelope is defined by
[TABLE]
see e.g. [3] for further details.
Assumption 2.1**.**
Let us assume that there exist constants and a map such that
- (N)
, for all and is convex, continuous, and satisfies , , , for .
Note that the doubling condition , is also known as -condition in the literature, see e.g. [44]. Our assumptions on shall be made precise further below.
Example 2.2**.**
Condition (N) is e.g. satisfied for the following choices of . For , let
**-Laplace, total variation flow: **
, ,
**Logarithmic diffusion: **
,
**Minimal surface flow: **
,
**Curve shortening flow: **
.
Remark 2.3*.*
Note that e.g. for , the subdifferential becomes multi-valued, and the l.s.c. envelope is a Radon measure on the space of functions of bounded variation, see [3] for details. At this point it is enough to recall that the above l.s.c. envelope exists in .
We continue with two technical lemmas needed later.
Lemma 2.4**.**
Assume that (N) holds. Then
[TABLE]
Proof.
By (N) and [51, Chapter II, Example 8.A], there exists a function such that is non-decreasing, right-continuous and satisfies
[TABLE]
Also, for , . Using (N) again, for ,
[TABLE]
However, it is easy to see that from , it follows that , for all and , cf. [51, Chapter II, Proposition 8.6]. Hence the claim follows. ∎
For the statement of the following lemma and for use further below, we shall introduce the notion of the so-called Moreau-Yosida approximation of , that is, the family of continuous convex functions defined by the following variational formula
[TABLE]
see [2, p. 266] or [9, p. 97] for further details.
Lemma 2.5**.**
Assume that (N) holds and let , be the Moreau-Yosida approximation of . Then, there exists a constant (not depending on ) such that
[TABLE]
Proof.
As above, denote . By the arguments preceding [30, Eq. (A.4) in Appendix A], we get in this slightly more general situation that
[TABLE]
By [28, Proof of Proposition 7.1] and (N), there exists a constant , such that
[TABLE]
Again, by (N) and by Lemma 2.4 there exists another constant such that
[TABLE]
Combining these inequalities finishes the proof. ∎
2.2. Basic hypotheses on the coefficient fields and
We shall formulate the main assumptions on and .
Assumption 2.6**.**
Assume that and that there exists a constant such that
[TABLE]
Assumption 2.7**.**
Assume that such that the rows , of satisfy
[TABLE]
Let . Let denote the Dirichlet operator associated to the Dirichlet form
[TABLE]
where denotes the application of matrix-multiplication for . For smooth functions , by -regularity of the coefficients, we have that , where denotes the matrix-adjoint of . Denote the associated Dirichlet operators by , , respectively, with Dirichlet forms , , respectively.
Remark 2.8*.*
Note that Assumption 2.6 implies that and that is uniformly elliptic with domain equal to . Note also that and thus may be degenerate.
Lemma 2.9**.**
Suppose that . Then condition (E) is equivalent to being linearly independent in for each .
Proof.
First note that is precisely the Gram matrix of the vectors . Consider the statement of being pointwise linearly independent. It is well-known that this is equivalent to the Gram matrix being positive definite, see e.g. [35, Chapter 10], which in turn is equivalent to (E) with , being the smallest eigenvalue of (all eigenvalues are strictly positive and real). By continuous dependence of the eigenvalues of on the space variable , noting that is compact, we argue by contradiction to see that . ∎
2.3. Curvature-dimension condition
In all of this subsection, assume condition (E) and (N). Set and set . By condition (E), is a Riemannian manifold which is quasi-isometric444See [32, p. 93] for this notion. to , when equipped with the flat metric. We write if we want to emphasize the choice of the metric. We denote the volume measure on by . Note that , equipped with the Lebesgue measure, is thus a weighted manifold with density with respect to . Note that by [32, Exercise 3.12], on the weighted manifold , where denotes the weighted Laplace-Beltrami operator of .
Definition 2.10**.**
Let . We say that satisfies a *Bakry-Émery curvature-dimension condition * if there exists with
[TABLE]
Proposition 2.11**.**
Suppose that (E) holds. Suppose that and that for all ,
[TABLE]
where . Then condition (BE) holds for some .
Proof.
See Appendix B. ∎
Note that condition (2.3) is symmetric in the sense that interchanging the indices and yields the same condition.
Theorem 2.12**.**
Let be the heat semigroup associated to . Then the following condition is equivalent to condition .
There exists a constant , such that
[TABLE]
Proof.
See [1, 5, 56, 57, 31, 58, 55] for the terminology and further results on equivalent curvature-dimension conditions as well as Ricci curvature bounds in weighted Riemannian manifolds.
Lemma 2.13**.**
Let condition (N) and (E) hold. Let be the heat semigroup associated to . Assume that holds for some . Let , be the resolvent associated to . Then
[TABLE]
for any and any , where is as in (2.3).
Proof.
Compare with [28, Proof of Example 7.11]. Let be the -semigroup associated to , . Then by the Trotter product formula [46, Ch. VIII.8], , . We get for and that
[TABLE]
where we have used Theorem 2.12. Since , is Markovian symmetric on , we get by an application of [40, Theorem 3] and condition (N) that for ,
[TABLE]
Altogether,
[TABLE]
density of and Lebesgue’s dominated convergence theorem completes the proof. ∎
Remark 2.14*.*
Note that
- (i)
the above statement remains true if holds for . If one investigates the proof carefully, one sees that then it even holds that
[TABLE]
for any and any ; 2. (ii)
the above statement remains true if