Hostname: page-component-cd9895bd7-q99xh Total loading time: 0 Render date: 2024-12-24T00:40:19.621Z Has data issue: false hasContentIssue false

On the Γ-convergence of the Allen–Cahn functional with boundary conditions

Published online by Cambridge University Press:  12 February 2024

Dimitrios Gazoulis*
Affiliation:
Department of Mathematics and Applied Mathematics, University of Crete, Heraklion 70013, Greece Institute of Applied and Computational Mathematics, FORTH, Heraklion 70013, Crete, Greece ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

We study minimizers of the Allen–Cahn system. We consider the $\varepsilon$-energy functional with Dirichlet values and we establish the $\Gamma$-limit. The minimizers of the limiting functional are closely related to minimizing partitions of the domain. Finally, utilizing that the triod and the straight line are the only minimal cones in the plane together with regularity results for minimal curves, we determine the precise structure of the minimizers of the limiting functional, and thus the limit of minimizers of the $\varepsilon$-energy functional as $\varepsilon \rightarrow 0$.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - ND
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by-nc-nd/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction

In this work we are concerned with the study of vector minimizers of the Allen–Cahn $\varepsilon$-functional,

(1.1)\begin{align} & J_{\varepsilon} (u, \Omega) := \int_{\Omega} \left( \frac{\varepsilon}{2} |\nabla u|^2 + \frac{1}{\varepsilon} W(u) \right) \,{\rm d}x , \nonumber\\ & u : \Omega \rightarrow \mathbb{R}^m , \end{align}

where $\Omega \subset {\mathbb {R}}^n$ is an open set and $W$ is a $N$-well potential with $N$ global minima.

Let

(1.2)\begin{equation} u_{\varepsilon} := \mathop{\rm argmin}\limits_{v \in W^{1,2}(\Omega ; {\mathbb{R}}^m)} \lbrace J_{\varepsilon}(v, \Omega) : v|_{\partial \Omega} = g_{\varepsilon}|_{\partial \Omega} \rbrace\ ,\ \textrm{where} \ g_{\varepsilon} \in W^{1,2} (\Omega ; {\mathbb{R}}^m). \end{equation}

Thus, $u_{\varepsilon } \in W^{1,2}(\Omega ; {\mathbb {R}}^m)$ is a weak solution of the system

(1.3)\begin{equation} \begin{cases} \varepsilon \Delta u_{\varepsilon} - \frac{1}{\varepsilon} W_u (u_{\varepsilon}) =0 \ ,\ \textrm{in} \ \Omega, \\ u_\varepsilon = g_\varepsilon,\ \textrm{on} \ \partial \Omega. \end{cases} \end{equation}

We study the asymptotic behaviour of $u_\varepsilon$ within the framework of $\Gamma$-convergence. Moreover, we analyse the relationship between minimizers of the Allen–Cahn system and minimizing partitions subject to Dirichlet boundary conditions. For some particular assumptions on the limiting boundary conditions, we will prove uniqueness for the limiting geometric problem and we will determine the structure of the minimizers of the limiting functional.

1.1 Main results

Hypothesis on $W$:

(H1) $W \in C^{1,\alpha }_{loc}({\mathbb {R}}^m ; [0,\,+ \infty )) \ ,\,\ \lbrace W=0 \rbrace = \lbrace a_1,\,a_2,\,...,\,a_N \rbrace \ ,\,\ N \in \mathbb {N} \ ,\, a_i$ are the global minima of $W$. Assume also that

\[ W_u(u) \cdot u > 0 \quad \textrm{and} \quad W(u)\geq c_1 |u|^2 \ ,\ \textrm{if}\ |u |>M . \]

Hypothesis on the Dirichlet data:

(H2)(i) $|g_{\varepsilon }| \leq M \ ,\, g_\varepsilon \stackrel {L^1(\Omega )}{\longrightarrow } g_0$ and $J_\varepsilon (g_\varepsilon,\, \Omega _{\rho _0} \setminus \Omega ) \leq C \ ,\,$ where $\partial \Omega$ is Lipschitz and $\Omega _{\rho _0}$ is a small dilation of $\Omega \ ,\,\ \rho _0 >1$, in which $g_\varepsilon$ is extended $(C ,\,\ M \ \textrm {indep. of}\ \varepsilon )$.

And either

(ii) $g_\varepsilon \in C^{1,\alpha }(\overline {\Omega }) \ ,\,\ | g_{\varepsilon } |_{1,\alpha } \leq \frac {M}{\varepsilon }$ and $\partial \Omega$ is $C^2$, where we denote with $| \cdot |_{1, \alpha }$ as the $C^{1, \alpha }$ norm.

Or (ii’) $g_\varepsilon \in H^1( \Omega )$ and $J_\varepsilon ( {u_\varepsilon },\, \Omega ) \leq C$.

For $i \neq j \ ,\,\ i,\,j \in \lbrace 1,\,2,\,...,\,N \rbrace$, let $U \in W^{1,2}( \mathbb {R};\mathbb {R}^m)$ be the 1D minimizer of the action

