On the effects of discontinuous capillarities for immiscible two-phase flows in porous media made of several rock-types

1.   Introduction

The aim of this work is to derive von-Kármán plate theory from nonlinear, three-dimensional, atomistic models in a certain energy scaling as the interatomic distance $\varepsilon$ and the thickness of the material $h$ both tend to zero.

The passage from atomistic interaction models to continuum mechanics (i.e., the limit $\varepsilon \to 0$) has been an active area of research over the last years. In particular, this limit has been well studied for three-dimensional elasticity, cf., e.g., [3,1,19,6,10,17,5,4]. At the same time, there have emerged rigorous results deriving effective thin film theories from three-dimensional nonlinear (continuum) elasticity in the limit of vanishing aspect ratio (i.e., the limit $h \to 0$), cf. [15,14,13,7,16]. First efforts to combine these passages and investigate the simultaneous limits $\varepsilon \to 0$ and $h\to0$ were made in [11,20,21] for membranes (whose energy scales as the thickness $h$) and in [18] for Kirchhoff plates (whose energy scales like $h^3$). In particular, this left open the derivation of the von-Kármán plate theory, which describes plates subject to small deflections with energy scale $h^5$ and might even be the most widely used model for thin structures in engineering. Though we do want to mention [2] for a result regarding discrete von-Kármán plate theory that is motivated numerically and not physically.

Our first aim is to close this gap. For thin films consisting of many atomic layers one expects the scales $\varepsilon$ and $h$ to separate so that the limit $\varepsilon,h\to 0$ along $\frac{h}{ \varepsilon} \to \infty$ is equivalent to first passing to the continuum limit $\varepsilon\to 0$ and reducing the dimension from 3d to 2d in the limit $h\to0$. We will show in Theorem 2.1a) that this is indeed true.

By way of contrast, for ultrathin films consisting of only a few atomic layers, more precisely, if $\varepsilon,h\to 0$ such that the number of layers $\nu = \frac{h}{ \varepsilon}+1$ remains bounded, the classical von-Kármán theory turns out to capture the energy only to leading order in $\frac{1}{\nu}$. The next aim is thus to derive a new finite layer version of the von-Kármán plate theory featuring additional explicit correction terms, see Theorem 2.1b). In view of the fabrication of extremely thin layers, such an analysis might be of some interest also in engineering applications. An interesting question related to such applications, which we do not address here, would be to extend our analysis to heterogeneous structures as in [9,8].

Our third aim concerns a more fundamental modelling point of view which is based on the very low energy of the von-Kármán scaling: If the the plate is not too thick (more precisely, if $\frac{h^5}{ \varepsilon^3} \to 0$), we strengthen the previous results to allow for a much wider range of interaction models, that allow for much more physically realistic atomic interactions (compared to [14,13]) as they can now be invariant under reflections and no longer need to satisfy growth assumptions at infinity, see Theorems 2.2 and 2.3. In particular, this includes Lennard-Jones-type interaction models, see Example 3.

Finally, on a technical note, the proof of the our main result set forth in Section 4 elucidates the appearance and structure of the correction terms in the ultrathin film regime. Both in [18] and the present contribution, at the core of the proof lies the identification of the limiting strain, which in the discrete setting can be seen as a $3 \times 8$ matrix rather than a $3 \times 3$ matrix. In [18] this has been accomplished with the help of ad hoc techniques that allowed to compare adjacent lattice unit cells. Now, for the proof of Proposition 4 we introduce a more general and flexible scheme to capture discreteness effects by splitting the deformation of a typical lattice unit cell into affine and non-affine contributions and passing to weak limits of tailor-made finite difference operators. While for $h \gg \varepsilon$ these operators will tend to a differential operator in the limit, if $h \sim \varepsilon$, finite differences in the $x_3$ direction will not become infinitesimal and lead to lower order corrections in $\frac{1}{\nu}$.

This work is organized as follows: In Section 2, we first describe the atomistic interaction model and then present our results. Our main theorem, Theorem 2.1, details the $\Gamma$-limits for both the thin ($\nu \to \infty$) and ultrathin ($\nu$ bounded) case. Theorems 2.2 and 2.3 then extend these results to more general and more physically realistic models. Section 3 contains a few technical tools to circumvent rigidity problems at the boundary and to compare continuous with discrete quantities. Using these tools we then prove our results in Section 4.

2.   Models and results

2.1.   Atomic model

Let $S \subset \mathbb{R}^2 = \mathbb{R}^2 \times \{ 0 \} \subset \mathbb{R}^3$ be an open, bounded, connected, nonempty set with Lipschitz boundary. To keep the notation simple we will only consider the cubic lattice. Let $\varepsilon>0$ be a small parameter describing the interatomic distance, then we consider the lattice $\varepsilon \mathbb{Z}^3$. We denote the number of atom layers in the film by $\nu \in \mathbb{N}$, $\nu \geq 2$ and the thickness of the film by $h = (\nu-1) \varepsilon$. In the following let us consider sequences $h_n, \varepsilon_n, \nu_n$, $n \in \mathbb{N}$, such that $\varepsilon_n, h_n \to 0$. The macroscopic reference region is $\Omega_n = S \times (0,h_n)$ and so the (reference) atoms of the film are $\Lambda_n = \overline{\Omega}_n \cap \varepsilon_n \mathbb{Z}^3$. We will assume that the energy can be written as a sum of cell energies.

More precisely, as in [18] we let $z^1, \dots, z^8$ be the corners of the unit cube centered at $0$ and write

Furthermore, by $\Lambda_n' = \big( \bigcup_{x \in \Lambda_n} (x + \varepsilon_n \{ z^1, \ldots, z^8 \}) \big) \cap \big( \mathbb{R}^2 \times (0, h_n) \big)$ we denote the set of midpoints of lattice cells $x + [- \varepsilon_n /2, \varepsilon_n /2]^3$ contained in $\mathbb{R}^2 \times [0, h_n]$ for which at least one corner lies in $\Lambda_n$. Additionally, let $\vec{w}(x) = \frac{1}{ \varepsilon_n}(w(x+ \varepsilon_n z^1), \dots, w(x+ \varepsilon_n z^8)) \in \mathbb{R}^{3 \times 8}$. Then, we assume that the atomic interaction energy for a deformation map $w \colon \Lambda_n \to \mathbb{R}^3$ can be written as

where $W(x,\cdot) : \mathbb{R}^{3\times 8} \to [0,\infty)$ only depends on those $\vec{w}_i$ with $x+ \varepsilon_n z^i \in \Lambda_n$, which makes (1) meaningful even though $w$ is only defined on $\Lambda_n$.

As a full interaction model with long-range interaction would be significantly more complicated in terms of notation and would result in a much more complicated limit for finitely many layers, we restrict ourselves to these cell energies.

In the following we will sometimes discuss the upper and lower part of a cell separately. We write $A = (A^{(1)}, A^{(2)})$ with $A^{(1)}, A^{(2)} \in \mathbb{R}^{3 \times 4}$ for a $3 \times 8$ matrix $A$.

If the full cell is occupied by atoms, i.e., $x+ \varepsilon_n z^i \in \Lambda_n$ for all $i$, then we assume that $W$ is is given by a homogeneous cell energy $W_{\rm cell}: \mathbb{R}^{3\times 8} \to [0,\infty)$ with the addition of a homogeneous surface energy $W_{\rm surf}: \mathbb{R}^{3\times 4} \to [0,\infty)$ at the top and bottom. That means,

Example 1. A basic example is given by a mass-spring model with nearest and next to nearest neighbor interaction:

$E_{\rm atom}$ can be written in the form (1) by setting

and

We will also allow for energy contributions from body forces $f_n \colon \Lambda_n \to \mathbb{R}^3$ given by

We will assume that the $f_n$ do not depend on $x_3$, that $f_n(x) = 0$ for $x$ in an atomistic neighborhood of the lateral boundary, see (17), and that there is no net force or first moment,

to not give a preference to any specific rigid motion. At last, we assume that after extension to functions $\bar{f}_n$ which are piecewise constant on each $x + (-\frac{ \varepsilon_n}{2}, \frac{ \varepsilon_n}{2})^2$, $x \in \varepsilon_n \mathbb{Z}^2$, $h_n^{-3} \bar{f}_n \to f$ in $L^2(S)$.

