Skip to main content Skip to navigation

Sample path large deviations and concentration

Minimisers of the Rate Function
We are considering minimisers of the rate function associated to the LDP for a Laplacian particle interaction and a Gaussian potential.

We consider minimising the rate function

 \Sigma^\Delta(h) = \frac{1}{2}\int_0^1 h''(t)^2 \text{d}t-\tau|\{t\in[0,1]:h(t)=0 \}|,

over the sets

 H^2_\bm{a} :=\{h\in H^2([0,1]):h(0)=a,h'(0)=\alpha\},
 H^2_\bm{a,b}:=\{h\in H^2([0,1]):h(0)=a,h'(0)=\alpha,h(1)=b,h'(1)=\beta\}.

These two cases are called the free boundary problem and Dirichlet boundary problem, respectively.  \tau is called the free energy for the pinning, and is treated as a constant throughout.

No wetting, Pinning

In the free boundary case, the set of minimisers is contained in the set \{\overline{h},\hat{h},\check{h} \}, where \overline{h},\hat{h},\check{h} correspond to the green, orange and blue curves in the figure below, respectively.

We note that these minimisers only exist at the same time for some particular choices of  \bm{a},\tau . Using Euler-Lagrange, we can see the minimisers will be cubics. We overview the different cases and give a description of the cubics.

Firstly, we split our interval up in to two. Namely  [0, t_1] , where we take  h^* to be some cubic, and  [t_1,1 ] , where  h^* is zero. For \tau= 0, we have the unique minimiser \overline{h }.

Secondly, by noting  \Sigma^\Delta(\bar{h} )=0 , we see that any minimiser must satisfy  \Sigma^\Delta(h^*)\leq 0 . A simple calculation reveals that

 \Sigma^\Delta(h^* ) = \frac{2(3a^2 + 3a\alpha t_1 + \alpha^2t_1^2)}{t_1^3} - (1-t_1)\tau.

We can then see suitable minima occur at the points

 t_1 = \frac{\pm \alpha + \sqrt{ \alpha^2 \pm 6\sqrt{2\tau}a}}{\sqrt{2\tau}}\, \text{ or } t_1 = \frac{\pm \alpha - \sqrt{\alpha^2 \pm 6\sqrt{2\tau}a}}{\sqrt{2\tau}}\,.

After substituting these values into  \Sigma^\Delta(h^* ) and equating to zero, we reveal critical cases corresponding to multiple minimisers, in terms of  \bm{a} and  \alpha only. We may then plot this curve to see regions where one, two and three minimisers occur. This plot is shown below. Inside the sausage the region the minimiser succumbs to the pinning strength. Most interestingly, we will get multiple minimisers when parameters are such that they lie on the boundary of this region. The most peculiar situation of three minimisers happen at the two points where blue and yellow lines cross.

The same analysis can be conducted for the Dirichlet boundary case, albeit more messy. This gives us 5 potential minimisers and in principle, the regions in parameter space where minimisers change or multiple minimisers occur can be calculated.

Wetting, No pinning

For the case of wetting, we consider minimising the rate function

\frac{1}{2}\int_0^1 h''(t)^2 \text{d}t

over the convex spaces

 H^{2,+}_{\bm{a}} := \{ h\in H^2_\bm{a} ([0,1]) : h\geq0 \},
 H^{2,+}_\bm{a,b}:=\{h\in H^2_\bm{a,b}([0,1]) : h\geq0\}

for the free and Dirichlet boundaries respectively.

Note: Here we have taken \tau= 0 so we have no pinning.

This problem is in fact equivalent to the biharmonic obstacle problem

 h\in H^{2,+}_\bm{a}:\int_0^1 h''(v'' -h'')\geq0 \forall v\in H^{2,+}_\bm{a},

where the function space depends whether we work with a free boundary or Dirichlet boundary.

It is well known this problem has a unique solution. Furthermore, the obstacle problem has an equivalent strong formulation known as the free boundary problem, given by

 h \in H^2_\bm{a}([0,1]):h\geq 0 \text{a.e. in} [0,1], h^{(iv)}\geq0 \text{a.e. in} [0,1], h^{(iv)}h=0 \text{a.e. in} [0,1].

where the fourth derivative is to be understood in the sense of distributions.

In our case, we note that on \Omega^+:= \{t \in [0,1] \, | \, h(t) > 0\}, u is a biharmonic function, which in 1D simply means that it's a cubic outside of the coincidence set \Omega_c = \{t \in [0,1] \, | \, h(t) = 0 \}.
In order to explicitly find a minimiser over H^{2,+}_{\bm{a} } we proceed as follows:

  1. Find conditions on boundary parameters for which \overline{h } (introduced in the previous section) is naturally positive on [0, 1] - call this set M_{+}. If \bm{a}\in M_+ then there is nothing to do and \overline{h } is the minimiser.
  2. If \bm{a} \not\in M_+, then assume there is some coincidence set [t_1,1] for t_1\in[0,1] on which the minimiser is a constant zero function. Let us denote such function by h_* and note the following:
    • on [0,t_1] we have a cubic, which necessarily satisfy h(t_1)= h'(t_1)= 0 (as h has to be differentiable to be in the set);
    • a potential minimiser would ideally like to coincide with the obstacle at t_1= 1, as then it has "more time" to counter its willingness to go negative: a simple calculation reveals that
       \Sigma^{\Delta}(h_* ) = \frac{12a^2}{t_1^3} + \frac{12a\alpha}{t_1^2} + \frac{4\alpha^2}{t_1},

      which, if treated as a function of  t_1 is non-increasing on  [0,1].
    • in many cases for  t_1 is large enough the aforementioned cubic will itself go negative before reaching  t_1, so one has to be careful.
  3. in many cases for t_1 is large enough the aforementioned cubic will itself go negative before reaching t_1 , so one has to be careful.

Using this recipe, explicit minimisers that occur in the case of wetting can be found. For the Dirchlet boundary we proceed similarly adding a mirror argument for the right-hand boundary, where we now assume that the coincidence set is [t_1,1-t_2].