Stochastic Cahn-Hilliard equation with double singular nonlinearities and two reflections
Arnaud Debussche (IRMAR), Ludovic Gouden\`ege (IRMAR)

TL;DR
This paper studies a complex stochastic PDE with double singular nonlinearities and reflections, establishing existence, uniqueness, ergodicity, and exponential mixing of solutions under challenging conditions.
Contribution
It introduces a novel approach to handle the bi-Laplacian's lack of maximum principle, proving well-posedness and ergodic properties for the stochastic Cahn-Hilliard equation with singular nonlinearities.
Findings
Existence and uniqueness of solutions for initial data in (-1,1)
Identification of the unique invariant measure as ergodic
Proof of exponential mixing for the system
Abstract
We consider a stochastic partial differential equation with two logarithmic nonlinearities, with two reflections at 1 and -1 and with a constraint of conservation of the space average. The equation, driven by the derivative in space of a space-time white noise, contains a bi-Laplacian in the drift. The lack of the maximum principle for the bi-Laplacian generates difficulties for the classical penalization method, which uses a crucial monotonicity property. Being inspired by the works of Debussche, Gouden\`ege and Zambotti, we obtain existence and uniqueness of solution for initial conditions in the interval . Finally, we prove that the unique invariant measure is ergodic, and we give a result of exponential mixing.
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 Cahn-Hilliard equation with double singular nonlinearities and two reflections
Arnaud Debussche, Ludovic Goudenège
ENS Cachan - Antenne de Bretagne, Campus de Ker-Lann, 35170 Bruz, France
[email protected], [email protected]
footnotetext: AMS 2000 subject classifications. 60H15, 60H07, 37L40.
Key words and phrases : Cahn-Hilliard, stochastic partial differential equations, integration by parts formulae, reflection measures, invariant measures, singular nonlinearity, double, two.
Résumé
On considère une équation aux dérivées partielles stochastique possédant deux non-linéarités de type logarithmique, avec deux réflexions en et sous la contrainte de conservation de masse. L’équation, dirigée par un bruit blanc en espace et en temps, contient un double Laplacien. L’absence de principe de maximum pour le double Laplacien pose des difficultés pour l’utilisation d’une méthode classique de pénalisation, pour laquelle une importante propriété de monotonie est utilisée. Etant inspiré par les travaux de Debussche, Goudenège et Zambotti, on démontre l’existence et l’unicité de solutions pour des données initiales entre et . Enfin, on démontre que l’unique mesure invariante est ergodique, et on énonce un résultat de mélange exponentiel.
Abstract
We consider a stochastic partial differential equation with two logarithmic nonlinearities, with two reflections at and and with a constraint of conservation of the space average. The equation, driven by the derivative in space of a space-time white noise, contains a bi-Laplacian in the drift. The lack of the maximum principle for the bi-Laplacian generates difficulties for the classical penalization method, which uses a crucial monotonicity property. Being inspired by the works of Debussche, Goudenège and Zambotti, we obtain existence and uniqueness of solution for initial conditions in the interval . Finally, we prove that the unique invariant measure is ergodic, and we give a result of exponential mixing.
Introduction and main results
The Cahn-Hilliard-Cook equation is a model to describe phase separation in a binary alloy (see [6], [7] and [8]) in the presence of thermal fluctuations (see [11] and [27]). It takes the form:
[TABLE]
where denotes the time variable and is the Laplace operator. Also represents the ratio between the two species and the noise term accounts for the thermal fluctuations. The nonlinear term has the double-logarithmic form:
[TABLE]
where and are temperatures with .
The study of this equation presents several difficulties. First, the singularities at have to be treated carefully. Also, since it is a fourth order equation, no comparison principle holds.
The deterministic equation where is replaced by a polynomial function have first been studied (see [7], [27] and [32]). Then non smooth have been considered (see [5] and [17]).
Phase separation have been analysed thanks to this model: see for example the survey [31], and the references therein, or others recent results on spinodal decomposition and nucleation in [1, 4, 24, 28, 29, 35, 36, 37].
In the case of a polynomial nonlinearity, some results have been obtained in the stochastic case (see [2, 3, 9, 10, 12, 19]).
Note that the solutions of the equation with polynomial nonlinearity do not remain in in general, and their physical interpretation is not clear.
To our knowledge, the case of the logarithmic nonlinearity in the presence of noise have never been studied. The presence of noise has a strong effect and equation (0.1) cannot have a solution. Indeed, a solution should remain in which is impossible with an additive noise. Two reflection measures have to be added to the model to remedy this problem. In this article, we propose to study:
[TABLE]
where the measures are subject to the contact conditions almost surely:
[TABLE]
The stochastic heat equation with reflection, i.e. when the fourth order operator is replaced by the Laplace operator, is a model for the evolution of random interfaces near a hard wall. It has been extensively studied in the literature (see [16], [21], [22], [33] [38], [39] and [40]). Essential tools in these articles are the comparison principle and the fact that the underlying Dirichlet form is symmetric so that the invariant measure is known explicitely.
In our case, we consider a noise which is obtained as the space derivative of the space-time white noise. In other words, the noise is the time derivative of a cylindrical Wiener process in . This is physically reasonable since the Cahn-Hilliard equation can be interpreted as a gradient system in this space. With such noise, the system is still symmetric and the invariant measure is known explicitely. As in the second order case, we use this fact in an essential way.
However, as already mentioned, no comparison principle holds and new techniques have to be developed. The equation (0.3) has been studied with a single reflection and when no nonlinear term is taken into account in [18]. The reflection is introduced to enforce positivity of the solution. Various techniques have been introduced to overcome this lack of comparison principle. Moreover, as in the second order case, an integration by part formula for the invariant measure has been derived. Then, in [23], a singular nonlinearity of the form or have been considered. Existence and uniqueness of solutions have been obtained and using the integration by parts formula as in [39], it has been proved that the reflection measure vanishes if and only if . In particular, for a logarithmic nonlinearity, the reflection is active.
Here, we consider the original Cahn-Hilliard-Cook model (0.1) with the double-logarithmic nonlinear term (0.2). The noise is as in the above mentioned articles and we still have an explicit invariant measure. Our method mixes ideas from [18], [23] and [39]. Additional difficulties are overcome, the main one being to understand how to deal with the nonlinear term. Indeed, in [23], the positivity of the nonlinear term was essential. We overcome this difficulty thanks to a delicate a priori estimate. Our main results state that equations (0.3), (0.4) together with an initial condition have a unique solution (see Proposition 2.3 and Theorem 1.1). As in [18], it is constructed thanks to the gradient structure of (0.3) and strong Feller property. Moreover, we prove that this solution is the limit of the solution of the Cahn-Hilliard-Cook equation with polynomial nonlinearity without reflections. This justifies the use of the polynomial models. We also prove that the invariant measure is unique and ergodic. Such property is very easy to obtain if is small (see [18]) or in the polynomial case (see [12]). Finally, a stronger result of exponential mixing is given in the last Theorem 3.1. It is based on coupling and arguments developped by Odasso in [34].
In future studies, we shall generalize the integration by part formula obtain in [18] to prove that the reflection measure does not vanish. The presence of two reflection measures introduces additional difficulties. In the second order case, this has been studied in [20].
1 Preliminaries
We denote by the scalar product in ; is the realization in of the Laplace operator with Neumann boundary condition, i.e.:
[TABLE]
where we use and to denote the Sobolev space and its associated norm. Remark that is self-adjoint on and we have a complete orthonormal system of eigenvectors in for the eigenvalues . We denote by the mean of :
[TABLE]
We remark that is invertible on the space of functions with average. In general, we define .
For , we define by classical interpolation. We set . It is endowed with the classical seminorm and norm :
[TABLE]
is associated to the scalar product . To lighten notations, we set and . The average can be defined in any by . It plays an important role and we often work with functions with a fixed average . We define for all .
We use the following regularization operators:
[TABLE]
It is defined on and can extended to any . Clearly converges to in if . Moreover, it is well known that if , then the converges holds in . Note also that is self-adjoint in and commutes with .
The covariance operator of the noise is the operator defined by
[TABLE]
Note that
[TABLE]
We denote by the space of all Borel bounded functions on . We set for with and , and for . Given a measure on and a continuous function on , we write
[TABLE]
For , we define:
[TABLE]
and the following antiderivative of :
[TABLE]
With these notations, we rewrite (0.3) in the abstract form:
[TABLE]
where is a cylindrical Wiener process on .
Definition 1.1
Let . We say that , defined on a filtered complete probability space , is a weak solution to (0.3) on for the initial condition if:
- (a)
a.s. and , 2. (b)
a.s. are two positive measures on , such that for all , 3. (c)
* is a cylindrical Wiener process on ,* 4. (d)
the process is -adapted, 5. (e)
a.s. , 6. (f)
for all and for all :
[TABLE] 7. (g)
*a.s. the contact properties hold : *
* and , that is,*
[TABLE]
Finally, a weak solution is a strong solution if the process is adapted to the filtration
Remark 1.1
In (f), the only term where we use the function is well defined. Indeed, by (e) we have and by Sobolev embedding . Hence the notation should be interpreted as a duality between and .
The solution of the linear equation with initial data is given by
[TABLE]
As easily seen this process is in (see [14]). In particular, the mean of is constant and the law of the process is the Gaussian measure:
[TABLE]
If we let , the law of converges to the Gaussian measure on :
[TABLE]
Notice that is concentrated on .
In order to solve equation (1.2), we use polynomial approximations of this equation. We denote by the sequence of polynomial functions which converges to the function on , defined for by:
[TABLE]
We use the following antiderivative of defined by:
[TABLE]
Then for , we study for the following polynomial approximation of (1.2) with an initial condition :
[TABLE]
This equation has been studied in [12] in the case . The results generalize immediately and it can be proved that for any , there exists a unique solution in . It is a solution in the mild or weak sense. Moreover the average of does not depend on .
For each , (1.2) defines a transition semigroup :
[TABLE]
Existence of an invariant measure can be proved as in [12].
Using Galerkin approximation and Bismut-Elworthy-Li formula, it can be seen that is Strong Feller. More precisely, for all , and :
[TABLE]
Irreducibility follows from a control argument. By Doob Theorem we deduce that there exists an unique and ergodic invariant measure .
It is classical that equation (1.3) is a gradient system in and can be rewritten as:
[TABLE]
where denotes the gradient in the Hilbert space , and:
[TABLE]
The measure is therefore given by:
[TABLE]
where is a normalization constant.
We prove in section 2 that, for , the sequence converges to the measure
[TABLE]
where
[TABLE]
and
[TABLE]
In section 2, we prove the following result.
Theorem 1.1
Let . Let such that , then there exists a continuous process denoted and two nonnegative measures and such that is the unique strong solution of (0.3) with a.s.
The Markov process is continous and has for transition semigroup which is strong Feller on .
For all and , is the limit in distribution of .
Finally is an invariant measure for .
In all the article, C denotes a constant which may depend on and its value may change from one line to another.
2 Proof of Theorem 1.1
2.1 Pathwise uniqueness
We first prove that for any pair , of weak solutions of (0.3) defined on the same probability space with the same driving noise and with , we have . This pathwise uniqueness will be used in the next subsection to construct stationary strong solutions of (0.3).
Proposition 2.1
Let . Let be two weak solutions of (0.3) with . Then .
Proof** : **We use the following Lemma from [23].
Lemma 2.1
Let be a finite measure on and . Suppose that there exists a positive continuous function such that :
- i)
for all , for all , such that , , 2. ii)
for all , with ,
then is the null measure.
Let , and , is the solution of the following equation:
[TABLE]
Taking the scalar product in with and integrating in time, we obtain since has zero average:
[TABLE]
For all ,
[TABLE]
where and are the classical norms on the space . The latter term converges to zero since converges uniformly to on . Since is nonincreasing,
[TABLE]
Taking the limit in (2.2) as grows to infinity, we obtain:
[TABLE]
We now write
[TABLE]
by the contact condition and the positivity of the measures. It follows:
[TABLE]
By Gronwall Lemma, and letting , we have for all . Since , we deduce for all . Moreover, with the definition of a weak solution, we see that :
[TABLE]
By density, we obtain that and satisfy the hypothesis of Lemma 2.1, and therefore is the null measure. And since and have disjoint supports, then and are the null measure, i.e. and .
2.2 Convergence of invariants measures
We know (see [18]) that is the law of , where is brownian motion. Then for , we remark the following inclusion :
[TABLE]
and we have a similar result for . Therefore with . Let us define the potential associated to the function :
[TABLE]
We have the following result :
Proposition 2.2
For ,
[TABLE]
where is a normalization constant.
Proof** : **Let . We want to prove that
[TABLE]
We first prove,
[TABLE]
Since , we can restrict our attention to . Then if there exists such that or , being the Lebesgue measure. Suppose , then we have since is positive and non increasing on
[TABLE]
And this latter term converges to zero as grows to infinity.
Now for , converges to almost everywhere as grows to infinity. Moreover , and by the dominated convergence Theorem, we deduce (2.4). Finally, (2.3) follows again by dominated convergence Theorem.
2.3 Existence of stationary solutions
In this section, we prove the existence of stationary solutions of equation (1.2) and that they are limits of stationary solutions of (1.3), in some suitable sense. Fix and consider the unique (in law) stationary solution of (1.3) denote in . We are going to prove that the laws of weakly converge as grows to infinity to a stationary strong solution of (0.3).
Proposition 2.3
Let and , converges in probability as grows to infinity to a process in . Moreover almost surely, and setting
[TABLE]
and
[TABLE]
then converges in probability to such that is a stationary strong solution of (0.3).
Proof** : **Proceeding exactly as in [18] (see Lemma 5.2), we prove that the laws of are tight in , . We have set , . We therefore can extract convergent subsequences. Let be such a subsequence. Using Skohorod theorem, one may find a probability space and a sequence of random variables on this probability space with the same laws as which converge almost surely.
Below, we show in Step 1 that its limit satisfies almost surely. Then in Step 2, we prove that the measures , defined as above with replaced by , converges to two positive measures and that is a weak solution in the probabilistic sense. It then remains to use pathwise uniqueness to conclude in Step 3. In this proof, we only treat the case . This assumption is not essential at all but lightens the computations. For , an extra term has to be taken into account. It is very easy to deal with.
**Step [.**1]
Applying Ito formula to , we obtain
[TABLE]
Note that
[TABLE]
We set
[TABLE]
and deduce
[TABLE]
Thus, for all , the laws of are tight. Therefore the laws of are tight and using Skohorod theorem on this sequence, we can assume that , and, for , converge almost surely. We have defined as above with instead of . In particular, is bounded almost surely:
[TABLE]
where is random. The first three terms are clearly also bounded almost surely. This uses the fact that is a bounded operator from to . Since has finite dimensional range, we obtain
[TABLE]
for a different random constant .
Let us choose and take such that
[TABLE]
and such that for
[TABLE]
Then, for all ,
[TABLE]
Moreover, if then and
[TABLE]
Similarly, if then and
[TABLE]
Finally, noticing that is uniformly bounded by a constant on , we deduce
[TABLE]
Thanks to (2.5), we obtain
[TABLE]
where the value of the random constant has again changed. It easily deduced from this uniform bound that almost everywhere with respect to and and by Fatou Lemma that almost surely.
1]
**Step [.**2]
Let now be the following measure on :
[TABLE]
and and the positive and negative parts:
[TABLE]
By step 1, and we can define the following measure:
[TABLE]
and the positive and negative parts:
[TABLE]
By (2.6), is bounded in . We deduce that has a subsequence which converges to a measure . Note that this subsequence may depend on the random parameter . We set .
Thanks to Fatou Lemma we have the following inequality for all nonnegative:
[TABLE]
Therefore converges to a positive measure. Similarly, converges to a positive measure. It follows:
[TABLE]
where and are the positive and negative parts of .
Let us now show that the contact conditions holds for and . Let us define the following measures for and .
[TABLE]
Clearly converges to , it follows
[TABLE]
Since converges uniformly to zero, we deduce:
[TABLE]
Letting , we obtain the first contact condition since the left hand side clearly converges to . The second is obtained similarly.
We now prove that does not have more than one limit point so that in fact the whole sequence converge to . Let , be two limit points.
For all and for all :
[TABLE]
We deduce
[TABLE]
And by density
[TABLE]
for any such that . Since by the contact condition
[TABLE]
We deduce from Lemma 2.1 that .
2]
**Step [.**3]
We use a result form [25] that allows to get the convergence of the approximated solutions in probability in any space in which these approximated solutions are tight.
Lemma 2.2
Let be a sequence of random elements on a Polish space endowed by its borel -algebra. Then converges in probability to an -valued random element if and any if from every pair of subsequences , one can extract a subsequence which converges weakly to a random element supported on the diagonal .
Assume and are two arbitrary subsequences. Clearly, the process is tight in a suitable space. By Skorohod’s theorem, we can find a probability space and a sequence of processes such that almost surely and has the same distribution as for all . In the Skorohod’s space, the approximated measures respectively converge to two contact measures and . By the second step, and are both weak solutions of (0.3). By uniqueness, necessarily and . Therefore the subsequence converges in distribution to a process supported on the diagonal. We use Lemma 2.2 to prove that the sequence converges in probability to a process . Clearly is stationary. Reproducing the argument of Step 1 and Step 2, we prove that it is a strong solution of (0.3) and the convergence of the contact measures.
3]
2.4 Convergence of the semigroup
First we state the following result which is a corollary of Proposition 2.3.
Corollary 2.1
Let .
- i)
There exists a continuous process with and a set dense in , such that for all there exists a unique strong solution of equation (0.3) given by . 2. ii)
The law of is a regular conditional distribution of the law of given .
Proof** : **By Proposition 2.3, we have a stationary strong solution in , such that and are independent. Conditioning on the value of , with , we obtain for -almost every a strong solution that we denote for all and for all . This process is the desired process. Indeed, since the support of is , we have a strong solution for a dense set in .
Notice that all processes with are driven by the same noise and are continuous with values in . Moreover, we have the following obvious identity:
[TABLE]
and by density we obtain a continuous process in for all .
Proposition 2.4
Let , for all and :
[TABLE]
Moreover the Markov process is strong Feller and its transition semigroup is such that:
[TABLE]
Proof** : **Since is tight in , then there exists an increasive sequence of compact sets in such that:
[TABLE]
Set . Since the support of is in and , then is dense in . Fix , by (1.4), for any :
[TABLE]
Let be any sequence in . With a diagonal procedure, by Arzelà-Ascoli Theorem, there exists a subsequence and a function such that:
[TABLE]
By density, can be extended uniquely to a bounded Lipschitz function on such that
[TABLE]
Note that the subsequence depends on . Therefore, we have to prove that the limit defines a semigroup and does not depend on the chosen subsequence.
By Proposition 2.3, we have for all :
[TABLE]
Thus, by Corollary 2.1, we have the following equality:
[TABLE]
Since and are continuous on , and , the equality (2.9) is true for all . Moreover the limit does not depend on the chosen subsequence, and we obtain (2.7). Letting in (1.4), we deduce (2.8).
2.5 End of the proof of Theorem 1.1
We have proved that there exists a continous process which is a strong solution of equation (0.3) for an in a dense space. In this section, we prove existence for an initial condition in with .
By Corollary 2.1 we have a process , such that for all in a set dense in we have a strong solution of (0.3) with initial condition . By Proposition 2.3, the Markov process has transition semigroup on .
The strong Feller property of implies that for all and the law of is absolutely continous with respect to the invariant measure . Indeed, if , then . So for -almost every and by continuity for all .
Therefore almost surely for all and . Fix , denote for all :
[TABLE]
and the measures such that for all , and for all :
[TABLE]
So we have a process and two measures and on which is finite on for all