Overall, the energy is given as the sum

Due to the factor $\frac{ \varepsilon_n^3}{h_n}$ this behaves like an energy per unit (undeformed) surface area.

Let us make some additional assumptions on the interaction energy. We assume that $W_{\rm cell}$, $W_{\rm surf}$, and all $W(x, \cdot)$ are invariant under translations and rotations, i.e., they satisfy

for any $A \in \mathbb{R}^{3 \times 8}$ or $A \in \mathbb{R}^{3 \times 4}$, respectively, and any $c \in \mathbb{R}^3$ and $R \in {\rm{SO}}(3)$. Furthermore, we assume that $W_{\rm cell}(Z) = W(x, Z) = 0$, which in particular implies $W_{\rm surf}(Z^{(1)}) = W_{\rm surf}(Z^{(2)}) = 0$, where $(Z^{(1)}, Z^{(2)}) = Z$. At last we assume that $W$ and $W_{\rm cell}$ are $C^2$ in a neighborhood of $Z$, while $W_{\rm surf}$ is $C^2$ in neighborhood of $Z^{(1)}$.

Since our model is translationally invariant, it is then equivalent to consider the discrete gradient

with

instead of $\vec{w}(x)$ for any $x$ with $x+ \varepsilon_n z^i \in \Lambda_n$ for all $i$. In particular, the discrete gradient satisfies

The bulk term is also assumed to satisfy the following single well growth condition.

(G) Assume that there is a $c_0>0$ such that

for all $A \in \mathbb{R}^{3 \times 8}$ with $\sum_{i = 1}^8 A_{\cdot i} = 0$.

2.2.   Rescaling and convergence of displacements