(1.4)\begin{align} & \sigma_{ij}:= \min \int_{-\infty}^{+\infty} \left( \frac{1}{2} |U'|^2 + W(U) \right) \,{\rm d}t <{+}\infty\ , \nonumber\\ & \lim_{t \rightarrow - \infty} U(t) = a_i \ ,\ \lim_{t \rightarrow + \infty} U(t) = a_j \ ,\ U( \mathbb{R}) \in \mathbb{R}^m \setminus \lbrace W=0 \rbrace \end{align}

where $U$ is a connection that connects $a_i$ to $a_j \ ,\,\ i,\,j \in \lbrace 1,\,2,\,...,\,N \rbrace$.

The existence of such geodesics has been proved under minimal assumptions on the potential $W$ in [Reference Zuniga and Sternberg38].

Let $J_\varepsilon$ defined in (1.1), we define

(1.5)\begin{align} \tilde{J}_{\varepsilon}(u, \Omega) := \begin{cases} J_{\varepsilon}(u,\Omega) \ ,\ \textrm{if} \ u = g_{\varepsilon} \ \textrm{on} \ \Omega_{\rho_0} \setminus \Omega\ ,\ u \in H^1_{loc}(\mathbb{R}^n;\mathbb{R}^m) \\ + \infty\ , \ \textrm{otherwise} \end{cases} \end{align}

where $\Omega \subset \Omega _{\rho _0}$ as in (H2)(i) and let

(1.6)\begin{align} J_0(u, \Omega) := \sum_{1 \leq i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1} (\partial^* \Omega_i \cap \partial^* \Omega_j \cap \Omega) = \sum_{1 \leq i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1} (S_{ij}(u) \cap \Omega), \end{align}

where $S_{ij}(u):= \partial ^* \lbrace u=a_i \rbrace \cap \partial ^* \lbrace u=a_j \rbrace \ ,\,\ u \in BV(\Omega ;\lbrace a_1,\,a_2,\,...,\,a_N \rbrace )$ and we denote as $\partial ^* \Omega _k$ the reduced boundary of $\Omega _k$.

Finally, we define the limiting functional subject to the limiting boundary conditions

(1.7)\begin{align} \tilde{J}_0 (u, \Omega) := \begin{cases} J_0(u,\Omega) \ ,\ \textrm{if} \ u \in BV(\Omega;\lbrace a_1,a_2,...,a_N \rbrace) \quad \textrm{and}\ u= g_0 \quad \textrm{on} \quad \Omega_{\rho_0} \setminus \Omega\ \\ + \infty\ ,\ \textrm{otherwise} \end{cases} \end{align}

We can write $J_\varepsilon,\,J_0,\, \tilde {J}_\varepsilon,\, \tilde {J}_0 : L^1(\Omega ; \mathbb {R}^n) \rightarrow \overline {\mathbb {R}}$, where $\overline {\mathbb {R}} = \mathbb {R} \cup \lbrace \infty \rbrace$ and the $\Gamma$-convergence will be with respect to the $L^1$ topology.

Our first main result is the following

Theorem 1.1 Let $J_{\varepsilon }$ be defined by (1.1) and $\tilde {J}_{\varepsilon } \ ,\,\ \tilde {J}_0$ defined in (1.5) and (1.7) respectively.

Then

(1.8)\begin{equation} \Gamma- \lim_{\varepsilon \rightarrow 0} \tilde{J}_{\varepsilon}(u,\Omega) = \tilde{J}_0(u,\overline{\Omega}) . \end{equation}

Remark 1.2 Note that the domain of $\tilde {J}_0$ is the closure of $\Omega$, which means that there is a boundary term (see also (2.9) in [Reference Owen, Rubinstein and Sternberg32] for the analogue in the scalar case). More precisely, by proposition 3.5 and theorem 5.8 in [Reference Evans and Gariepy16] we can write

(1.9)\begin{align} & \tilde{J}_0(u, \overline{\Omega}) = \frac{1}{2} \sum_{i=1}^N \int_{\overline{\Omega}} | D ( \phi_i \circ u) | \ = \frac{1}{2} \sum_{i=1}^N \int_{\Omega} | D ( \phi_i \circ u) |\nonumber\\ & \qquad + \frac{1}{2} \sum_{i=1}^N \int_{\partial \Omega} | T(\phi_i \circ u) - T(\phi_i \circ g_0) | \ d \mathcal{H}^{n-1} \nonumber\\ & \quad \textrm{where} \ {\phi_i \ \textrm{defined in} \ (3.2) \ } \textrm{and} \ T \ \textrm{is the trace operator for} \ BV \ \textrm{functions}. \end{align}

The overview of the strategy of the proof of theorem 1.1 is as follows. First we observe that the $\Gamma$-limit established in [Reference Baldo7], in particular theorem 2.5, holds also without the mass constraint (see theorem 2.2 in Preliminaries section). Next, we apply a similar strategy to that of [Reference Ansini, Braides and Chiadò Piat6, Theorem 3.7] in which there is a $\Gamma$-convergence result with boundary conditions in the scalar case which states that we can incorporate the constraint of Dirichlet values in the $\Gamma$-limit, provided that this $\Gamma$-limit is determined. Since by theorem 2.2 we have that $J_\varepsilon \Gamma$-converges to $J_0$, we establish the $\Gamma$-limit of $\tilde {J}_\varepsilon$, that is, the $\Gamma$-limit of the functional $J_\varepsilon$ with the constraint of Dirichlet values. For the proof of the $\Gamma$-limit we can assume either (H2)(ii) or (H2)(ii’).

Next, we study the solution of the geometric minimization problem that arises from the limiting functional.

In order to obtain precise information about the minimizer of the limiting functional $\tilde {J}_0 (u,\, \overline {B}_1) \ ,\, \ B_1 \subset \mathbb {R}^2$, we impose that the limiting boundary conditions $g_0$ have connected phases. So we assume,

(H2) (iii) Let $g_0 = \sum _{i=1}^3 a_i \ \chi _{I_i}(\theta ) \ ,\,\ \theta \in [0,\, 2 \pi ) \ ,\,\ I_i \subset [0,\, 2 \pi ) \ ,\,\ \cup _{i=1}^3 I_i = [0,\, 2 \pi )$ be the limit of $g_\varepsilon$. Assume that $I_i$ are connected and that

\begin{align*} & \theta_0 < \frac{2 \pi}{3} \ ,\ \textrm{where} \ \theta_0 \ \textrm{is the largest angle of the points} \ p_i = \partial I_k \cap \partial I_l \\ & k \neq l \ ,\ i \in \lbrace 1,2,3 \rbrace \setminus \lbrace k,l \rbrace . \end{align*}

The assumption $\theta _0 < \frac {2 \pi }{3}$ arises from Proposition 3.2 in [Reference Morgan30] that we utilize for the proof (see proposition 2.5 in Preliminaries section) and guarantees that the boundary of the partition defined by the minimizer will be line segments meeting at a point inside $B_1$.

Our second main result is the following

Theorem 1.3 Let $u_0 = a_1 \chi _{\Omega _1} +a_2 \chi _{\Omega _2} + a_3 \chi _{\Omega _3}$ be a minimizer of $\tilde {J}_0(u,\,\overline {B}_1)$ subject to the limiting Dirichlet values (H2)(iii).

Then the minimizer is unique and in addition,

(1.10)\begin{equation} \partial \Omega_i \cap \partial \Omega_j \ \textrm{are line segments meeting at} \ 120^o \ \textrm{in a point in} \ B_1 \ (i \neq j). \end{equation}

For proving theorem 1.3, we first prove that the partition defined by $u_0$ is $(M,\,0,\, \delta )$-minimal as in Definition 2.1 in [Reference Morgan30] (see definition 2.4). This is proved by a comparison argument by defining a Lipschitz perturbation of the partition of the minimizer with strictly less energy. Then, by utilizing a uniqueness result for $(M,\, 0 ,\,\delta )$-minimal sets in [Reference Morgan30] (see proposition 2.5), we can conclude that the minimizer of the limiting energy is unique and the boundaries of the partition that the minimizer defines are line segments meeting at $120^o$ degrees in an interior point of the unit disc.

In the last subsection, we note that the result in theorem 1.3 can be extended also to the mass constraint case (see [Reference Baldo7]). However, in this case the uniqueness will be up to rigid motions of the disc (see Theorems 3.6 and 4.1 in [Reference Canete and Ritore10]).

1.2 Previous fundamental contributions

We will now briefly introduce some of the well-known results in the scalar case. The notion of $\Gamma$-convergence was introduced by De Giorgi and Franzoni in [Reference De Giorgi and Franzoni14] and in particular relates phase transition-type problems with the theory of minimal surfaces. One additional application of $\Gamma$-convergence is the proof of existence of minimizers of a limiting functional, say $F_0$, by utilizing an appropriate sequence of functionals $F_\varepsilon$ that we know they admit a minimizer and the $\Gamma$-limit of $F_\varepsilon$ is $F_0$. And also vice versa ([Reference Kohn and Sternberg25]), we can obtain information for the $F_\varepsilon$ energy functional from the properties of minimizers of the limiting functional $F_0$. We can think of this notion as a generalization of the Direct Method in the Calculus of Variations i.e. if $F_0$ is lower semicontinuous and coercive we can take $F_\varepsilon = F_0$ and then $\Gamma$-lim $F_\varepsilon =F_0$.

There are many other ways of thinking of this notion, such as a proper tool in finding the limiting functional among a sequence of functionals.

Let $X$ be the space of the measurable functions $u : \Omega \subset \mathbb {R}^n \rightarrow \mathbb {R}$ endowed with the $L^1$ norm and

\begin{align*} & F_\varepsilon(u,\Omega) := \begin{cases} \int_\Omega \frac{\varepsilon}{2} | \nabla u|^2 + \frac{1}{\varepsilon}W(u) \,{\rm d}x \ ,\ u \in W^{1,2} (\Omega ; \mathbb{R}) \cap X \ \\ + \infty\ , \ \textrm{elsewhere in} \ X \end{cases} \\ & F_0(u,\Omega) := \begin{cases} \sigma \mathcal{H}^{n-1} (Su) \ ,\ u \in SBV(\Omega ; \lbrace -1,1 \rbrace) \cap X \ \\ + \infty\ , \ \textrm{elsewhere in} \ X \end{cases} \ \\ & \quad \textrm{where} \ W : \mathbb{R} \rightarrow [0, + \infty ) \ ,\ \lbrace W =0 \rbrace = \lbrace -1 , 1 \rbrace\ ,\ \sigma = \int_{{-}1}^1 \sqrt{2 W(u)}\,{\rm d}u \\ & \qquad {and} \ Su \ \textrm{is the singular set of the SBV function} \ u. \ \end{align*}

Let now $u_\varepsilon$ be a minimizer of $F_\varepsilon$ subject to a mass constraint, that is, $\int _{\Omega }u = V \in (0,\, |\Omega |)$. The asymptotic behaviour of $u_\varepsilon$ was first studied by Modica and Mortola in [Reference Modica and Mortola28] and by Modica in [Reference Modica27, Reference Modica29]. Also, later Sternberg [Reference Sternberg34] generalized these results for minimizers with volume constraint. Furthermore, Owen et al. in [Reference Owen, Rubinstein and Sternberg32] and Ansini et al. in [Reference Ansini, Braides and Chiadò Piat6], among others, studied the asymptotic behaviour of the minimizers subject to Dirichlet values for the scalar case.

As mentioned previously, one of the most important outcomes of $\Gamma$-convergence in the scalar phase transition-type problems is the relationship with minimal surfaces. More precisely, the well-known theorem of Modica and Mortola states that the $\varepsilon$-energy functional of the Allen–Cahn equation $\Gamma$-converges to the perimeter functional that measures the perimeter of the interface between the phases (i.e. $\Gamma$-$\textrm {lim} \ F_\varepsilon = F_0$). So the interfaces of the limiting problem will be minimal surfaces.

This relationship is deeper as indicated in the De Giorgi conjecture (see [Reference De Giorgi15]) which states that the level sets of global entire solutions of the scalar Allen–Cahn equation that are bounded and strictly monotone with respect to $x_n$ are hyperplanes if $n \leq 8$. The relationship with the Bernstein problem for minimal graphs is the reason why $n \leq 8$ appears in the conjecture. The $\Gamma$-limit of the $\varepsilon$-energy functional of the Allen–Cahn equation is a possible motivation behind the conjecture.

In addition, Baldo in [Reference Baldo7] and Fonseca and Tartar in [Reference Fonseca and Tartar18] extended the $\Gamma$-convergence analysis for the phase transition-type problems to the vector case subject to a mass constraint and the limiting functional measures the perimeter of the interfaces separating the phases, and thus there is a relationship with the problem of minimizing partitions. In § 5 we analyse this in the set up of Dirichlet boundary conditions. Furthermore, the general vector-valued coupled case has been thoroughly studied in the works of Borroso–Fonseca and Fonseca–Popovici in [Reference Cristina Barroso and Fonseca12] and [Reference Fonseca and Popovici17] respectively.

There are many other fundamental contributions on the subject, such as the works of Gurtin [Reference Gurtin21, Reference Gurtin22], Gurtin and Matano [Reference Gurtin and Matano23] on the Modica–Mortola functional and its connection with materials science, the work of Hutchingson and Tonegawa on the convergence of critical points in [Reference Hutchinson and Tonegawa24], the work of Bouchitté [Reference Bouchitté8] and of Cristoferi and Gravina [Reference Cristoferi and Gravina13] on space-dependent wells and extensions on general metric spaces in the work of Ambrosio in [Reference Ambrosio5]. Several extensions to the non-local case and fractional setting have also been studied by Alberti-Bellettini in [Reference Alberti and Bellettini1], by Alberti-Bouchitté-Seppecher in [Reference Alberti, Bouchitté and Seppecher2] and by Savin-Valdinoci in [Reference Savin and Valdinoci33] among others.

2. Preliminaries

2.1 Specialized definitions and theorems for the $\Gamma$-limit

First, we will define the supremum of measures that allow us to express the limiting functional in an alternative way. Let $\mu$ and $\nu$ be two regular Borel measures on $\Omega$ we denote by $\mu \bigvee \nu$ the smallest regular positive measure which is greater than or equal to $\mu$ and $\nu$ on all borel subsets of $\Omega$, for $\mu \ ,\,\ \nu$ being two regular positive Borel measures on $\Omega.$ We have

\begin{align*} (\mu \bigvee \nu)(\Omega) := \sup \lbrace \mu (A) + \nu(B) : A \cap B = \emptyset, \ A\cup B \subset \Omega, \ A \quad \textrm{and} \\ B \ \textrm{are open sets in} \ \Omega \rbrace. \end{align*}

Now let

\begin{align*} \bigvee_{k=1}^N \int_{\Omega} | D (\phi_k \circ u_0 ) | := \sup \Bigg\{ \sum_{k=1}^N \int_{A_k} | D(\phi_k \circ u_0 ) | : \cup_{k=1}^N A_k \subset \Omega ,\\ A_i \cap A_j = \emptyset\ ,\ i \neq j, A_i \ \textrm{open sets in} \ \Omega \Bigg\}. \end{align*}

We will now provide a lemma from [Reference Ansini, Braides and Chiadò Piat6] that is crucial in the description of the behaviour of the $\Gamma$-limit with respect to the set variable. Let $\Omega \subset \mathbb {R}^n$ be an open set. We denote by $\mathcal {A}_{\Omega }$ the family of all bounded open subsets of $\Omega$.

Lemma 2.1 ([Reference Ansini, Braides and Chiadò Piat6]) Let $J_\varepsilon$ defined in (1.1). Then for every $\varepsilon >0$, for every bounded open set $U \ ,\,\ U' \ ,\,\ V$, with $U \subset \subset U'$, and for every $u,\,v \in L^1_{loc}(\mathbb {R}^n)$, there exists a cut-off function $\phi$ related to $U$ and $U'$, which may depend on $\varepsilon \ ,\,\ U \ ,\,\ U' \ ,\,\ V \ ,\,\ u \ ,\,\ v$ such that

\[ J_\varepsilon ( \phi u + (1- \phi) v, U \cup V) \leq J_\varepsilon (u, U') + J_\varepsilon (v,V) + \delta_\varepsilon (u,v,U,U',V), \]

where $\delta _\varepsilon : L^1_{loc}(\mathbb {R}^n)^2 \times \mathcal {A}^3_{\Omega } \rightarrow [0,\, + \infty )$ are functions depending only on $\varepsilon$ and $J_\varepsilon$ such that

\[ \lim_{\varepsilon \rightarrow 0} \delta_{\varepsilon} (u_\varepsilon, v_\varepsilon, U,U',V) =0, \]

whenever $U \ ,\,\ U' \ ,\,\ V \in \mathcal {A}_{\Omega } \ ,\,\ U \subset \subset U'$ and $u_\varepsilon \ ,\, \ v_\varepsilon \in L^1_{loc}(\mathbb {R}^n)$ have the same limit as $\varepsilon \rightarrow 0$ in $L^1( (U' \setminus \overline {U}) \cap V )$ and satisfy

\[ \sup_{\varepsilon>0}( J_\varepsilon (u_\varepsilon,U') + J_\varepsilon(v_\varepsilon,V)) <{+} \infty . \]

The above result is Lemma 3.2 in [Reference Ansini, Braides and Chiadò Piat6] and has been proved in the scalar case. The proof also works in the vector case with minor modifications. In [Reference Ansini, Braides and Chiadò Piat6], there is an assumption on $W$, namely $W \leq c( | u |^\gamma +1)$ with $\gamma \geq 2$ (see (2.2) in [Reference Ansini, Braides and Chiadò Piat6]). This assumption however is only utilized in the proof of lemma 2.1 above to apply the dominated convergence theorem in the last equation. In our case, this assumption is not necessary since $W(u_\varepsilon )$ and $W (g_\varepsilon )$ are uniformly bounded (see (H2)(i) and lemma 3.1). In fact, the only reason we assume in (H1) that $W(u) \geq c_1 | u |^2$ for $| u |>M$ is to apply the above lemma.

In [Reference Baldo7] it has been proved that $J_{\varepsilon } \ \Gamma$-converges to $J_0$ with mass constraint, but it also holds without mass constraint (see theorem 2.5). We will point out this more clearly in the proof of theorem 1.1. In particular, it holds

Theorem 2.2 [Reference Baldo7]

Let $J_\varepsilon$ defined in (1.1) and $J_0$ defined in (1.6). Then $\Gamma$-$\lim _{\varepsilon \rightarrow 0} J_{\varepsilon }(u,\,\Omega ) = J_0(u,\, \Omega )$ in $L^1(\Omega ; \mathbb {R}^m).$ That is, for every $u \in L^1(\Omega ; \mathbb {R}^m)$, we have the following two conditions:

  1. (i) If $\lbrace v_\varepsilon \rbrace \subset L^1(\Omega ; \mathbb {R}^m)$ is any sequence converging to $u$ in $L^1$, then

    (2.1)\begin{equation} \liminf_{\varepsilon \rightarrow 0} J_\varepsilon (v_\varepsilon, \Omega) \geq J_0 (u, \Omega), \end{equation}
    and
  2. (ii) There exists a sequence $\lbrace w_\varepsilon \rbrace \subset L^1(\Omega ; \mathbb {R}^m)$ converging to $u$ in $L^1$ such that

    (2.2)\begin{equation} \lim_{\varepsilon \rightarrow 0} J_\varepsilon (w_\varepsilon, \Omega) = J_0 (u, \Omega). \end{equation}

Remark 2.3 We note that in [Reference Baldo7], there is also a technical assumption for the potential $W$ (see (1.2) in p.70). However, for the proof of the $\Gamma$-limit this assumption is only utilized for the proof of the liminf inequality in order to obtain the equiboundedness of the minimizers $u_\varepsilon$ (see proof of (2.8) in [Reference Baldo7]). However, in our case we obtain equiboundedness from lemma 3.1 in the following section. Therefore, in our case this assumption is dismissed.

2.2 Specialized definitions and theorems for the geometric problem

In addition, we introduce the notion of $(M,\, 0 ,\,\delta )$-minimality as defined in [Reference Morgan30] together with a proposition that certifies the shortest network connecting three given points in $\mathbb {R}^2$ as uniquely minimizing in the context of $(M,\, 0 ,\,\delta )$-minimal sets. This characterization is one of the ingredients for the solution of the geometric minimization problem in the last section. In fact, in [Reference Morgan30] the more general notion of $(M,\, \varepsilon,\,\delta )$-minimality (or $(M,\, c r^\alpha,\,\delta )$-minimality) is introduced and regularity results for such sets are established. Particularly, $(M,\, 0 ,\,\delta )$-minimality implies $(M,\, c r^\alpha,\,\delta )$-minimality (see [Reference Morgan30]).

Definition 2.4 [Reference Morgan30]

Let $K \subset \mathbb {R}^n$ be a closed set and fix $\delta >0$. Consider $S \subset \mathbb {R}^n \setminus K$ be a nonempty bounded set of finite $m$-dimensional Hausdorff measure. $S$ is $(M,\, 0 ,\,\delta )$-minimal if $S = spt( \mathcal {H}^m \lfloor S) \setminus K$ and

\begin{align*} \mathcal{H}^m (S \cap W) \leq \mathcal{H}^m ( \phi(S \cap W)), \end{align*}

whenever

\begin{align*} & (a) \ \phi: \mathbb{R}^n \rightarrow \mathbb{R}^n \ \textrm{is lipschitzian}, \\ & (b) \ W = \mathbb{R}^n \cap \lbrace z \ : \ \phi(z) \neq z \rbrace , \\ & (c) \ \textrm{diam}(W \cup \phi(W)) < \delta , \\ & (d) \ \textrm{dist}(W \cup \phi(W), K) > 0. \end{align*}

Proposition 2.5 [Reference Morgan30]

Let $K = \lbrace p_1,\,p_2 ,\,p_3 \rbrace$ be the vertices of a triangle in the open $\delta$-ball $B(0,\,\delta ) \subset \mathbb {R}^2$, with largest angle $\theta$ for some fixed $\delta >0$. Then there exists a unique smallest $(M,\,0 ,\, \delta )$-minimal set in $B(0,\,\delta )$ with closure containing $K$, in particular:

\begin{align*} & (a) \ \textrm{if} \ \theta \geq 120^o,\ \textrm{the two shortest sides of the triangle} ; \\ & (b) \ \textrm{if} \ \theta < 120^o ,\ \textrm{segments from three vertices meeting at} \ 120^o . \end{align*}

Here by the ‘unique smallest’ we mean any other such $(M,\,0 ,\, \delta )$-minimal set $S$ has larger one-dimensional Hausdorff measure.

We now state a well-known Bernstein-type theorem in $\mathbb {R}^2$.

Theorem 2.6 [Reference Alikakos4]

Let $A$ be a complete minimizing partition in $\mathbb {R}^2$ with $N=3$ (three phases), with surface tension coefficients satisfying

(2.3)\begin{equation} \sigma_{ik} < \sigma_{ij} + \sigma_{jk} ,\ \textrm{for} \ j \neq i,k \ \textrm{with} \quad i,j,k \in \lbrace 1,2,3 \rbrace. \end{equation}

Then $\partial A$ is a triod.

For a proof and related material we refer to [Reference White37] and the expository [Reference Alikakos4].

3. Basic lemmas

Lemma 3.1 For every critical point $u_{\varepsilon } \in W^{1,2}(\Omega ;\mathbb {R}^m)$, satisfying (1.3) weakly together with assumptions (H1) and (H2)(i),(ii), it holds

\begin{align*} || u_{\varepsilon}||_{L^{\infty}} < M \quad \textrm{and} \quad || \nabla u_{\varepsilon}||_{L^{\infty}} < \frac{\tilde{C}}{\varepsilon}. \end{align*}

Proof. By linear elliptic theory, we have that $u_{\varepsilon } \in C^2(\Omega ;{\mathbb {R}}^m)$ (see e.g. Theorem 6.13 in [Reference Gilbarg and Trudinger19]). Set $v_{\varepsilon }(x) = |u_{\varepsilon }(x)|^2$, then

\[ \Delta v_{\varepsilon} = 2W_u(u_{\varepsilon}) \cdot u_{\varepsilon} + 2 | \nabla u_{\varepsilon}|^2 >0 \ \text{for} \ | u_{\varepsilon}| >M, \]

Hence $\max _{\Omega } | u_{\varepsilon }|^2 \leq M^2$.

On the other hand (from (H2)), $\max _{\partial \Omega } | u_{\varepsilon } | \leq M$. Thus, $\max _{\overline {\Omega }} | u_{\varepsilon }| \leq M .$

For the gradient bound, consider the rescaled problem $y= \frac {x}{\varepsilon }$, denote by $\tilde {u} \ ,\, \tilde {g}$ the rescaled $u_{\varepsilon } \ ,\, g_{\varepsilon }$, so by elliptic regularity (see e.g. Theorem 8.33 in [Reference Gilbarg and Trudinger19]),

\begin{align*} & | \tilde{u} |_{1,\alpha} \leq C (||\tilde{u}||_{L^{\infty}} + | \tilde{g}|_{1,\alpha}) \leq 2CM \\ & \Rightarrow || \nabla \tilde{u} ||_{L^{\infty}} \leq 2CM \Rightarrow | \nabla u_{\varepsilon} | \leq \frac{\tilde{C}}{\varepsilon} . \end{align*}

Lemma 3.2 Let $u_\varepsilon$ defined in (1.2), then

\[ J_{\varepsilon} (u_{\varepsilon}) = \int_{\Omega} \left( \frac{\varepsilon}{2} |\nabla u_{\varepsilon}|^2 + \frac{1}{\varepsilon} W(u_{\varepsilon}) \right) \,{\rm d}x \leq C \ , \]

$C$ independent of $\varepsilon >0$, if $\Omega$ is bounded.

Proof. Without loss of generality we will prove lemma 3.2 for $\Omega =B_1$ (or else we can cover $\Omega$ with finite number of unit balls and the outside part is bounded by (H2)(i)).

Substituting $y = \frac {x}{\varepsilon }$,

\[ J_{\varepsilon} (u_{\varepsilon}) = \int_{B_{\frac{1}{\varepsilon}}} \left( \frac{\varepsilon}{2} |\nabla_y \tilde{u}_{\varepsilon}|^2 \ \frac{1}{\varepsilon^2} + \frac{1}{\varepsilon} W(\tilde{u}_{\varepsilon}) \right) \varepsilon^n \,{\rm d}y , \]

where $\tilde {u}_{\varepsilon } = u_{\varepsilon } (\varepsilon y)$ and for $\varepsilon = \frac {1}{R}$,

\begin{align*} \Rightarrow J_{\varepsilon} (u_{\varepsilon}) & = \varepsilon^{n-1} \int_{B_{\frac{1}{\varepsilon}}} \left( \frac{1}{2} |\nabla_y \tilde{u}_{\varepsilon}|^2 + W(\tilde{u}_{\varepsilon}) \right) \,{\rm d}y\\ & = \frac{1}{R^{n-1}} \int_{B_R} \left( \frac{1}{2} |\nabla_y \tilde{u}_R|^2 + W(\tilde{u}_R) \right) \,{\rm d}y \ = \frac{1}{R^{n-1}} \tilde{J}_R (\tilde{u}_R) . \end{align*}

So, $\tilde {u}_R$ is minimizer of $\tilde {J}_R (v) = \int _{B_R} (\frac {1}{2} |\nabla v|^2 + W(v) ) \,{\rm d}x$.

By lemma 3.1 applied in $u_\varepsilon$, it holds that $| \tilde {u}_R | ,\,| \nabla \tilde {u}_R |$ are uniformly bounded independent of $R$ and via the comparison function (see [Reference Alikakos, Fusco and Smyrnelis3] p.135), for $R>1$

\begin{align*} v(x) := \begin{cases} a_1 , \ \textrm{for} \ | x | \leq R-1 \\ (R- | x|)a_1 + (|x| -R +1)\tilde{u}_R(x) ,\ \textrm{for} \ | x | \in (R-1,R] \\ \tilde{u}_R(x) , \ \textrm{for} \ |x| >R \end{cases} , \end{align*}

we have

\[ \tilde{J}_R (\tilde{u}_R) \leq J(v) \leq C R^{n-1} , \ C \ \textrm{independent of} \ R. \]

Thus

\[ J_{\varepsilon}(u_{\varepsilon}) = \frac{1}{R^{n-1}} \tilde{J}_R (\tilde{u}_R) \leq C \ ( C \ \textrm{independent of} \ \varepsilon >0 ).\]

Lemma 3.3 Let $u_\varepsilon$ defined in (1.2), then $u_{\varepsilon } \stackrel {L^1}{\longrightarrow } u_0$, along subsequences and $u_0 \in BV(\Omega ;\mathbb {R}^m)$. Moreover, $u_0 = \sum _{i=1}^N a_i \chi _{\Omega _i} \ ,\,\ \mathcal {H}^{n-1}(\partial ^* \Omega _i) < \infty$ and $| \Omega \setminus \cup _{i=1}^N \Omega _i | =0$.

Proof. By lemma 3.1 we have that $u_\varepsilon$ is equibounded. Now arguing as in the proof of Proposition 4.1 in [Reference Baldo7] (see also remark 2.3), we obtain that $|| u_{\varepsilon } ||_{BV(\Omega ;\mathbb {R}^m)}$ is uniformly bounded, $u_\varepsilon \rightarrow u_0$ in $L^1$ along subsequences and also $u_0 \in BV(\Omega ; \mathbb {R}^m) .$

From lemma 3.2, it holds

\[ \frac{1}{\varepsilon} \int_{\Omega} W(u_{\varepsilon}(x))\,{\rm d}x \leq C \ ( C \ \textrm{independent of} \ \varepsilon >0). \]

Since $| u_{\varepsilon } | \leq M$ and $W$ is continuous in $\overline {B}_M \subset {\mathbb {R}}^m \ \Rightarrow W( u_{\varepsilon }) \leq \tilde {M}$, by the dominated convergence theorem we obtain

\[ \int_{\Omega} W(u_0(x)) \,{\rm d}x =0 \Rightarrow u_0 \in \lbrace W=0 \rbrace\ a.e. \ \Rightarrow u_0 = \sum_{i=1}^N a_i \chi_{\Omega_i} \]

where $\chi _{\Omega _i}$ have finite perimeter since $u_0 \in BV( \Omega ;\mathbb {R}^m)$ (see [Reference Evans and Gariepy16]).

The proof of lemma 3.3 is complete.

Also, $g_0$ takes values on $\lbrace W=0 \rbrace$.

Lemma 3.4 Let $g_0$ be the limiting boundary condition of $g_\varepsilon$.

Then

\[ g_0 = \sum_{i=1}^N a_i \chi_{I_i} \ , \ \textrm{where} \ I_i \ \textrm{have finite perimeter and} \ | \partial \Omega \setminus \cup_{i=1}^N I_i|=0. \]

Proof. By (H2)(i),

\begin{align*} J_\varepsilon ( g_\varepsilon , \Omega_{\rho_0} \setminus \Omega) \leq C \\ \Rightarrow \frac{1}{\varepsilon} \int_{\Omega_{\rho_0} \setminus \Omega} W(g_\varepsilon) \,{\rm d}x \leq C \end{align*}

So, arguing as in the proof of lemma 3.3, we have that $g_0 \in \lbrace W=0 \rbrace$ and we conclude.

Proposition 3.5 It holds that

(3.1)\begin{align} & \int_{\Omega'} | D (\phi_k \circ u_0 ) | = \sum_{i=1, i\neq k}^N \sigma_{ik} \mathcal{H}^{n-1} ( \partial^* \Omega_k \cap \partial^* \Omega_i \cap \Omega' ) \nonumber\\ & k =1,2,..,N , \ \textrm{for every open} \ \Omega' \subset \Omega , \end{align}

where $\phi _k (z) = \,{\rm d}(z,\,a_k) \ ,\,\ k=1,\,2,\,...,\,N,\,$ and $a_k$ are the zeros of $W$ and $d$ is the Riemannian metric derived from $W^{1/2}$, that is

(3.2)\begin{align} & {\rm d}(z_1,z_2)\nonumber\\ & \quad := \inf \left\lbrace \int_0^1 \sqrt{2} W^{1/2} ( \gamma (t)) | \gamma' (t)| \,{\rm d}t : \gamma \in C^1 ([0,1];{\mathbb{R}}^2), \gamma(0) = z_1, \ \gamma(1)=z_2 \right\rbrace . \end{align}

Proof. The proof can be found in proposition 2.2 in [Reference Baldo7].

Furthermore, reasoning as in the proof of proposition 2.2 in [Reference Baldo7] we have,

(3.3)\begin{equation} \bigvee_{k=1}^N \int_{\Omega} | D (\phi_k \circ u_0 ) | = \sum_{1 \leq i< j \leq N} \sigma_{ij} \mathcal{H}^1 (\partial^* \Omega_i \cap \partial^* \Omega_j \cap \Omega) = J_0(u_0, \Omega). \end{equation}

The above equation is an alternative way to express the limiting functional.

4. Proof of the $\Gamma$-limit

Throughout the proof of the $\Gamma$-limit we will assume (H1) and (H2)(i),(ii). The proof if we assume (H2)(ii’) instead of (H2)(ii) is similar with minor modifications.

Proof Proof of theorem 1.1

We begin by proving the $\Gamma$-$\liminf$ inequality.

Let $u_\varepsilon \in L^1(\Omega ; \mathbb {R}^m)$ such that $u_\varepsilon \rightarrow u$ in $L^1(\Omega ; \mathbb {R}^m)$. If $u_\varepsilon \notin H^1_{loc}$ or $u_\varepsilon \neq g_\varepsilon$ on $\Omega _{\rho _0} \setminus \Omega$, where $\Omega \subset \Omega _{\rho _0}$ as in (H2)(i), then $\tilde {J}_\varepsilon (u_\varepsilon,\, \Omega ) = + \infty$ and the liminf inequality holds trivially. So, let $u_\varepsilon \in H^1_{loc}(\Omega ; \mathbb {R}^m)$ such that $u_\varepsilon \rightarrow u$ in $L^1$ and $u_\varepsilon = g_\varepsilon$ on $\Omega _{\rho _0} \setminus \Omega$.

Let $\rho >1$ such that $\rho < \rho _0$ in (H2)(i), we have

(4.1)\begin{equation} \tilde{J}_{\varepsilon}(u_{\varepsilon},\Omega) = J_{\varepsilon}(u_{\varepsilon}, \Omega_{\rho}) - J_{\varepsilon}(g_{\varepsilon}, \Omega_{\rho} \setminus \Omega), \end{equation}

where $\partial \Omega _\rho \in C^2$ since it is a small dilation of $\Omega$ and there is a unique normal vector $\nu \perp \partial \Omega _{\rho }$, such that each $x \in \partial \Omega$ can be written as $x = y + \nu (y) d \ ,\,\ d = dist(x,\, \partial \Omega _{\rho })$ (see the Appendix in [Reference Gilbarg and Trudinger19]).

So,

(4.2)\begin{equation} J_{\varepsilon}(g_{\varepsilon}, \Omega_{\rho} \setminus \Omega) = \int_1^\rho \int_{\partial \Omega_r} \left( \frac{\varepsilon}{2} | \nabla g_\varepsilon|^2 + \frac{1}{\varepsilon} W(g_\varepsilon) \right) dS dr \leq C ( \rho -1), \end{equation}

by Fubini's Theorem and (H2)(i).

Hence, by (4.1), for every $u_{\varepsilon }$ converging to $u$ in $L^1$ such that $u_{\varepsilon }= g_{\varepsilon }$ on $\Omega _{\rho _0} \setminus \Omega$ and $\liminf _{\varepsilon \rightarrow 0} \tilde {J}_{\varepsilon } (u_{\varepsilon },\, \Omega ) < + \infty$, we have that

(4.3)\begin{equation} \liminf_{\varepsilon \rightarrow 0} \tilde{J}_{\varepsilon}(u_{\varepsilon}, \Omega) \geq \liminf_{\varepsilon \rightarrow 0} J_{\varepsilon}(u_{\varepsilon},\Omega_{\rho}) - O(\rho-1). \end{equation}

Also, by the liminf inequality for $J_{\varepsilon }$ (see theorem 2.2 and (3.3)), we can obtain

(4.4)\begin{equation} \liminf_{\varepsilon \rightarrow 0} J_{\varepsilon}(u_{\varepsilon}, \Omega_\rho) \geq \sum_{1 \leq i< j \leq N} \sigma_{ij} \mathcal{H}^1 (\partial^* \Omega_i \cap \partial^* \Omega_j \cap \Omega_{\rho} ) = J_0(u, \Omega_{\rho}). \end{equation}

Thus, by (4.3) and (4.4), passing the limit as $\rho$ tends to $1$ we have the liminf inequality

(4.5)\begin{equation} \liminf_{\varepsilon \rightarrow 0} \tilde{J}_{\varepsilon}(u_{\varepsilon}, \Omega) \geq J_0 (u, \overline{\Omega}), \end{equation}

utilizing also the continuity of measures on decreasing sets.

We now prove the $\Gamma$-limsup inequality. Let $u \in BV(\Omega ; \lbrace a_1,\,a_2,\,...,\,a_N \rbrace )$ be such that $u = g_0$ on $\Omega _{\rho _0} \setminus \Omega. \\$ a) We first assume that $u = g_0$ on $\Omega \setminus \Omega _{\rho _1}$ with $\rho _1 < 1$ and $|\rho _1 -1|$ small.

As we observe in the proof of Theorem 2.5 in [Reference Baldo7] the $\Gamma$-limsup inequality for $J_\varepsilon$ also holds without the mass constraint, see in particular the proof of Lemma 3.1 in [Reference Baldo7]. Since the $\Gamma$-liminf inequality holds, the $\Gamma$-limsup inequality is equivalent with

(4.6)\begin{equation} J_0 (u, \Omega) = \lim_{\varepsilon \rightarrow 0} J_{\varepsilon}(u_{\varepsilon},\Omega ), \end{equation}

for some sequence $u_{\varepsilon }$ converging to $u$ in $L^1(\Omega ;\mathbb {R}^m)$. So let $u_\varepsilon$ be a sequence converging to $u$ in $L^1(\Omega _{\rho _1};\mathbb {R}^m)$ such that (4.6) is satisfied. In particular, $u_\varepsilon$ converges to $g_0$ on $\Omega \setminus \Omega _{\rho _1}$, where $\Omega _{\rho _1}$ is a small contraction of $\Omega$.

Now, utilizing the sequence $u_\varepsilon$ obtained from (4.6), we will modify it by a cut-off function so that the boundary condition is satisfied. By lemma 2.1, there exists a cut-off function $\phi$ between $U= \Omega _{\frac {1+\rho _1}{2}}$ and $U'= \Omega$ such that

(4.7)\begin{equation} J_{\varepsilon}(u_\varepsilon \phi + (1- \phi )g_\varepsilon, \Omega) \leq J_{\varepsilon}(u_\varepsilon, \Omega) + J_{\varepsilon} (g_\varepsilon, V) + \delta_{\varepsilon} (u_\varepsilon, g_\varepsilon, U,U',V), \end{equation}

where $V = \Omega \setminus \overline {\Omega }_{\rho _1}$ and $g_\varepsilon$ is extended in $V$ trivially.

By the assumptions on $u_\varepsilon$ and (H2) we also have

\[ u_\varepsilon \rightarrow g_0 ,\quad g_\varepsilon \rightarrow g_0 \ \textrm{in} \ L^1(V). \]

Hence, again by lemma 2.1 we get

\[ \lim_{\varepsilon \rightarrow 0} \delta_{\varepsilon} (u_\varepsilon, g_\varepsilon, U,U',V) =0. \]

Note that the condition $\sup _{\varepsilon >0}( J_\varepsilon (u_\varepsilon,\,U') + J_\varepsilon (g_\varepsilon,\,V)) < + \infty$ in lemma 2.1 is satisfied. To be more precise, from lemma 3.2 it holds

\[ \sup_{\varepsilon>0} J_{\varepsilon} (u_\varepsilon, U') <{+} \infty ,\ \textrm{where} \ U' = \Omega, \]

and by (H2)(i),

\[ \sup_{\varepsilon>0} J_{\varepsilon} (g_\varepsilon, V) <{+} \infty ,\ \textrm{where} \ V = \Omega \setminus \overline{\Omega}_{\rho_1}. \]

So, by (4.1), (4.2) and (4.7)

\begin{align*} \Gamma - \limsup_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon ( \tilde{u}_{\varepsilon}, \Omega) \leq \tilde{J}_0(u, \Omega), \end{align*}

where $\tilde {u}_{\varepsilon }= u_\varepsilon \phi + (1- \phi )g_\varepsilon$ and $\tilde {u}_{\varepsilon }= g_\varepsilon$ in $\Omega _{\rho _0} \setminus \Omega. \\$ b) In the general case we consider $\rho _1 <1$ and we define $u_{\rho _1}(x) = u(\frac {1}{\rho _1} x)$ and without loss of generality we may asume that the origin of $\mathbb {R}^n$ belongs in $\Omega$.

By the previous case (a) and (1.6),

(4.8)\begin{align} & \Gamma - \limsup_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon (u_{\rho_1}, \Omega) \leq \tilde{J}_0(u_{\rho_1}, \Omega) = \sum_{1 \leq i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1} (S_{ij}(u_{\rho_1}) \cap \Omega) \nonumber\\ & \quad\leq \sum_{1 \leq i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1} (S_{ij}(u) \cap \overline{\Omega}) + O(1 - \rho_1^{n-1}) \ = \tilde{J}_0 (u, \overline{\Omega}) + O(1- \rho_1^{n-1}). \end{align}

Since $u_{\rho _1}$ converges to $u$ as $\rho _1$ tends to $1$, if we denote

\[ J'(u_{\rho_1}, \Omega) := \Gamma - \limsup_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon (u_{\rho_1}, \Omega), \]

then by the lower semicontinuity of the $\Gamma$-upper limit (see e.g. Proposition 1.28 in [Reference Braides9]) and (4.8),

(4.9)\begin{equation} \Gamma - \limsup_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon (u_{\rho_1}, \Omega) \leq \liminf_{\rho_1 \rightarrow 1} J'(u_{\rho_1}, \Omega) \leq \tilde{J}_0 (u, \overline{\Omega}). \end{equation}

Hence, by (4.5) and (4.9) we get the required equality (1.8).

5. Minimizing partitions and the structure of the minimizer

In this section, we begin with the basic definitions of minimizing partitions. Then we underline the relationship of minimizing partitions in $\mathbb {R}^2$ with the minimizers of the functional $\tilde {J}_0$ and we analyse the structure of the minimizer of $\tilde {J}_0$ that we obtain from the $\Gamma$-limit. Utilizing a Bernstein-type theorem for minimizing partitions, we can explicitly compute the energy of the minimizer in proposition 5.5, and by regularity results in [Reference Morgan30], we can determine the precise structure of a minimizer subject to the limiting boundary conditions in theorem 1.3 and prove uniqueness. In subsection 5.2, we make some comments for the limiting minimizers in dimension three. Finally, in the last subsection, we note that we can extend these results to the mass constraint case.

Let $\Omega \subset \mathbb {R}^n$ open, occupied by $N$ phases. Associated to each pair of phases $i$ and $j$, there is a surface energy density $\sigma _{ij}$, with $\sigma _{ij}>0$ for $i \neq j$ and $\sigma _{ij} = \sigma _{ji}$, with $\sigma _{ii}=0$. Hence, if $A_i$ denoted the subset of $\Omega$ occupied by phase $i$, then $\Omega$ is the disjoint union

\[ \Omega = A_1 \cup A_2 \cup ... \cup A_N \]

and the energy of the partition $A = \lbrace A_i \rbrace _{i=1}^N$ is

(5.1)\begin{equation} E(A) = \sum_{1 \leq i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1}(\partial^* A_i \cap \partial^* A_j), \end{equation}

where $\mathcal {H}^{n-1}$ is the $(n-1)$-Hausdorff measure in $\mathbb {R}^n$ and $A_i$ are sets of finite perimeter. If $\Omega$ is unbounded, for example $\Omega =\mathbb {R}^n$ (we say then that $A$ is complete), the quantity above in general will be infinity. Thus, for each $W$ open, with $W \subset \subset \Omega$, we consider the energy

(5.2)\begin{equation} E(A;W)= \sum_{0< i < j \leq N} \sigma_{ij} \mathcal{H}^{n-1}(\partial^* A_i \cap \partial^* A_j \cap W). \end{equation}

Definition 5.1 The partition $A$ is a minimizing $N$-partition if given any $W \subset \subset \Omega$ and any $N$-partition $A'$ of $\Omega$ with

(5.3)\begin{equation} \bigcup_{i=1}^N (A_i \triangle A_i') \subset{\subset} W, \end{equation}

we have

\begin{align*} E(A;W) \leq E(A';W). \end{align*}

The symmetric difference $A_i \triangle A_i'$ is defined as their union minus their intersection, that is, $A_i \triangle A_i' = (A_i \cup A_i') \setminus (A_i \cap A_i')$.

To formulate the Dirichlet problem, we assume that $\partial \Omega$ is $C^1$ and given a partition $C$ of $\partial \Omega$ up to a set of $\mathcal {H}^{n-1}$-measure zero, we may prescribe the boundary data for $A$:

\begin{align*} (\partial_{\Omega} A)_i = \partial A_i \cap \partial \Omega = C_i,\quad i=1,...,N. \end{align*}

Now the energy is minimized subject to such a prescribed boundary.

Remark 5.2 Note that the minimization of the functional $\tilde {J}_0(u,\, \Omega )$ is equivalent to minimizing the energy $E(A;\Omega )$ under the appropriate Dirichlet conditions.

In figure 1 we show a triod with angles $\theta _1,\, \theta _2,\, \theta _3$, and the corresponding triangle with their supplementary angles $\hat {\theta }_i = \pi - \theta _i$ . For these angles Young's law holds, that is,

(5.4)\begin{equation} \frac{\sin \hat{\theta}_1}{\sigma_{23}} = \frac{\sin \hat{\theta}_2}{\sigma_{13}} = \frac{\sin \hat{\theta}_3}{\sigma_{12}}. \end{equation}

Definition 5.3 Let $\mathcal {A}_{x_0} = \lbrace A_1,\,A_2,\,A_3 \rbrace$ be a $3$-partition of $\mathbb {R}^2$ such that $A_i$ is a single infinite sector emanating from the point $x_0 \in \mathbb {R}^2$ with three opening angles $\theta _i$ that satisfy (5.4). We call as a triod $C_{tr}(x_0)$ the boundary of the partition $\mathcal {A}_{x_0}$, that is, $C_{tr}(x_0) = \lbrace \partial A_i \cap \partial A_j \rbrace _{1 \leq i < j \leq 3}$.

Figure 1. In the left we show a triod with angles $\theta_1, \theta_2, \theta_3$. In the right there is the corresponding triangle with supplementary angles $\hat{\theta}_1, \hat{\theta}_2, \hat{\theta}_3$ that satisfy the Young's law.

So, in other words, the triod is consisted of three infinite lines meeting at a point $x_0$ and their angles between the lines satisfy the Young's law (5.4) (see figure 1). As we see in theorem 2.6, the triod is the unique locally $3$-minimizing partition of $\mathbb {R}^2$. The point $x_0$, i.e. the centre of the triod, is often called a triple junction point.

5.1 The structure of the minimizer in the disk

Throughout this section, we will assume that $\sigma _{ij} = \sigma >0$ for $i \neq j$, therefore we have by Young's law $\theta _i = \frac {2 \pi }{3} \ ,\,\ i=1,\,2,\,3$. As a result of theorem 2.6, we expect that, by imposing the appropriate boundary conditions, the minimizer $u_0$ of $\tilde {J}_0(u,\, \overline {B}_1) \ ,\, B_1 \subset \mathbb {R}^2$ which we obtain from the $\Gamma$-limit will be a triod with angles $\frac {2 \pi }{3}$ restricted in $B_1$ and centred at a point $x \in B_1$.

We now recall Steiner's problem that gives us some geometric intuition about this fact.

Let us take three points $A ,\, \ B$ and $C$, arranged in any way in the plane. The problem is to find a fourth point $P$ such that the sum of distances from $P$ to the other three points is a minimum; that is, we require $AP + BP + CP$ to be a minimum length.

If the triangle $ABC$ possesses internal angles which are all less than $120^o$, then $P$ is the point such that each side of the triangle, i.e. $AB ,\,\ BC$ and $CA$, subtends an angle of $120^o$ at $P$. However, if one angle, say $A \hat {C} B$, is greater than $120^o$, then $P$ must coincide with $C .$

The Steiner's problem is a special case of the geometric median problem and has a unique solution whenever the points are not collinear. For more details and proofs, see [Reference Gilbert and Pollak20].

The problem of minimizing partitions subject to boundary conditions, in contrast to the mass constraint case, might not always admit a minimum, we provide an example in figure 2 below.

Figure 2. The geometric problem subject to such boundary conditions does not admit a minimum. However, the limiting functional admits a minimizer that forms a boundary layer.

However, a minimizer will exist for the minimization problem $\min _{u\in BV(\Omega ; \lbrace W=0 \rbrace )} \tilde {J}_0(u,\, \overline {\Omega })$, for instance the one we obtain from the $\Gamma$-limit, which will form a ‘boundary layer’ in the boundary of the domain instead of internal layer (i.e. the interface separating the phases). Particularly, in figure 2 above, $u_0 = a_1$, a.e. will be a minimizer of $\tilde {J}_0$ and

\[ \tilde{J}_0(u_0,\overline{\Omega}) = \frac{1}{2} \sum_{i=1}^3 \int_{\partial \Omega} | T(\phi_i \circ u_0) - T ( \phi_i \circ g_0) | d\mathcal{H}^1 = \sigma \mathcal{H}^1(\partial \Omega_{AB}) , \]

where $\partial \Omega _{AB}$ is the part of the boundary of $\Omega$ in which $g_0 = a_2$. When there are no line segments in the boundary of the domain or when $g_0$ does not admit jumps nearby such line segments, then we expect that there are no boundary layers and the boundary term in the energy of $\tilde {J}_0$ vanishes (see remark 1.2), otherwise we could find a minimizer with strictly less energy. In the cases where the boundary term vanishes we can write $\tilde {J}_0(u_0,\, \overline {\Omega }) = \tilde {J}_0(u_0,\,\Omega )$. This can be proved rigorously in the case where $\Omega = B_1$ and assuming (H2)(iii), utilizing also proposition 2.5 as we will see in the proof of theorem 1.3.

Remark 5.4 For the mass constraint case, by classical results of Almgren's improved and simplified by Leonardi in [Reference Leonardi26] for minimizing partitions with surface tension coefficients $\sigma _{ij}$ satisfying the strict triangle inequality (see (2.3)), $\Omega _j$ can be taken open with $\partial \Omega _j$ real analytic except possibly for a singular part with Hausdorff dimension at most $n-2$. Therefore, $\partial ^* \Omega _i \cap \partial ^* \Omega _j = \partial \Omega _i \cap \partial \Omega _j \ ,\,\ \mathcal {H}^{n-1}$-a.e., where $u_0 = \sum _{i=1}^N a_i \chi _{\Omega _i}$ is the minimizer of $J_0$ with a mass constraint. These regularity results have been stated by White in [Reference White36] but without providing a proof. Also, Morgan in [Reference Morgan31] has proved regularity of minimizing partitions in the plane subject to mass constraint. However, we deal with the problem with boundary conditions, so we cannot apply these regularity results.

Notation: We set as $x_0 \in B_1$ the point such that the line segments starting from $p_i = \partial I_k \cap \partial I_l \ ,\,\ k \neq l \ ,\,\ i \in \lbrace 1,\,2,\,3 \rbrace \setminus \lbrace k,\,l \rbrace$ and ending at $x_0$ meet all at angle $\frac {2 \pi }{3}$ (see (H2)(iii) and proposition 2.5). Also we denote by $C_0$ the sum of the lengths of these line segments. The following proposition measures the energy of the limiting minimizer.

Proposition 5.5 Let $(u_\varepsilon )$ be a minimizing sequence of $\tilde {J}_\varepsilon (u,\,B_1)$. Then $u_\varepsilon \rightarrow u_0$ in $L^1$ along subsequence with $u_0 \in BV(B_1; \lbrace a_1,\,a_2,\,a_3 \rbrace )$ and $u_0$ is a minimizer of $\tilde {J}_0(u,\, \overline {B}_1)$ subject to the limiting Dirichlet values (H2)(iii), where we extend $u$ by setting $u=g_0$ on $\mathbb {R}^2 \setminus B_1$.

In addition, we have

(5.5)\begin{equation} \sum_{1 \leq i < j \leq 3} \mathcal{H}^1 (\partial^* \Omega_i \cap \partial^* \Omega_j \cap \overline{B}_1 ) = C_0 \ , \end{equation}

where $u_0 = a_1 \chi _{\Omega _1} +a_2 \chi _{\Omega _2} + a_3 \chi _{\Omega _3}.$

Proof. From lemmas 3.2 and 3.3, it holds that if $u_\varepsilon$ is a minimizing sequence for $\tilde {J}_\varepsilon (u,\,B_1)$, then $\tilde {J}_\varepsilon (u_\varepsilon,\,B_1) \leq C$ and thus $u_\varepsilon \rightarrow u_0$ in $L^1$ along subsequence. The fact that $u_0$ is a minimizer of $\tilde {J}_0$ is a standard fact from the theory of $\Gamma$-convergence. It can be seen as follows.

Let $w \in BV(\overline {B_1},\, \lbrace a_1,\,a_2,\,a_3 \rbrace )$ such that $w = g_0$ on $\mathbb {R}^2 \setminus B_1$, then from the limsup inequality in theorem 1.1, we have that there exists $w_\varepsilon \in H^1_{loc}(\mathbb {R}^2 ;\mathbb {R}^m) \ ,\,\ w_\varepsilon = g_\varepsilon$ on $\mathbb {R}^2 \setminus B_1$ such that $w_\varepsilon \rightarrow w$ in $L^1$ and $\limsup _{\varepsilon \rightarrow 0} \tilde {J}_\varepsilon (w_\varepsilon,\,B_1) \leq \tilde {J}_0 (w,\, \overline {B}_1)$. Now since $u_\varepsilon$ is a minimizing sequence for $\tilde {J}_\varepsilon (u,\,B_1)$ and from the liminf inequality in theorem 1.1, we have

(5.6)\begin{align} \tilde{J}_0 (u_0 , \overline{B}_1) & \leq \liminf_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon (u_\varepsilon, B_1) \leq \liminf_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon ( w_\varepsilon,B_1) \nonumber\\ & \leq \limsup_{\varepsilon \rightarrow 0} \tilde{J}_\varepsilon ( w_\varepsilon,B_1) \leq \tilde{J}_0 (w, \overline{B}_1) \end{align}

For proving (5.5), we utilize theorem 2.6 (i.e. Theorem 2 in [Reference Alikakos4]). Since the triod is a minimizing 3-partition in $\mathbb {R}^2$ we have that for any $W \subset \subset \mathbb {R}^2$ and any partition it holds that $E(A,\,W) \leq E(V,\,W)$, where suppose that $A = \lbrace A_1,\,A_2,\,A_3 \rbrace$ is the partition of the triod and $V = \lbrace V_1,\,V_2,\,V_3 \rbrace$ is a 3-partition in $\mathbb {R}^2 .$

We have $u_0 = a_1 \chi _{\Omega _1} +a_2 \chi _{\Omega _2} + a_3 \chi _{\Omega _3}$ such that $u_0 = g_0$ on $\partial B_1$ and extend $u_0$ in $\mathbb {R}^2$, being the triod with $\theta _i = \frac {2 \pi }{3}$ in $\mathbb {R}^2 \setminus B_1$ centred at $x_0$. This defines a 3-partition in $\mathbb {R}^2$, noted as $\tilde {\Omega } = \lbrace \tilde {\Omega }_i \rbrace _{i=1}^3$. Since the triod is a minimizing 3-partition in the plane, we take any $W \subset \subset \mathbb {R}^2$ such that $B_2 \subset \subset W$ and $\bigcup _{i=1}^3 (A_i \triangle \tilde {\Omega }_i) \subset \subset W$, so we have

(5.7)\begin{equation} E(A, W) = E(A,\overline{B}_1) + E(A , W \setminus \overline{B}_1) \leq E( \tilde{\Omega},W) = E( \tilde{\Omega} , \overline{B}_1) + E( \tilde{\Omega} ,W \setminus \overline{B}_1) \end{equation}

where $A$ is the partition of the triod.

Now since

\[ E(A , W \setminus \overline{B}_1) = E( \tilde{\Omega} ,W \setminus \overline{B}_1) \]

from the way we extended $u_0$ in $\mathbb {R}^2$ and

\[ E(A , \overline{B}_1) = \sigma \sum_{1 \leq i < j \leq 3} \mathcal{H}^1 ( \partial A_i \cap \partial A_j \cap \overline{B}_1) = C_0 \sigma \]

since $\partial A_i \cap \partial A_j \cap \overline {B}_1$ are line segments inside $B_1$ with sum of their lengths equals $C_0$, we conclude

(5.8)\begin{align} C_0 \sigma & \leq E( \tilde{\Omega} , \overline{B}_1) = \tilde{J}_0 (u_0, \overline{B}_1) \nonumber\\ & \Leftrightarrow C_0 \leq \sum_{1 \leq i < j \leq 3} \mathcal{H}^1(\partial^* \Omega_i \cap \partial^* \Omega_j \cap \overline{B}_1) \end{align}

For the upper bound inequality $\sum _{1 \leq i < j \leq 3} \mathcal {H}^1(\partial ^* \Omega _i \cap \partial ^* \Omega _j \cap \overline {B}_1) \leq C_0$, we consider as a comparison function $\tilde {u}= a_1 \chi _{A_1} +a_2 \chi _{A_2} + a_3 \chi _{A_3}$, where $C_{tr}(x_0)= \lbrace A_1,\,A_2,\,A_3 \rbrace$ is the partition of the triod centred at $x_0 \in B_1$ and angles $\theta _i = \frac {2 \pi }{3}$ (see definition 5.3).

Then $\tilde {u}$ satisfies the boundary condition $\tilde {u} = g_0$ on $\mathbb {R}^2 \setminus B_1$ and therefore by the minimality of $u_0$ we have

(5.9)\begin{align} \tilde{J}_0 (u_0, \overline{B}_1) & \leq \tilde{J}_0 ( \tilde{u}, \overline{B}_1) = C_0 \sigma \nonumber\\ & \Rightarrow \sum_{1 \leq i < j \leq 3} \mathcal{H}^1(\partial^* \Omega_i \cap \partial^* \Omega_j \cap \overline{B}_1) \leq C_0. \end{align}

Corollary 5.6 Assume for simplicity that $x_0$ in proposition 5.5 above is the origin of $\mathbb {R}^2$. Then for every $R>0$ the energy of the limiting minimizer will satisfy

(5.10)\begin{equation} \tilde{J}_0 (u_0, \overline{B}_R) = 3 \sigma R. \end{equation}

In addition, there exists an entire minimizer in the plane and the partition that defines is a minimal cone.

Proof. Since $x_0$ is the origin of $\mathbb {R}^2$, it holds that $C_0$ in (5.5) equals $3$. Arguing as in proposition 5.5 above we can similarly obtain a minimizer of $\tilde {J}_0 (u_0,\, \overline {B}_R)$ that satisfies (5.10). By a diagonal argument the minimizer can be extended in the entire plane and will also satisfy

\[ \frac{\mathcal{H}^1(\partial \Omega_i \cap \partial \Omega_j \cap B_R)}{\omega_1 R} = C ,\quad \forall\ R >0. \]

Thus, the partition that it defines is a minimal cone (see [Reference White37] or [Reference Alikakos4]).

Finally, we will prove that the minimizer of $\tilde {J}_0$ in $\overline {B}_1$ is unique, that is, the only minimizer is the triod restricted to $B_1$ centred at a point in $B_1 .$ In figure 3 below we provide the structure of the minimizer $u_0$ obtained in theorem 1.3.

Proof Proof of theorem 1.3

Firstly, we show that the minimizing partition of $B_1$ with respect to the boundary conditions defined from $g_0$ is a $(M,\,0,\,\delta )$-minimal for $\delta >0$ (see definition 2.4). If not, let $S$ be the partition defined from $u_0$, we can find a Lipschitz function $\phi : \mathbb {R}^2 \rightarrow \mathbb {R}^2$ such that

\[ \mathcal{H}^1 (S \cap W) > \mathcal{H}^1( \phi(S \cap W)), \]

with

\begin{align*} W & = \mathbb{R}^2 \cap \lbrace x : \phi(x) \neq x \rbrace\ ,\ \textrm{diam}(W \cup \phi(W)) < \delta \\ & \quad \textrm{and dist}(W \cup \phi(W), \mathbb{R}^2 \setminus B_1) >0. \end{align*}

So if we consider the partition

\[ \tilde{S} := \begin{cases} S \ ,\ S \cap W = \emptyset \\ \phi(S \cap W) \ ,\ S \cap W \neq \emptyset \end{cases} , \]

then the boundary of the partition defined by $\tilde {S}$ will satisfy the boundary conditions (since dist$(W \cup \phi (W),\, \mathbb {R}^2 \setminus B_1) >0$) and also $\mathcal {H}^1 (\tilde {S}) < \mathcal {H}^1(S)$ which contradicts the minimality of $S .$

Thus, by (H2)(iii) we apply proposition 2.5 and we have that the unique smallest $(M,\,0,\, \delta )$-minimal set consists of three line segments from the three vertices defined from $g_0$ (i.e. the jump points in $\partial B_1$) meeting at $\frac {2 \pi }{3}$. The meeting point is unique and belongs in the interior of $B_1$. Thus, $\partial \Omega _i \cap \partial \Omega _j = \partial ^* \Omega _i \cap \partial ^* \Omega _j$ are line segments meeting at $\frac {2 \pi }{3}$ in an interior point of $B_1$.

Corollary 5.7 Let $u_0 = a_1 \chi _{\Omega _1} +a_2 \chi _{\Omega _2} + a_3 \chi _{\Omega _3}$ be a minimizer of $\tilde {J}_0(u,\,\overline {B}_1)$ subject to the limiting Dirichlet values $g_0( \theta ) = a_1 \chi _{(0, \frac {2 \pi }{3})} + a_2 \chi _{(\frac {2 \pi }{3} ,\, \frac {4 \pi }{3})} +a_3 \chi _{(\frac {4 \pi }{3} ,\, 2 \pi )}, \theta \in (0,\, 2\pi ) .$ Then $\partial \Omega _i \cap \partial \Omega _j$ are radi of $B_1 \ ,\,\ | \Omega _i |= \frac {1}{3}| B_1 |$ and the minimizer is unique.

Figure 3. Here is an example of a minimizer that we obtain in theorem 1.3.

In figure 4 above we illustrate the structure of the minimizer $u_0$ obtained in corollary 5.7.

Figure 4. The singular set of the minimizer obtained in corollary 5.7 is consisted of three radii of the ball.

5.2 Minimizers in dimension three

In this subsection, we will briefly make some comments for the structure of minimizers in $\mathbb {R}^3$. If we impose the appropriate boundary conditions in $B_R \subset \mathbb {R}^3$ and $\lbrace W=0 \rbrace = \lbrace a_1,\,a_2,\,a_3 \rbrace \ ,\,\ g_\varepsilon \rightarrow g_0 \ \textrm {in} \ L^1(B_R; \mathbb {R}^3)$ such that the partition in $\partial B_R$ defined by $g_0$ is equal to the partition of $(C_{tr} \times \mathbb {R}) \cap \partial B_R$, where $C_{tr}$ is the triod as in figure 1 (with equal angles), then by Theorem 3 in [Reference Alikakos4], arguing as in proposition 5.5 (see also corollary 5.6), we can obtain

\[ \tilde{J}_0 (u, B_R) = \frac{3}{2} \sigma \pi R^2 \ , \]

which gives

\[ \frac{ \mathcal{H}^2 ( \partial \Omega_i \cap \partial \Omega_j \cap B_R)}{\omega_2 R^2} = \frac{3}{2} \ , \]

where $\omega _2$ is the volume of the 2-dimensional unit ball (see [Reference White37]). That is, the partition that the minimizer defines can be extended to a minimal cone in $\mathbb {R}^3$. Now since the only minimizing minimal cones are the triod and the tetrahedral cone (see [Reference Taylor35]), then the minimizer of $\tilde {J}_0$ is such that $u_0 = \sum _{i=1}^3 a_i \chi _{\Omega _i}$, where $\Omega = \lbrace \Omega _i \rbrace _{i=1}^3$ is the partition of $(C_{tr} \times \mathbb {R} ) \cap B_R$.

Similarly, if $\lbrace W=0 \rbrace = \lbrace a_1,\,a_2,\,a_3,\,a_4 \rbrace$ and we impose the Dirichlet conditions such that $g_0$ defines the partition of the tetrahedral cone intersection with $\partial B_R$, then again $u_0 = \sum _{i=1}^4 a_i \chi _{\Omega _i}$, where $\Omega = \lbrace \Omega _i \rbrace _{i=1}^4$ is the partition of the tetrahedral cone restricted in $B_R$.

5.3 Minimizers in the disc for the mass constraint case

Throughout this subsection, we will assume that $a_i \ ,\, i=1,\,2,\,3,\,$ are affinely independent, that is, they are not contained in a single line. This can also be expressed as

(5.11)\begin{equation} \textrm{whenever} \ \sum_{i=1}^3 a_i \lambda_i =0 \ \textrm{with} \ \sum_{i=1}^3 \lambda_i =0 \ ,\ \textrm{then} \ \lambda_i =0 ,\quad i=1,2,3. \end{equation}

In addition, we consider that $m =(m_1,\,m_2) \in \mathbb {R}^2$ such that $m_1,\,m_2 >0$ (as in [Reference Baldo7]).

Let $u_0$ be a minimizer of $J_0(u,\,B_1) \ ,\,\ B_1 \subset \mathbb {R}^2$ defined in (1.6) subject to the mass constraint

(5.12)\begin{equation} \int_{B_1} u(x) \,{\rm d}x = m, \end{equation}

(i.e. the minimizer $u_0$ of Theorem p.70 in [Reference Baldo7]) and $\lbrace W=0 \rbrace = \lbrace a_1,\,a_2,\,a_3 \rbrace$. Then $u_0 = \sum _{i=1}^3 a_i \chi _{\Omega _i}$, where $\Omega _1,\,\Omega _2,\,\Omega _3$ is a partition of $B_1$ which minimizes the quantity

(5.13)\begin{equation} \sum_{1 \leq i< j \leq 3} \sigma \mathcal{H}^1 (\partial^* \Omega_i \cap \partial^* \Omega_j), \end{equation}

among all other partitions of $B_1$ such that $\sum _{i=1}^3 | \Omega _i | a_i = m .$

Theorem 5.8 Let $u_0$ be a minimizer of $J_0(u,\,B_1)$ as above and assume that

(5.14)\begin{equation} m = \sum_{i=1}^3 c_i a_i ,\ \textrm{where} \ c_i >0 ,\ \textrm{with} \ \sum_{i=1}^3 c_i = | B_1 |. \end{equation}

Then

(5.15)\begin{align} & | \Omega_i | = c_i ,\ i=1,2,3 , \partial^* \Omega_i \cap \partial^* \Omega_j = \partial \Omega_i \cap \partial \Omega_j \ \textrm{are piecewise smooth} \nonumber\\ & \quad \textrm{and the minimizer is unique up to a rigid motion of the disc}. \end{align}

In particular, the boundary of the partition is consisted of three circular arcs or line segments meeting at an interior vertex at $120$ degrees angles, reaching orthogonally $\partial B_1$ and so that the sum of geodesic curvature is zero.

Proof. We have that $u_0 = \sum _{i=1}^3 a_i \chi _{\Omega _i}$, where $\Omega _i$ are such that $\sum _{i=1}^3 | \Omega _i | = | B_1 |$ and $u_0$ minimizes the quantity (5.13).

By assumption (5.14), since $u_0$ satisfies (5.12), we have

(5.16)\begin{align} \sum_{i=1}^3 a_i | \Omega_i | & = \sum_{i=1}^3 c_i a_i \ \textrm{and} \ \sum_{i=1}^3 ( | \Omega_i | -c_i) = 0 \nonumber\\ & \Rightarrow | \Omega_i | = c_i \ ,\ i=1,2,3 \ ,\ \textrm{and} \ c_i \in (0, | B_1 |), \end{align}

since $a_i$ are affinely independent.

Now by Theorem 4.1 in [Reference Canete and Ritore10] we conclude that the minimizer is a standard graph i.e. it is consisted of three circular arcs or line segments meeting at an interior vertex at $120$ degrees angles, reaching orthogonally $\partial B_1$ and so that the sum of geodesic curvature is zero. So, $\partial ^* \Omega _i \cap \partial ^* \Omega _j = \partial \Omega _i \cap \partial \Omega _j$ are piecewise smooth.

Finally, the minimizer is unique up to rigid motions of the disc by Theorem 3.6 in [Reference Canete and Ritore10].

Note that in the case where $m = \frac {1}{3} | B_1 | \sum _{i=1}^3 c_i a_i$, it holds that $| \Omega _i | = \frac {1}{3} | B_1 | \ ,\,\ i=1,\,2,\,3 \ ,\,\ \textrm {and} \partial \Omega _i \cap \partial \Omega _j$ are line segments meeting at the origin and the minimizer is unique up to rotations.

Acknowledgements

I wish to thank my advisor Professor Nicholas Alikakos for his guidance and for suggesting this topic as a part of my thesis for the Department of Mathematics and Applied Mathematics at the University of Crete. Also, I would like to thank Professor P. Sternberg and Professor F. Morgan for their valuable comments on a previous version of this paper, which let to various improvements. Finally, I would like to thank the anonymous referee for their valuable suggestions, which not only enhanced the presentation but also significantly improved the quality of the paper by relaxing some of the assumptions in our results.

References

Alberti, G. and Bellettini, G.. A non-local anisotropic model for phase transitions: asymptotic behaviour of rescaled energies. Eur. J. Appl. Math. 9 (1998), 261284.Google Scholar
Alberti, G., Bouchitté, G. and Seppecher, P.. Un résultat deperturbations singulières avec la norme $H^1/2$. C.R.Acad. Sci. Paris – Serie I 319 (1994), 333338.Google Scholar
Alikakos, N. D., Fusco, G. and Smyrnelis, P.. Elliptic systems of phase transition type. Progress in Nonlinear Differential Equations and their Applications, (Birkhauser, 2018).CrossRefGoogle Scholar
Alikakos, N. D.. On the structure of phase transition maps for three or more coexisting phases, Geometric Partial Differential Equations proceedings (2013), pp. 1–31.Google Scholar
Ambrosio, L.. Metric space valued functions of bounded variation. Ann. Sc. Norm. Super. Pisa Cl. Sci. 4 17 (1990), 439478.Google Scholar
Ansini, N., Braides, A. and Chiadò Piat, V.. Gradient theory of phase transitions in composite media. P. Roy. Soc. Edinb. 133A (2003), 265296.Google Scholar
Baldo, S.. Minimal interface criterion for phase transitions in mixtures of Cahn–Hilliard fluids. Ann. Inst. Henri Poincare (C) Anal. Non Lineaire 7 (1990), 6790.Google Scholar
Bouchitté, G.. Singular perturbations of variational problems arising from a two-phase transition model. Appl. Math. Optim. 21 (1990), 289314.CrossRefGoogle Scholar
Braides, A.. $\Gamma$-convergence for beginners, Oxford Lecture Series in Mathematics and its Applications.Google Scholar
Canete, A. and Ritore, M.. Least-perimeter partitions of the disk into three regions of given areas. Indiana Univ. Math. J. 53 (2004), 883904.Google Scholar
Chan, C. C.. Structure of the singular set in energy-minimizing partitions and area-minimizing surfaces in $\mathbb {R}^N$, PhD thesis, Stanford University, 1995.Google Scholar
Cristina Barroso, A. and Fonseca, I.. Anisotropic singular perturbations—the vectorial case. Proc. Roy. Soc. Edinb. Sec. A. Math. 124 (2011), 527571.Google Scholar
Cristoferi, R. and Gravina, G.. Sharp interface limit of a multi-phase transitions model under nonisothermal conditions. Calc. Var. Partial. Differ. Equ. 60 (2021), 142.CrossRefGoogle Scholar
De Giorgi, E. and Franzoni, T.. Su un tipo di convergenza variazionale. Atti della Accad. Nazionale dei Lincei. Cl. Sci. Fis. Mat. Natur. Rend. Ser. 58 (1975), 842850.Google Scholar
De Giorgi, E.. Convergence problems for functionals and operators, Proceedings of the International Meeting on Recent Methods in Non linear Analysis (1978), pp.131–188.Google Scholar
Evans, L. and Gariepy, R.. Measure Theory and Fine Properties of Functions (New York: Chapman and Hall/CRC, 2015).Google Scholar
Fonseca, I. and Popovici, C.. Coupled singular perturbations for phase transitions. Asymptot. Anal. 44 (2005), 299325.Google Scholar
Fonseca, I. and Tartar, L.. The gradient theory of phase transitions for systems with two potential wells. Proce. Roy. Soc. Edinb. Sec. A. Math. 111 (1989), 89102.Google Scholar
Gilbarg, D. and Trudinger, N.. Elliptic Partial Differential Equations of Second Order, 2nd ed, Classics in Mathematics (Springer).Google Scholar
Gilbert, E. N. and Pollak, H. O.. Steiner minimal trees. SIAM. J. Appl. Math. 16 (1968), 129.Google Scholar
Gurtin, M. E.. On phase transitions with bulk, interfacial, and boundary energy. Arch. Ration. Mech. Anal. 96 (1986), 243264.Google Scholar
Gurtin, M. E.. Some results and conjectures in the gradient theory of phase transitions. Metastability Incompletely Posed Probl. 3 (1987), 135146.CrossRefGoogle Scholar
Gurtin, M. E. and Matano, H.. On the structure of equilibrium phase transitions within the gradient theory of fluids. Q. Appl. Math. 46 (1988), 301317.CrossRefGoogle Scholar
Hutchinson, J. and Tonegawa, Y.. Convergence of phase interfaces in the van der Waals–Cahn–Hilliard theory. Calc. Var. 10 (2000), 4984.Google Scholar
Kohn, R. V. and Sternberg, P.. Local minimisers and singular perturbations, Proceedings of the Royal Society of Edinburgh. A, 111, pp. 69–84.Google Scholar
Leonardi, G. P.. Infiltrations in immiscible fluids systems. Proc. Roy. Soc. Edinb. A 131 (2007), 0.Google Scholar
Modica, L.. Gradient theory of phase transitions with boundary contact energy. Ann. Inst. Henri Poincare (C) Anal. Non Lineaire 4 (1987), 487512.Google Scholar
Modica, L. and Mortola, S.. Un esempio di $\Gamma$–convergenza. Boll. Un. Mat. Ital. B(5), 14 (1977), 285299.Google Scholar
Modica, L.. $\Gamma$-convergence to minimal surfaces problem and global solutions of $\Delta u = 2(u^3 - u)$, Proceedings of the International Meeting on Recent Methods in Nonlinear Analysis (Rome, 1978, Pitagora, Bologna, 1979), pp. 223–244.Google Scholar
Morgan, F.. $(M,\, \varepsilon,\, \delta )$–Minimal Curve Regularity. Proc. Am. Math. Soc. 120 (1994), 677686.Google Scholar
Morgan, F.. Soap bubbles in $\mathbb {R}^2$ and in surfaces. Pac. J. Math. 165 (1994), 347361.Google Scholar
Owen, N. C., Rubinstein, J. and Sternberg, P.. Minimizers and gradient flows for singular perturbated bi-stable potentials with a Dirichlet condition. Proc. R. Soc. Lond. A 429 (1990), 505532.Google Scholar
Savin, O. and Valdinoci, E.. $\Gamma$-convergence for nonlocal phase transitions. Ann. Inst. Henri Poincare (C) Anal. Non Lineaire 29 (2012), 479500.CrossRefGoogle Scholar
Sternberg, P.. The effect of a singular perturbation on nonconvex variational problems. Arch. Rational Mech. Anal. 101 (1988), 209260.Google Scholar
Taylor, J. E.. The structure of singularities in soap–bubble–like and soap–film–like minimal surfaces. Ann. Math. 103 (1976), 489539.Google Scholar
White, B.. Existence of least energy configurations of immiscible fluids. Geom. Anal. 6 (1996), 151161.Google Scholar
White, B.. Lecture notes on flat chains and geometric measure theory, Stanford 2012 (notes by Otis Chodosh).Google Scholar
Zuniga, A. and Sternberg, P.. On the heteroclinic connection problem for multi-well gradient systems. J. Differ. Equ. 261 (2016), 39874007.Google Scholar
Figure 0

Figure 1. In the left we show a triod with angles $\theta_1, \theta_2, \theta_3$. In the right there is the corresponding triangle with supplementary angles $\hat{\theta}_1, \hat{\theta}_2, \hat{\theta}_3$ that satisfy the Young's law.

Figure 1

Figure 2. The geometric problem subject to such boundary conditions does not admit a minimum. However, the limiting functional admits a minimizer that forms a boundary layer.

Figure 2

Figure 3. Here is an example of a minimizer that we obtain in theorem 1.3.

Figure 3

Figure 4. The singular set of the minimizer obtained in corollary 5.7 is consisted of three radii of the ball.