In the same way as in a pure continuum approach, it is convenient to rescale the reference sets to the fixed domain $\Omega = S \times (0,1)$. For $x \in \mathbb{R}^3$ let us always write $x = (x', x_3)^T$ with $x'\in \mathbb{R}^2$. We define $\tilde{\Lambda}_n = H_n^{-1} \Lambda_n$ and $\tilde{\Lambda}_n' = H_n^{-1} \Lambda_n'$ with the rescaling matrix

A deformation $w \colon \Lambda_n \to \mathbb{R}^3$ can be identified with the rescaled deformation $y \colon \tilde{\Lambda}_n \to \mathbb{R}^3$ given by $y(x) = w(H_n x)$. We then write $E_n(y)$ for $E_n(w)$. The rescaled discrete gradient is then given by

for $x \in \tilde{\Lambda}_n'$, where now

For a differentiable $v \colon \Omega \to \mathbb{R}^k$ we analogously set $\nabla_n v : = \nabla v H_n^{-1} = (\nabla'v, \frac{1}{h_n} \partial_3 v)$.

In Section 3 we will discuss a suitable interpolation scheme with additional modifications at $\partial S$ to arrive at a $\tilde{\tilde{y}}_{n} \in W^{1,2}(\Omega; \mathbb{R}^3)$ corresponding to $y_n$. Furthermore, for sequences in the von-Kármán energy scaling we will expect $y_n$ and $\tilde{\tilde{y}}_{n}$ to be close to a rigid motion $x \mapsto R^*_n(x+c_n)$ for some $R^*_n, c_n$ and will therefore be interested in the normalized deformation

which would then be close to the identity. The von-Kármán displacements in the limit will then be found as the limit objects of

2.3.   The $\Gamma$-convergence result

To describe the limit energy, let $Q_{\rm cell}(A) = D^2W_{\rm cell}(Z)[A,A]$ for $A\in \mathbb{R}^{3 \times 8}$ and $Q_{\rm surf}(A) = D^2W_{\rm surf}(Z^{(1)})[A,A]$ for $A \in \mathbb{R}^{3 \times 4}$. As the reference configuration is stress free, frame indifference implies

for all $A \in \mathbb{R}^{3 \times 8}$, $A' \in \mathbb{R}^{3 \times 4}$ and all skew symmetric $B \in \mathbb{R}^{3 \times 3}$. As in continuum elasticity theory this just follows from looking at $0 = \frac{\partial}{\partial t}\frac{\partial}{\partial A} W_{\rm cell}(e^{tB}A) |_{t = 0,A = Z}$.

In particular,

for all $c \in \mathbb{R}^3$ and all skew symmetric $B \in \mathbb{R}^{3 \times 3}$.

We introduce a relaxed quadratic form on $\mathbb{R}^{3 \times 8}$ by

By Assumption (G) $Q_{\rm cell}$ is positive definite on $( \mathbb{R}^3 \otimes e_3) Z$. Therefore, for each $A \in \mathbb{R}^{3 \times 8}$ there exists a (unique) $b = b(A)$ such that

Here we used (7) to arrive at the symmetric version. Furthermore, the mapping $A \mapsto b(A)$ is linear. (If $((v_i \otimes e_3) Z)_{i = 1,2,3}$ is a $Q_{\rm cell}$-orthonormal basis of $( \mathbb{R}^3 \otimes e_3) Z$, then $b(A) = - \sum_{i = 1}^3 Q_{\rm cell}[(v_i \otimes e_3) Z, A]$, where $Q_{\rm cell}[\cdot, \cdot]$ denotes the symmetric bilinear form corresponding to the quadratic form $Q_{\rm cell}(\cdot)$.)

At last, let us write

for any $A \in \mathbb{R}^{2 \times 2}$.

We are now in place to state our main theorem in its first version.

Theorem 2.1. (a) If $\nu_n \to \infty$, then $\frac{1}{h_n^4} E_n\stackrel{\Gamma}{\longrightarrow} E_{\rm vK}$ with

where $G_1(x') = {\rm sym} \nabla' u(x') + \tfrac{1}{2} \nabla' v(x') \otimes \nabla' v(x')$ and $G_2(x') = - (\nabla')^2 v(x')$. More precisely, for every sequence $y_n$ with bounded energy $\frac{1}{h_n^4} E_n(y_n) \leq C$, there exists a subsequence (not relabeled), a choice of $R^*_n \in {\rm{SO}}(3), c_n \in \mathbb{R}^3$, and maps $u\in W^{1,2}(S; \mathbb{R}^2)$, $v \in W^{2,2}(S)$ such that $(u_n, v_n)$ given by (5), (6) and (4) satisfy $u_n \rightharpoonup u$ in $W^{1,2}_{\rm loc}(S; \mathbb{R}^2)$, $v_n \to v$ in $W^{1,2}_{\rm loc}(S)$, $R^*_n \to R^*$, and

On the other hand, this lower bound is sharp, as for every $u\in W^{1,2}(S; \mathbb{R}^2)$, $v \in W^{2,2}(S)$, and $R^* \in {\rm{SO}}(3)$ there is a sequence $y_n$ such that $u_n \rightharpoonup u$ in $W^{1,2}_{\rm loc}(S; \mathbb{R}^2)$, $v_n \to v$ in $W^{1,2}_{\rm loc}(S)$ (where we can take $R^*_n = R^*$, $c_n = 0$ without loss of generality) and

(b) If $\nu_n \equiv \nu \in \mathbb{N}$, then $\frac{1}{h_n^4} E_n\stackrel{\Gamma}{\longrightarrow} E^{(\nu)}_{\rm vK}$, to be understood in exactly the same way as in a), where

Here,

In the following we use the notation $E_{\rm vK}(u,v)$, respectively, $E^{(\nu)}_{\rm vK}(u,v)$, for the functionals without the force term.

Example 2. Theorem 2.1 applies to the interaction energy of Example 1 if $W_{\rm cell}$ is augmented by an additional penalty term $+ \chi(\vec{w})$ which vanishes in a neighborhood of ${\rm{SO}}(3) Z$ but is $\ge c > 0$ in a neighborhood of ${\rm O}(3)Z \setminus {\rm{SO}}(3) Z$, so as to guarantee orientation preservation.

Remark 1. 1. The result in a) is precisely the functional one obtains by first applying the Cauchy-Born rule (in 3d) in order to pass from the discrete set-up to a continuum model and afterwards computing the (purely continuum) $\Gamma$-limit on the energy scale $h^5$ as $h\to0$ as in [13]. Indeed, the Cauchy-Born rule associates the continuum energy density

to the atomic interaction $W_{\rm cell}$, and so $Q_{\rm cell}(AZ) = D^2W_{\rm CB}(Z)[A,A] = : Q_{\rm CB}(A)$ for $A\in \mathbb{R}^{3 \times 3}$, in particular,

2. In contrast, for finite $\nu$ non-affine lattice cell deformations of the form $A Z_- + a M$, $A \in \mathbb{R}^{3 \times 3}$, $a \in \mathbb{R}$ need to be taken into account. While $A Z_-$ is non-affine in the out-of-plane direction, $aM$ distorts a lattice unit cell in-plane in a non-affine way.

3. Suppose that in addition $W_{\rm cell}$ and $W_{\rm surf}$ satisfy the following antiplane symmetry condition:

where $P$ is the reflection $P(x',x_3) = (x',-x_3)$. This holds true, e.g., in mass-spring models such as in Example 1. As both terms in $G_3$ switch sign under this transformation, while the affine terms with $G_1$ and $G_2$ remain unchanged, one finds that the quadratic terms in $E^{(\nu)}_{\rm vK}$ decouple in this case and we have

4. Standard arguments in the theory of $\Gamma$-convergence show that for a sequence $(y_n)$ of almost minimizers of $E_n$ the in-plane displacement $u_n$, the out-of-plane displacement $v_n$ and the overall rotation $R^*_n$ converge (up to subsequences) to a minimizer $(u, v, R^*)$ of $E_{\rm vK}$, respectively, $E_{\rm vK}^{(\nu)}$.

5. For the original sequence $y_n$ near the lateral boundary there can be lattice cells for which only a subset of their corners belong to $\Lambda_n$. As a consequence these deformation cannot be guaranteed to be rigid on such cells and the scaled in-plane and out-of-plane displacements may blow up. We thus chose to modify $y_n$ in an atomistic neighborhood of the lateral boundary so as to pass to the globally well behaved quantities $\tilde{y}_n$, see Section 3. For the original sequence $y_n$, Theorem 2.1 implies a $\Gamma$-convergence result with respect to weak convergence in $W^{1,2}_{\rm loc}$.

2.4.   The $\Gamma$-convergence result under weaker assumptions

One physically unsatisfying aspect of Theorem 2.1 is the strong growth assumption (G) which is in line with the corresponding continuum results [13]. The problem is actually two-fold. First, typical physical interaction potentials, like Lennard-Jones potentials, do not grow at infinity but converge to a constant with derivatives going to $0$. And second, (G) also implies that $W_{\rm cell}(-Z)>W_{\rm cell}(Z)$. In particular, the atomistic interaction could not even be ${\rm{O}}(3)$-invariant.

Contrary to the continuum case, it is actually possible to remove these restrictions in our atomistic approach. Indeed, if one assumes $\nu_n^5 \varepsilon_n^2 \to 0$ or equivalently $h_n^5/ \varepsilon_n^3 \to 0$, then the von-Kármán energy scaling implies that the cell energy at every single cell must be small. In terms of the number of atom layers $\nu$, this condition includes the case of fixed $\nu$, as well as the case $\nu_n \to \infty$ as long as this divergence is sufficiently slow, namely $\nu_n \ll \varepsilon_n^{-2/5}$.

In this case, growth assumptions at infinity should no longer be relevant. In fact, we can replace (G) by the following much weaker assumption with no growth at infinity and full ${\rm{O}}(3)$-invariance.

(NG) Assume that $W_{\rm cell}(A) = W_{\rm cell}(-A)$ and that there is some neighborhood $U$ of ${\rm{O}}(3)Z$ and a $c_0>0$ such that

for all $A \in U$ with $\sum_{i = 1}^8 A_{\cdot i} = 0$ and

for all $A \notin U$ with $\sum_{i = 1}^8 A_{\cdot i} = 0$.

One natural problem arising from this is that atoms that are further apart in the reference configuration can end up at the same position after deforming. In particular, due to the full ${\rm{O}}(3)$-symmetry, neighboring cells can be flipped into each other without any cost to the cell energies, which completely destroys any rigidity that one expects in this problem.

As a remedy, whenever we assume (NG), we will add a rather mild non-penetration term to the energy that can be thought of as a minimal term representing interactions between atoms that are further apart in the reference configuration. To make this precise, for small $\delta, \gamma > 0$ let $V \colon \mathbb{R}^3\times \mathbb{R}^3 \to [0,\infty]$ be any function with $V(v,w) \geq \gamma$ if $\lvert v -w \rvert < \delta$ and $V(v,w) = 0$ if $\lvert v -w \rvert \geq 2\delta$. Then define

Then, $\gamma >0$ ensures that there is a positive energy contribution whenever two atoms are closer than $\delta \varepsilon$.

The overall energy is then given by

Theorem 2.2. Assume that $\nu_n^5 \varepsilon_n^2 \to 0$, that $f_n = 0$, that $E_n$ is given by (13), and that (G) is replaced by (NG). Then all the statements of Theorem 2.1 remain true, where now $R^*_n, R^* \in {\rm{O}}(3)$.

Note that in this version, we assume $f_n = 0$. Indeed, if one were to include forces, one can typically reduce the energy by moving an atom infinitely far away in a suitable direction. Without any growth assumption in the interaction energy this can easily lead to $\inf E_n = - \infty$ and a loss of compactness. However, this is just a problem about global energy minimization. Not only should there still be well-behaved local minima of the energy, but the energy barrier in between should become infinite in the von-Kármán energy scaling.

In the spirit of local $\Gamma$-convergence, we can thus consider the set of admissible functions

where ${\Lambda'_n}^{\circ}$ labels 'interior cells' away from the lateral boundary, cf. Section 3. This leads us to the total energy

We then have a version of the $\Gamma$-limit that does allow for forces.

Theorem 2.3. Assume that $\nu_n^5 \varepsilon_n^2 \to 0$, that $E_n$ is given by (14) with $\delta>0$ sufficiently small, and that (G) is replaced by (NG). Then all the statements of Theorem 2.1 remain true. Furthermore, there is an infinite energy barrier in the sense that

Remark 2. 1. For $n$ large enough, the energy barrier implies that minimizers of the restricted energy (13) correspond to local minimizers of the unrestricted energy (3). The results thus implies convergence of local minimizers of (3) in $\mathcal{S}_\delta$.

2. To formulate it differently, if a sequence $(w_n)$ is not separated by a diverging (unrestricted) energy barrier from the reference state id, i.e. each $w_n$ can be connected by a continuous path of deformations $(w_n^t)_{t \in [0,1]}$ with equibounded energy $E_{\rm atom}(w_n^t) + E_{\rm body}(w_n^t)$, then $w_n \in \mathcal{S}_\delta$ for large $n$. This implies convergence of minimizers of the unrestricted energy under the assumption that a diverging energy barrier cannot be overcome.

3. As the energy only has to be prescribed in $\mathcal{S}_\delta$, Theorem 2.3 also describes local minimizers of energy functionals which are invariant under particle relabeling for point configurations which after labeling with their nearest lattice site by $\{w(x) : x \in \Lambda_n\}$ belong to $\mathcal{S}_{\delta}$, where their energy can be written in the form (14).

Example 3. In the setting of Theorems 2.2 and 2.3, Example 2 can be generalized to energies of the form

where $V_1, V_2$ are pair interaction potentials with $V_i(0) = 0$, $V_i$ $C^2$ in a neighborhood of $0$ and $V_i(r) \ge c_0 \min\{ r^2, 1 \}$ for some $c_0 > 0$. (This is satisfied, e.g., for the Lennard-Jones potential $r \mapsto (1+r)^{-12} - 2 (1+r)^{-6} + 1$.) Due to the non-penetration term in (13) no additional penalty terms for orientation preservation are necessary. Most notably, it is not assumed that $V_i(r) \to \infty$ as $r \to \infty$.

3.   Preparations

We first extend a lattice deformation slightly beyond $\Lambda_n$ and in doing so possibly modify it near the lateral boundary $\partial S \times [0,h_n]$ where lattice cells might not be completely contained in $\bar{\Omega}_n$. Then we interpolate so as to obtain continuum deformations to which the continuum theory set forth in [12,13] applies.

For $x \in \Lambda_n'$, with $\Lambda_n'$ as defined at the beginning of Section 2, we set

and also write $Q_n(\xi) = Q_n(x)$ whenever $\xi \in Q_n(x)$.

3.1.   Modification and extension

On a cell that has a corner outside of $\Lambda_n$ there is no analogue to (G) (or (NG)) and hence no control of $\vec{w}(x)$ in terms of $W(x,\vec{w}(x))$. For this reason we modify our discrete deformations $w : \Lambda_n \to \mathbb{R}^3$ near the lateral boundary of $\Omega_n$.

Let $S_n = \{ x \in S : {\rm{dist}}(x, \partial S) > \sqrt{2} \varepsilon_n \}$ and note that, for $\varepsilon_n > 0$ sufficiently small, $S_n$ is connected with a Lipschitz boundary. (This follows from the fact that $\partial S$ can be parameterized with finitely many Lipschitz charts.) If $x \in \Lambda_n'$ is such that $\overline{Q_n(x)} \cap (S_n \times \mathbb{R}) \ne \emptyset$, we call $Q_{n}(x)$ an inner cell and write $x \in {\Lambda_n'}^{\circ}$. The corners of these cells are the interior atom positions $\Lambda_n^{\circ} = {\Lambda_n'}^{\circ} + \varepsilon_n \{z^1, \ldots, z^8\}$ and the part of the specimen made of such inner cells is denoted

Recall the definition of $\Lambda_n'$ from Section 2 and set

The (lateral) boundary cells $Q_{n}(x)$ are those for which

Later we will also use the rescaled versions of these sets which are denoted $\tilde{\Lambda}_n = H_n^{-1} \Lambda_n$, $\tilde{\bar{\Lambda}}_n = H_n^{-1} \bar{\Lambda}_n$, ${\tilde{\Lambda}_n}^{\circ} = H_n^{-1} \Lambda_n^{\circ}$, $\tilde{\Lambda}'_n = H_n^{-1} \Lambda_n'$, $(\tilde{\Lambda}_n')^{\circ} = {\Lambda_n'}^{\circ}$. The rescaled lattice cells are $\tilde{Q}_n(x) = H_n^{-1} Q_n(H_n x)$.

If $w : \Lambda_n \to \mathbb{R}^3$ is a lattice deformation, following [19] we define a modification and extension $w' : \bar{\Lambda}_n \to \mathbb{R}^3$ as follows. First we set $w'(x) = w(x)$ if $x \in \Lambda_n^{\circ}$. Now partition $\partial \Lambda_n'$ into the $8$ sublattices $\partial \Lambda_{n,i}' = \partial \Lambda_n' \cap \varepsilon_n ( z^i + 2 \mathbb{Z}^3 )$. We apply the following extension procedure consecutively for $i = 1, \ldots, 8$:

For every cell $Q = Q_{n}(x)$ with $x \in \partial \Lambda_{n,i}'$ such that there exists a neighboring cell $Q' = Q_{n}(x')$, i.e. sharing a face with $Q$, on the corners of which $w'$ has been defined already, we extend $w'$ to all corners of $Q$ by choosing an extension $w'$ such that ${\rm{dist}}^2({\bar{\nabla} w(x)}, {\rm{SO}}(3)Z)$ is minimal.

As a result of this procedure, $w'$ will be defined on every corner of each cell neighboring an inner cell. Now we repeat this procedure until $w'$ is extended to $\bar{\Lambda}_n$, i.e., to every corner of all inner and boundary cells. Since $S$ is assumed to have a Lipschitz boundary, the number of iterations needed to define $w'$ on all boundary cells is bounded independently of $\varepsilon$.

Our modification scheme guarantees that the rigidity and displacements of boundary cells can be controlled in terms of the displacements, respectively, rigidity of inner cells, see [19,Lemmas 3.2 and 3.4]¹:

¹We apply these lemmas without a Dirichlet part of the boundary, i.e., $\partial \mathcal{L}'_{ \varepsilon}(\Omega)_* = \emptyset$ in the notation of [19]. Note also that there is a typo in the statement of these lemmas. The set $\mathcal{B}_{ \varepsilon}$ should read $\{ \bar{x} \in \mathcal{L}_{ \varepsilon}'(\Omega)^{\circ} \cup \partial \mathcal{L}_{ \varepsilon}'(\Omega)_* : \bar{x} \notin V_{ \varepsilon} \}$, which in our notation (and without Dirichlet part of the boundary) is a subset of ${\Lambda_n'}^{\circ}$.

Lemma 3.1. There exist constants $c, C > 0$ (independent of $n$) such that for any $w : \Lambda_n \to \mathbb{R}^3$ and $R^* \in {\rm{SO}}(3)$

as well as

For the sake of notational simplicity, we will sometimes write $w$ instead of $w'$.

3.2.   Interpolation

Let $w : \bar{\Lambda}_n \to \mathbb{R}^3$ be a (modified and extended) lattice deformation. We introduce two different interpolations: $\tilde{w}$ and $\bar{w}$. $\tilde{w} \in W^{1,2}(\Omega^{\rm out}_n; \mathbb{R}^3)$ is obtained by a specific piecewise affine interpolation scheme as in [18,19] which in particular associates the exact average of atomic positions to the center and to the faces of lattice cells. This will allow for a direct application of the results in [13] on continuum plates. By way of contrast, $\bar{w}$ is a piecewise constant interpolation on the lattice Voronoi cells of $\bar{\Lambda}_n$. The advantage of this interpolation will be that a discrete gradient of $w$ translates into a continuum finite difference operator acting on $\bar{w}$.

Let $x \in \Lambda_n'$. In order to define $\tilde{w}$ on the cube $\overline{Q(x)}$ we first set $\tilde{w}(x) = \frac{1}{8} \sum_{i = 1}^8 w(x + \varepsilon_n z^i)$. Next, for the six centers $v^1, \ldots, v^6$ of the faces $F^1, \ldots, F^6$ of $[-\frac{1}{2}, \frac{1}{2}]^3$ we set $\tilde{w}(x + \varepsilon_n v^i) = \frac{1}{4} \sum_j w(x + \varepsilon_n z^j)$, where the sum runs over those $j$ such that $z^j$ is a corner of the face with center $v^i$. Finally, we interpolate linearly on each of the 24 simplices

with $|z^i - z^j| = 1$, $|z^i - v^k| = |z^j - v^k| = \tfrac{1}{\sqrt{2}}$, i.e., whose corners are given by the cube center and the center and two neighboring vertices of one face. Note that for this interpolation

for every face $x + \varepsilon_n F^k$ of $Q(x)$.

For the second interpolation we first let $V_n^{\rm out} : = \big( \bigcup_{x \in \bar{\Lambda}_n} ( x + [-\tfrac{ \varepsilon_n}{2}, \tfrac{ \varepsilon_n}{2}]^3 ) \big)^{\circ}$ and then define $\bar{w} \in L^2(V_n^{\rm out}; \mathbb{R}^3)$ by $\bar{w}(\xi) = w(x)$ for all $\xi \in x + (-\tfrac{ \varepsilon_n}{2}, \tfrac{ \varepsilon_n}{2})^3$, $x \in \bar{\Lambda}_n$. Note that

with $\langle \bar{w} \rangle = \frac{1}{8} \sum_{i = 1}^8 \bar{w}(x+ \varepsilon_n z^i)$ defines a piecewise constant mapping on $\Omega^{\rm out}_{n}$ such that

It is not hard to see that the original function controls the interpolation and vice versa.

Lemma 3.2. There exist constants $c, C > 0$ such that for any (modified, extended and interpolated) lattice deformation $\tilde{w} : \Omega^{\rm out}_{n} \to \mathbb{R}^3$ and any cell $Q = Q_n(x)$, $x \in \Lambda_n'$,

Proof. After translation and rescaling we may without loss assume that $\varepsilon_n = 1$ and $Q = (0,1)^3$, hence $x = (\tfrac{1}{2}, \tfrac{1}{2}, \tfrac{1}{2})^T$. The claim then is an immediate consequence of the fact that both

are norms on the finite dimensional space of continuous mappings $\tilde{w}$ which are affine on each $\operatorname{co} ( x, x + v^k, x + z^i, \varepsilon_n z^j )$ with $|z^i - z^j| = 1$, $|z^i - v^k| = |z^j - v^k| = \tfrac{1}{\sqrt{2}}$, and which have $\int_Q \tilde{w}(\xi) \, d\xi = 0$.

Lemma 3.3. There exist constants $c, C > 0$ such that for any (modified, extended and interpolated) lattice deformation $\tilde{w} : \Omega^{\rm out}_{n} \to \mathbb{R}^3$ and any cell $Q = Q_n(x)$, $x \in \Lambda_n'$,

This is in fact [19,Lemma 3.6]. We include a simplified proof.

Proof. After translation and rescaling we may without loss assume that $\varepsilon_n = 1$ and $Q = (0,1)^3$. The geometric rigidity result [14,Theorem 3.1] (indeed, an elementary version thereof) yields

By definition also

The claim then follows from applying Lemma 3.2 to $\xi \mapsto \tilde{w}(\xi) - R\xi$ for each $R \in {\rm{SO}}(3)$.

For a sequence $w_n$ of (modified and extended) lattice deformations $w_n : \bar{\Lambda}_n \to \mathbb{R}^3$ with interpolations $\tilde{w}_n : \Omega^{\rm out}_{n} \to \mathbb{R}^3$ and $\bar{w}_n : V_n^{\rm out} \to \mathbb{R}^3$ we consider the rescaled deformations $\tilde{\tilde{{y}}}_n : \tilde{\Omega}^{\rm out}_n \to \mathbb{R}^3$ defined by

and $\bar{ \bar{{y}}}_n : \tilde{V}^{\rm out}_n \to \mathbb{R}^3$ defined by

(Later we will normalize by a rigid change of coordinates to obtain $\tilde{y}_n$ and $\bar{y}_n$.) Their rescaled (discrete) gradients are

for all $x \in\tilde{\Omega}^{\rm out}_{n}$. Finally, the force $f_n$ after extension to $\bar{\Lambda}_n$ is assumed to satisfy

and its the piecewise constant interpolation is $\bar{f}_n : \tilde{V}^{\rm out}_n \to \mathbb{R}^3$.

Remark 3. Suppose $\nu_n = \nu$ constant. We note that for a sequence of mappings $y_n : \Lambda_n \to \mathbb{R}^3$, if $\tilde{\tilde{{y}}}_n \to y$ in $L^2(\Omega; \mathbb{R}^3)$ then $y$ is continuous in $x_3$ and affine in $x_3$ on the intervals $(\frac{i-1}{\nu-1}, \frac{i}{\nu-1})$, $i = 1, \ldots, \nu$. Similarly, if $\bar{\bar{{y}}}_n \to y^*$ in $L^2(S \times (\frac{-1}{2(\nu-1)}, \frac{2\nu -1}{2(\nu-1)}); \mathbb{R}^3)$, then $y^*$ is constant in $x_3$ on the intervals $(\frac{2i-1}{2(\nu-1)}, \frac{2i + 1}{2(\nu-1)})$, $i = 0, \ldots, \nu -1$.

Suppose $y, y^* \in L^2(\Omega; \mathbb{R}^3)$ are piecewise affine, respectively, constant in $x_3$ as detailed above with $y^*(x',x_3) = y(x',\frac{i}{\nu-1})$ if $x_3 \in (\frac{2i-1}{2(\nu-1)}, \frac{2i + 1}{2(\nu-1)})$, $i = 0, \ldots, \nu -1$. It is not hard to see that the following are equivalent.

● $\tilde{\tilde{{y}}} \to y$ in $L^2(\Omega; \mathbb{R}^3)$.

● $\bar{\bar{{y}}} \to y^*$ in $L^2(S \times (\frac{-1}{2(\nu-1)}, \frac{2\nu -1}{2(\nu-1)}); \mathbb{R}^3)$.

● $\frac{ \varepsilon_n^3}{h_n} \sum_{x \in \tilde{\Lambda}_n} | y_n(x) - \mathit{{\rlap{-} \smallint }}_{x + (-\frac{ \varepsilon_n}{2},\frac{ \varepsilon_n}{2})^2 \times (-\frac{ \varepsilon_n}{2h_n},\frac{ \varepsilon_n}{2h_n})} y^*(\xi) \, d\xi|^2 \to 0$.

The same is true in case $\nu_n \to \infty$ for $y = y^*$ if in the second statement $S \times (\frac{-1}{2(\nu-1)}, \frac{2\nu -1}{2(\nu-1)})$ is replaced by $\Omega$.

In particular, limiting deformations do not depend on the interpolation scheme.

4.   Proofs

4.1.   Compactness

For the compactness we will heavily use the corresponding continuum rigidity theorem from [12,Theorem 3] and [13,Theorem 6]:

Theorem 4.1. Let $y \in W^{1,2}(\Omega; \mathbb{R}^3)$ and set $\mathcal{I} = \mathcal{I}(y) = \int_\Omega {\rm{dist}}^2(\nabla_n y, {\rm{SO}}(3))\,dx$. Then there exists maps $R:S \to {\rm{SO}}(3)$ and $\tilde{R} \in W^{1,2}(S; \mathbb{R}^{3 \times 3})$ with $\lvert \tilde{R} \rvert \leq C$, and a constant $R^* \in {\rm{SO}}(3)$ such that

Crucially, none of the constants depend on $n$, $y$, or $\mathcal{I}$.

Furthermore, we will also use the continuum compactness result [12,Lemmas 4 and 5] and [13,Lemma 1,Eq. (96),and Lemma 2] based on the previous rigidity result applied to some sequence $(\hat{y}_n)$.

Theorem 4.2. Let $\hat{y}_n \in W^{1,2}(\Omega; \mathbb{R}^3)$ with $\mathcal{I}(\hat{y}_n)\leq Ch_n^4$. Then there are $R^*_n\in {\rm{SO}}(3)$, $c_n \in \mathbb{R}^3$ as well as a $u \in W^{1,2}(S; \mathbb{R}^2)$ and a $v \in W^{2,2}(S)$ such that $y_n = {R^*_n}^T \hat{y}_n -c_n$ satisfies

And, up to extracting subsequences,

where the upper left $2 \times 2$ submatrix $G''$ of $G$ is given by

with

The following proposition allows us to apply these continuum results.

Proposition 1. In the setting of Theorem 2.1, consider a sequence $w_n$ with

Then,

Here, $\tilde{\tilde{{y}}}_n\in W^{1,2}(\Omega; \mathbb{R}^3)$ is the rescaled, modified, and interpolated version of $w_n$ according to Section 3.

In the setting of Theorem 2.3 the statement remains is true as well, while in the setting of Theorem 2.2 (36) is still true but now $\tilde{ \tilde{{y}}}_n$ is the rescaled, modified, and interpolated version of either $w_n$ or $-w_n$ where the correct sign does depend on $w_n$.

Proof. Rescaling the $w_n$ and applying the modification and interpolation steps from Section 3, we have sequences $\tilde{\tilde{{y}}}_n \in W^{1,2}(\Omega; \mathbb{R}^3)$ and $\bar{\bar{{y}}}_n \in L^{2}(\Omega; \mathbb{R}^3)$. In particular, we can use Theorem 4.1 for this sequence.

Take $R^*_n$ according to Theorem 4.2. Then by Lemmas 3.1 and 3.3,

A standard discrete Poincaré-inequality then shows

for a suitable $\bar{c}_n \in \mathbb{R}^3$. Now $f_n$ does not depend on $x_3$, vanishes close to $\partial S$ where the modification takes place, and satisfies $\sum_{x \in \Lambda_n} f_n = 0$, as well as $\sum_{x \in \Lambda_n} f_n \otimes x' = 0$. Hence, we see that

Using $\lVert \bar{f}_n \rVert_{L^2(S)} \leq C h_n^3$ and abbreviating $\mathcal{I}(\tilde{\tilde{{y}}}_n) = \mathcal{I}_n$, we thus find

On the other hand, due to ${\bf (G)}$ and Lemmas 3.1 and 3.3 we have

Hence,

We thus have

All these statements remain true in the setting of Theorem 2.3 as the Assumptions ${\bf (G)}$ and ${\bf (NG)}$ are equivalent on $\mathcal{S}_\delta$.

Now, consider the setting of Theorem 2.2 with Assumption ${\bf (NG)}$ instead of ${\bf (G)}$, as well as $f_n = 0$ and $\nu_n^5 \varepsilon_n^2 \to 0$ with the energy given by (13). Using (35), we find

for every $x \in {\Lambda'_n}^{\circ}$ and

for all $\bar{x}, \bar{ \bar{{x}} } \in \Lambda_n$. As $\frac{h_n^5}{ \varepsilon_n^3} \to 0$, for $n$ large enough, the right hand side is strictly smaller then $c_0$ or $\gamma$, respectively. Therefore, for all $n$ large enough we have

and

for all $\bar{x}, \bar{ \bar{{x}}} \in \Lambda_n$.

$\bar{\nabla} w_n(x) \in U$ implies $W_{\rm cell}(\bar{\nabla} w_n(x)) \geq c_0\;{\rm{dist}}^2(\bar{\nabla} w_n(x), {\rm{O}}(3)Z)$. In particular, we thus find

Again, for $n$ large enough, this means that every $x\in {\Lambda'_n}^{\circ}$ the discrete gradient $\bar{\nabla} w_n(x)$ is arbitrarily close to ${\rm{O}}(3)Z$ and thus very close to $\sigma_n(x){\rm{SO}}(3)Z$ with a unique $\sigma_n(x) \in \{\pm 1\}$. We now want to show that the sign $\sigma_n(x)$ is the same for all $x$ in the interior cells. As the interior of the union of all these cells is connected, it suffices to show that $\sigma_n$ is the same on any two cells that share a $(d-1)$-face. Indeed, if that were false, we would have some $x,x'$ in cells that share a $(d-1)$-face such that

and

with $Q,Q' \in {\rm{SO}}(3)$. Without loss of generality assume $x = x' + \varepsilon_n e_3$. Then

for all $b\in \mathbb{R}^4$ with $\sum_i b_i = 0$. In particular choosing $b = (-1,+1,+1,-1)$ and $b = (-1,-1,+1,+1)$, we get $\lvert (Q+ Q') e_i \rvert \leq C \frac{h_n^5}{ \varepsilon_n^3}$ for $i = 1,2$. As $Q,Q' \in {\rm{SO}}(3)$, we find $\lvert (Q- Q') e_3 \rvert \leq C \frac{h_n^5}{ \varepsilon_n^3}$. Overall, we see that both deformed cells are almost on top of each other. More specifically,

for $n$ large enough. This is a contradiction to the non-penetration condition (37).

That means, we have

for an $x$-independent $\sigma_n \in \{\pm 1\}$. Applying the modification and interpolation procedure from Section 3 to $\sigma_n w_n$ as in the case (G) above, we find

Now we can directly apply Theorems 4.1 and 4.2 for the continuum objects $\tilde{\tilde{{y}}}_n$. In particular, for $\tilde{y}_n = {R^*_n}^T \tilde{ \tilde{{y}} }_n - c_n$ as defined in (4) and corresponding $u_n$ and $v_n$ as in (5), respectively, (6), after extracting a subsequence from (28) and (29) we get that

For later we also introduce $\bar{y}_n = {R^*_n}^T \bar{ \bar{{y}} }_n - c_n$.

We will also use the following finer statement.

Proposition 2. In the setting of Theorem 4.2, applied to $\tilde{ \tilde{{y}}}_n$ and with $\tilde{y}_n = {R^*_n}^T \tilde{\tilde{{y}}}_n - c_n$, we have

where

Proof. According to Korn's inequality

According to Theorem 4.2, ${\rm{sym}} \nabla' \hat{u}_n$ is bounded in $L^2$ by (23) and (31). Furthermore, $\int {\rm{skew}} \nabla' \hat{u}_n \,dx = 0$ by (27), and $\int \hat{u}_n \,dx$ is bounded due to (28). As

$i = 1,2$, this term is bounded in $L^2$ as well. This shows compactness. To identify the limit and thus show convergence of the entire sequence, note that

by (28) and

for $i = 1,2$ by (30).

(26) and (29) in Theorem 4.2 also show that $\hat{v}_n$ is bounded in $W^{1,2}(\Omega)$ with $\frac{\partial \hat{v}_n}{\partial x_3} \to 1$ and

As a first consequence, we will now describe the limiting behavior of the force term $E_{\rm body}(w_n) = E_{\rm body}(y_n) = \sum_{x \in \tilde{\bar{\Lambda}}_n} f_n(x) \cdot y_n(x)$, where $f_n(x) = f_n(x')$ satisfies (2), (17) and $h_n^{-3} \bar{f}_n \to f$ in $L^2(S)$.

Note that the forces considered are a bit more general than in [13].

Proposition 3. Let $y_n$ be a sequence with $E_n(y_n) \le C h_n^4$ and suppose that (38) holds true for $\tilde{y}_n$, $u_n$, $v_n$ as defined in (4), (5), (6). Assume that $R^*_n \to R^*$. Then

as $n \to \infty$.

Proof. In terms of the extended and interpolated force density we have

By Proposition 2, $h_n^{-1} \Big( \tilde{y}_n - \begin{pmatrix} x' \\ 0 \end{pmatrix} \Big) \to \hat{v} e_3$ in $L^2(\Omega; \mathbb{R}^3)$ with $\hat{v}$ as in (40) and so Remark 3 shows that

if $\nu_n \to \infty$, where in the last step we have used that (2) together with $f_n(x) = f_n(x')$ also implies that $\sum_{x \in \Lambda_n} x_3 f_n(x) = 0$. If $\nu_n = \nu$ constant, then Remark 3 gives

with an analogous argument for the last step.

4.2.   Lower bounds

To show the lower bounds in our $\Gamma$-convergence results, we have to understand the limit of the discrete strain. Let $(y_n)$ satisfy $E_n(y_n) \leq C h_n^4$ and set

By Proposition 1 $(\tilde{ \tilde{{y}} }_n)$ satisfies the assumptions of Theorem 4.2 and Proposition 2 so that, after a rigid change of coordinates, $\tilde{y}_n$ satisfies (23)–(34) and (39)–(42). In particular, by (32) we know that for a subsequence the continuum strain converges as

where $G$ satisfies (33) and (34).

For the discussion of discrete strains, recall that we defined

We define a projection $P$ acting on maps via

in case $\nu_n \equiv \nu < \infty$ and $P = {\rm{id}}$ in case $\nu_n \to \infty$.

Proposition 4. Let $(y_n)_n$ satisfy $E_n(y_n) \leq C h_n^4$ with $\frac{1}{h_n^2} (R_n^T \nabla_n \tilde{y}_n - {\rm{Id}}) \rightharpoonup G$ in $L^2(\Omega; \mathbb{R}^{3 \times 3})$. Then,

in $L^2(\Omega; \mathbb{R}^{3 \times 8})$, where $G_3$ is as in Theorem 2.1.

Proof. The compactness follows from Theorem 4.2. On a subsequence (not relabeled) we thus find $\bar{G}_n \rightharpoonup \bar{G}$. As $R_n \to {\rm{Id}}$ in $L^2$ while being uniformly bounded, we also find

We have

weakly in $L^2(\Omega; \mathbb{R}^{3 \times 3})$ where $G$ satisfies (33) and (34).

In order to discuss the discrete strains in more detail, we separate affine and non-affine contributions. We say that a $b \in \mathbb{R}^8$ is affine if it is an element of the linear span of $b^0, b^1, b^2, b^3$, where $b^0 = (1, \ldots, 1)$ and $b^i = Z^T e_i$, $i = 1,2,3$. Any $b \in \mathbb{R}^8$ which is perpendicular to all affine vectors is called non-affine. I.e., a non-affine $b$ is characterized by $\sum_{i = 1}^8 b_i = 0$ and $Zb = 0$.

We begin by identifying the easier to handle affine part of the limiting strain. By construction we have $R_n \bar{G}_n b^0 \equiv 0$ and so $\bar{G} b^0 = 0 = G Z b^0$. For $i \in \{1, 2, 3\}$ we use that on any $\tilde{Q}_n(x)$, $x \in \tilde{\Lambda}'_n$,

where $y_i = \tilde{y}_n(x' + \varepsilon_n (z^{i})', x_3 + \tfrac{ \varepsilon_n}{h_n} z^i_3)$. So, using (16) for $\tilde{y}_n$,

Analogous arguments yield

By $P_n$ we denote the projection which maps functions to piecewise constant functions via $P_n f(x) = \mathit{{\rlap{-} \smallint }}_{\tilde{Q}_n(x)} f(\xi) \, d\xi$ on $\tilde{Q}_n(x)$. Then $P_n [ R_n\bar{G}_n ] \rightharpoonup \bar{G}$. On the other hand, observing that $Z Z^T = 2 {\rm{Id}}_{3\times 3}$, we find

and

In summary we get that for every affine $b \in \mathbb{R}^8$

For the discussion of the non-affine part of the strain we fix a non-affine $b \in \mathbb{R}^8$, i.e., a $b$ satisfying $\sum_{i = 1}^8 b_i = 0$, $Zb = 0$, and write $b^T = ((b^{(1)})^T, (b^{(2)})^T)$, where $b^{(1)}, b^{(2)} \in \mathbb{R}^4$. Let $Z^{\rm 2dim} : = ((z^1)', (z^2)', (z^3)', (z^4)') \in \mathbb{R}^{2 \times 4}$ be the matrix of two-dimensional directions. Then $Z^{\rm 2dim} (b^{(1)} + b^{(2)}) = 0$ and $\sum_{i = 1}^4 b^{(1)}_i = \sum_{i = 1}^4 b^{(2)}_i = 0$. We introduce the difference operator

The idea is now to separate differences into in-plane and out-of-plane differences, as all in-plane differences are infinitesimal, while out-of-plane differences stay non-trivial if $\nu_n \equiv \nu$ and have to be treated more carefully.

Using

we find

where we have used that $\sum_{i = 1}^4 b^{(1)}_i = \sum_{i = 1}^4 b^{(2)}_i = 0$.

First consider the term (45). Since $\bar{\nabla}_n^{\rm 2dim} \bar{{\rm{id}}}(x - \frac{ \varepsilon_n}{2h_n}e_3) = Z^{\rm 2dim}$ and $Z^{\rm 2dim} (b^{(1)}+b^{(2)}) = 0$, for any $\varphi \in C_c^\infty(\Omega)$ and $i = 1,2$ by (39) and Remark 3 we have

where, either $\tilde{x} = x$ (if $\nu_n \to \infty$), or $\tilde{x} = (x', \frac{\lfloor(\nu-1)x_3\rfloor}{\nu-1})$ (if $\nu_n = \nu$ is constant).

For the third component, we instead have

Now,

uniformly. Therefore, (40) gives

if $\nu_n \to \infty$. For $\nu_n = \nu$ constant however, using (40) and (42) we find

where we have used that $\sum_{i = 1}^8 b_i = 0$.

We still need to find the limit of (44). For any test function $\varphi \in C_c^\infty(\Omega; \mathbb{R}^3)$ we find

Here the penultimate step is true by our specific choice of interpolation to define $\tilde{y}_n$, whereas the last step follows from (30) and $\bar{\nabla}_n^{\rm 2dim} \frac{1}{h_n} e_3 = 0$. If $\nu_n \to \infty$ this converges to $0$. In case $\nu_n = \nu$ constant we obtain from (30)

Summarizing (46), (47), (48), and (49), we see that for non-affine $b$ we have $\bar{G} b = 0$ in case $\nu_n \to \infty$ and

as $\sum_{j = 1}^8 b_j = 0$, if $\nu_n \equiv \nu$.

Elementary computations show that for the affine basis vectors $b^k$, $k \in \{0,1,2,3\}$,

and also

Thus combining with (43), for every $b \in \mathbb{R}^8$ we get

if $\nu_n \to \infty$ and

if $\nu_n = \nu$ is constant. So $\bar{G} = GZ$ if $\nu_n \to \infty$ and

with $Z_-$ as in (12) if $\nu_n = \nu$ is constant. Noting that

with $M$ as in (11) this can be written as

Last, we note that subsequences were indeed not necessary, as the limit is characterized uniquely.

Having established convergence of the strain, the $\liminf$ inequality in Theorems 2.1, 2.2 and 2.3 can now be shown by a careful Taylor expansion of $W(x, \cdot)$, cf. [14,13,18].

Proof of the $\liminf$ inequality in Theorems 2.1, 2.2 and 2.3. The $\liminf$ inequality in Theorem 2.3 is an immediate consequence of the $\liminf$ inequality in Theorem 2.1 applied to a cell energy $W_{\rm cell}'$ of the form

Furthermore, in view of Proposition 3 it suffices to establish the lower bound for $f_n = 0$.

Assume that $(y_n)$ is a sequence of atomistic deformations such that

so that by Proposition 1 its modification and interpolation $(\tilde{y}_n)$ verifies the assertions of Theorem 4.2. Set

By frame indifference and nonnegativity of the cell energy we have

First assume that $\nu_n \to \infty$ as $n \to \infty$. Due to nonnegativity of $W_{\rm surf}$ we can estimate

where $\chi_n$ is the characteristic function of $\{ x \in \Omega^{\rm in}_n : \bar{G} \le h_n^{-1} \} \subset \Omega$ and

so that $t^{-2} \omega(t) \to 0$ as $t \to 0$. Since $\bar{G}_n^2$ is bounded in $L^1(\Omega; \mathbb{R}^{3 \times 8})$ and

uniformly,

Moreover, $\chi_n \to 1$ boundedly in measure and so by Proposition 4 $\chi_{n} \bar{G}_n \rightharpoonup \bar{G} = GZ$, where $G$ satisfies (33) and (34). By lower semicontinuity it follows that

Integrating the last expression over $x_3 \in (0,1)$ and noting that the integral of the cross terms vanish we obtain

Now suppose that $\nu_n \equiv \nu \in \mathbb{N}$. We let $\chi_n$ as above but now define

so that still $t^{-2} \omega(t) \to 0$ as $t \to 0$. With $\bar{G}(x) = (\bar{G}^{(1)}(x), \bar{G}^{(2)}(x))$ we have

where we have used that $\bar{G}$ is constant on $S \times (0,\frac{1}{\nu-1})$ and on $S \times (\frac{\nu-2}{\nu-1},1)$. Here (see Eq. (10) for $G_3$),

The bulk part is estimated as

where we have used that $\sum_{k = 1}^{\nu-1}\tfrac{(2k - \nu)^2}{(2\nu-2)^3} = \frac{\nu(\nu-2)}{24(\nu-1)^2}$.

For the surface part first note that by (8), for any $A = (a_{ij}) \in \mathbb{R}^{3 \times 3}$ and $B \in \mathbb{R}^{3 \times 4}$ we have

where $a_{\cdot 3}$ denotes the third column, $a_{3 \cdot}$ the third row and $A'' = (a_{ij})_{1 \le i,j \le 2}$ the upper left $2 \times 2$ part of $A$. Thus also

It follows that

and so

Adding bulk and surface contributions and integrating over $x'$ we arrive at

Note that in the Theorem 2.1 the skew symmetric part of $G_1$ is then just set to be zero as it does not impact the energy.

4.3.   Upper bounds

Without loss of generality we assume that $R^* = {\rm{Id}}$. (For general $R^*$ one just considers the sequence $R^* y_n$ with $y_n$ as in (50) and $R^*_n = R^*$ below.

If $u : S \to \mathbb{R}^2$ and $v : S \to \mathbb{R}$ are smooth up to the boundary, we choose a smooth extension to a neighborhood of $S$ and define the lattice deformations $y_n : \tilde{\bar{\Lambda}}_n \to \mathbb{R}^3$ by restricting to $\tilde{\bar{\Lambda}}_n$ the mapping $y_n : \overline{\tilde{\Omega}_n^{\rm out}} \to \mathbb{R}^3$, defined by

for all $x \in \overline{\tilde{\Omega}_n^{\rm out}}$. Here $d : \overline{\tilde{\Omega}_n^{\rm out}} \to \mathbb{R}^3$ will be determined later, see (55) and (56) for films with many, respectively, a bounded number of layers. In both cases, $d$ is smooth and bounded in $W^{1,\infty}(\overline{\tilde{\Omega}_n^{\rm out}}; \mathbb{R}^3)$ uniformly in $n$.

We let $R^*_n = {\rm{Id}}$ and $c_n = 0$ for all $n$ and define $\tilde{y}_n \in W^{1,2}(\tilde{\Omega}_n^{\rm out}; \mathbb{R}^3)$ as in (4) by interpolating as in Section 3 (more precisely, descaling to $w_n$ and then interpolating and rescaling) to obtain $\tilde{y}_n = \tilde{ \tilde{{y}}}_n$. Analogously we let $\bar{y} = \bar{\bar{{y}}}$. We define $u_n$ and $v_n$ as in (5) and (6), respectively. It is straightforward to check that indeed $u_n \to u$ in $W^{1,2}(S; \mathbb{R}^2)$ and $v_n \to v$ in $W^{1,2}(S)$.

In order to estimate the energy of $y_n$ we need to compute its discrete gradient. Instead of directly calculating $\bar{\nabla} \bar{y}_n = (\bar{\partial}_1 \bar{y}_n, \ldots, \bar{\partial}_8 \bar{y}_n)$ it is more convenient to first determine $\bar{D} y_n = (\bar{D}_1 y_n, \ldots, \bar{D}_8 y_n)$ which for each $x \in \tilde{\Lambda}'_n - (\frac{ \varepsilon_n}{2}, \frac{ \varepsilon_n}{2}, \frac{ \varepsilon_n}{2 h_n})$ is defined by

where for $x \in \tilde{\Omega}_n^{\rm out}$ we have set

so that $\tilde{Q}_n(x) = \hat{x} + (0, \varepsilon_n)^2 \times (0, (\nu_n-1))$. We set $a^i = \frac{1}{2}(1,1,1)^T + z^i \in \{0,1\}^3$ and write $A : = (a^1, \ldots, a^8) = Z + \frac{1}{2}(1,1,1)^T \otimes (1, 1, 1, 1, 1, 1, 1, 1)^T$. Note that

In particular, if $\bar{D} y_n(x)$ is affine, i.e., $\bar{D} y_n(x) = F A$ for some $F \in \mathbb{R}^{3 \times 3}$, then

and so $\bar{\nabla} \bar{y}_n(x) = F Z$.

For $x$ in a fixed cell $\tilde{Q}_n(x) = \hat{x} + (0, \varepsilon_n)^2 \times (0, (\nu_n-1))$, Taylor expansion of $y_n$ (restricted to $\overline{\tilde{Q}}_n(x)$) yields

for some $\zeta_{ \varepsilon_n} \in \hat{x} + [0, \varepsilon_n]^2\times[0, \varepsilon_n h_n^{-1}]$. Plugging in (50) we get

It follows that

We define the skew symmetric matrix $B(\hat{x}) = B_n(\hat{x})$ by

where we have written $d' = (d_1, d_2)^T$ for $d = (d_1, d_2, d_3)^T$, and consider the special orthogonal matrix

Now compute

Here, the error term is uniform in $\hat{x}$.

We can now conclude the proof of Theorems 2.1, 2.2 and 2.3.

Proof of the $\limsup$ inequality in Theorems 2.1, 2.2 and 2.3. As the discrete gradient $\bar{\nabla}_n \bar{y}_n$ is uniformly close to ${\rm{SO}}(3) Z$, the following arguments apply to show that $y_n$ defined by (50) serves as a recovery sequence in all three theorems. Moreover, in view of Proposition 3 it suffices to construct recovery sequences for $f_n = 0$.

We first specialize now to the case $\nu_n \to \infty$. For

choosing $d(x) = x_3 d_0(x') + \frac{x_3^2 - x_3}{2} d_1(x')$ with

according to (9), from (52) and (53) we obtain

and, Taylor expanding $W_{\rm cell}$, we see that due to the smoothness of $u$ and $v$ the piecewise constant mappings $x \mapsto h_n^{-4} W_{\rm cell} (\bar{\nabla}\bar{y}_n(x)) = h_n^{-4} W_{\rm cell} (e^{-B(\hat{x})} \bar{\nabla} \bar{y}_n(x))$ converge uniformly to

This shows that

and thus finishes the proof in case $\nu_n \to \infty$.

Now suppose that $\frac{ \varepsilon_n}{h_n} \equiv \frac{1}{\nu-1}$. Abbreviating $(\nabla')^2 v(\hat{x}') = -G_2(\hat{x}') = -G_2 = (f_{ij}) \in \mathbb{R}^{2 \times 2}$, we observe that

and hence, with $b = b(\hat{x}') = \big( (\partial_{11}+\partial_{12}) v(\hat{x}'), (\partial_{21}+\partial_{22}) v(\hat{x}'), 0 \big)^T = (f_{11} + f_{12}, f_{21} + f_{22}, 0)^T$,

This shows that

We define the affine part of the strain $G(x) = G_1(x') + (x_3 - \tfrac{1}{2}) G_2(x')$ as in (54). The non-affine part is abbreviated by $\frac{1}{2(\nu - 1)} G_3(x')$ as in (10). Then using (53) we can write

where we have used (52) and (51).

We set

according to (9) and define $d : S' \times [0,1] \to \mathbb{R}$, $S'$ a neighborhood of $S$, inductively by $d(x,0) = 0$ and

for $j = 1, \ldots, \nu-1$. Then $d$ is smooth in $x'$ and piecewise linear in $x_3$, more precisely, affine in $x_3$ in between two atomic layers: On $S' \times [\frac{j-1}{\nu-1}, \frac{j}{\nu-1}]$, $j \in \{ 1, \ldots, \nu-1 \}$, it satisfies

since $\hat{x}_3 = \hat{x}_3(x) = \tfrac{j-1}{\nu-1}$. Taylor expanding $W_{\rm cell}$, we see that the piecewise constant mappings $x \mapsto h_n^{-4} W_{\rm cell} (\bar{\nabla} \bar{y}_n(x)) = h_n^{-4} W_{\rm cell} (e^{-B(\hat{x})} \bar{\nabla} \bar{y}_n(x))$ converge uniformly on $S' \times [\frac{j-1}{\nu-1}, \frac{j}{\nu-1}]$ to

for each $j \in \{ 1, \ldots, \nu-1 \}$. Since $\tfrac{1}{\nu-1} \sum_{j = 1}^{\nu-1} \tfrac{2j-\nu}{2(\nu-1)} = 0$ and $\tfrac{1}{\nu-1} \sum_{j = 1}^{\nu-1} \big( \tfrac{2j-\nu}{2(\nu-1)} \big)^2 = \tfrac{\nu(\nu-2)}{12(\nu-1)^2}$, this shows

For the surface part we write $\bar{\nabla}\bar{y}_n = ([\bar{\nabla}\bar{y}_n]^{(1)}, [\bar{\nabla}\bar{y}_n]^{(2)})$ and use that the piecewise constant mappings $S' \times [0, \frac{1}{\nu-1}] \to \mathbb{R}$,

converge uniformly to

Similarly, the mappings $S' \times [\frac{\nu-2}{\nu-1}, 1] \to \mathbb{R}$,

converge uniformly to

So with $S^{\rm out}_n$ such that $\tilde{\Omega}^{\rm out}_n = S^{\rm out}_n \times (0,1)$,

Summarizing (58) and (57), we have shown that

as $n \to \infty$, where we have also used that the contribution of the lateral boundary cells $\varepsilon_n^3 h_n^{-5} \sum_{x \in \partial \tilde{\Lambda}_n'} W(x, \bar{\nabla}y_n(x))$ is negligible in the limit $n \to \infty$.

Proof of the energy barrier in Theorem 2.3. If a sequence of $w_n\in \mathcal{S}_{\delta}$ satisfies the energy bound $E_n(w_n) \leq C h_n^4$, then the proof of Proposition 1 shows $\frac{ \varepsilon_n^3}{h_n} E_{\rm atom} (w_n) \leq C h_n^4$. Hence,

which tends to $0$ by assumption. This implies that $w_n \in \mathcal{S}_{\delta/2}$ for $n$ large enough.

Acknowledgments

This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under project number 285722765, as well as the Engineering and Physical Sciences Research Council (EPSRC) under the grant EP/R043612/1.