\section{Introduction}

The formation of wrinkles and their localization into folds is observed when a layered material is progressively compressed. As such, they occur in a variety of systems at different scales in Nature from the folding of geological strata \cite{Pollard2005} to  the morphogenesis of biological tissues such as the cerebral cortex\cite{VanEssen1997, Toro2005, Tallinen2016} or fingerprints \cite{Kucken2004} as well as in the growth of biofilms \cite{Trejo2013}. Just as nature uses these mechanical instabilities to generate intricate patterns, wrinkling instabilities are also finding a number of technological applications  from flexible electronic devices \cite{Rogers2010, Kim2010a} to controlled patterned surfaces \cite{Bowden1998, Schweikart2009} that have improved light harvesting efficiency \cite{Kim2012b} or surface hydrophobicity \cite{Lin2015}. Because of this ubiquity and the renewed interest in taking control of elastic instability, the deformation of layered elastic materials has been the object of renewed interest recently, with a number of fundamental studies in the past decade \cite{Cerda2003, Bo2012, Reis2009, Brau2011, Kim2011}.

The simplest example of elastic instability is the classic Euler buckling \cite{Levien2008,Landau} in which a compressed beam buckles over its entire length. To use elastic instability for pattern formation usually requires the selection of a length scale other than the system size, however. Perhaps the simplest example of such a system is provided by an elastic film (of bending stiffness per unit width $B$) floating on a liquid of density $\rho$ and subject to an axial compression. This system quickly forms wrinkles with a wavelength $\lambda \sim (B/\rho g)^{1/4}$ --- a result that expresses the compromise between the bending rigidity of the film (which prefers large-amplitude wrinkles) and the weight of the liquid ($g$ being the gravitational acceleration), which prefers many small-amplitude wrinkles\cite{Pocivavsek2008, Holmes2010, King2012, Pineirua2013}. This wrinkling instability can also be observed in other types of floating materials such as thin layers of nanoparticles\cite{Leahy2010}, surfactants \cite{Lu2002, Boatwright2010} or particle rafts \cite{Vella2004}. However, at some point most of these wrinkles begin to shrink with a very small number growing to form much larger `folds'; although other options (such as delamination\cite{Wagner2011}) exist, folds are more generically observed and so we focus on them in this paper.

The nature of the wrinkle-to-fold transition in this system has been somewhat controversial: it was initially suggested (on the basis of experiments) that folds form at a critical compression. Several numerical and analytical works have studied the wrinkle-to-fold transition in idealized situations (in which  the elastic film is infinitely long and weightless \cite{Diamant2011, Diamant2013, Brau2013, Rivetti2013}) and have shown that in these scenarios the transition is not sharp: i.e.~folds emerge continuously from the wrinkled state. Recently, a more realistic model (incorporating the finite length of the film) demonstrated that this transition emerges at a finite compression making a second-order transition \cite{Rivetti2014, Oshri2015}. Despite the amount of work on the wrinkle-to-fold transition, several questions remain unanswered. For example, existing models suggest that downwards and upwards folds should be equally favourable, while experiments report exclusively downward folds. Moreover, what happens beyond the transition between these two states has not been described experimentally but only via numerical simulations\cite{Demery2014} with the weight of the film discussed in scaling terms. In the first part of this paper we investigate the effect of the film weight on the wrinkle-to-fold transition while in the second part we focus on the evolution of the fold upon further compression. We observe that the fold becomes a closed loop after self-contact, encapsulating a volume of the upper fluid. As the applied compression is increased, we observe that the fold either sinks deeper without limit or eventually bends upward, depending on the relative buoyancy of the drop-film combination. We predict the regime diagram for this behaviour theoretically.

\section{Experimental section}
\textbf{Elastic film production:} Thin elastic films are produced by spincoating vinylpolysiloxane VPS (elite double 32, \texttt{Zhermack}) of base density $\rho_{VPS}=1.20\mathrm{~g\,cm^{-3}}$. To produce films of various densities, we add iron powder (97\%, -325 mesh, \texttt{Sigma Aldricht}) $\rho_{Fe}=7.87 \mathrm{~g\,cm^{-3}}$ to the non cross-linked polymer. This allows us to work with films of densities $\rho_s=1.2-2.6\mathrm{~g\,cm^{-3}}$ (see table, SI). The spin-coated polymer/iron mixture viscosity changes with the iron powder concentration so that at a fixed rotation speed the film thickness, $t$, obtained varies. We adjust the rotation speed to narrow the thickness range ($50\mathrm{~\mu m}<t<120 \mathrm{~\mu m}$). Two sets of film lengths $L_0$ and width $W$ are produced: $L_0=90\mathrm{~mm}$, $W=60\mathrm{~mm}$ and $L_0=75\mathrm{~mm}$, $W=50\mathrm{~mm}$.

\noindent\textbf{Compression experiment:}

\begin{figure}
\centering
  \includegraphics[width=8.3cm]{Fig1.pdf}
  \caption{Schematic of the experimental setup defining the length $L$ between the boundaries, the wrinkle wavelength $\lambda$ and the wrinkle amplitude $A$ as well as the Cartesian and intrinsic coordinates used in the model: ($x$,$y$) and ($s$,$\theta$).}
  \label{fgr:manip}
\end{figure}

The experiments are conducted in a glass tank ($12 \times 11 \times 6.5\mathrm{~cm^3}$) with two parallel polymethylmethacrylate (PMMA) plates with horizontal protruding ends (Fig.~\ref{fgr:manip}). The tank is filled with tap water to a level higher than the protruding ends. The elastic film is then carefully placed at the air/water interface between the plates. The water level is then lowered (using a syringe) until the edges of the film comes into contact with the protruding ends of the PMMA plates. The film naturally adheres strongly to the PMMA achieving a clamped boundary condition at the protruding ends. The water level is adjusted with the syringe so that the film is completely flat. A light mineral oil ($\rho_o=0.836\mathrm{~g\,cm^{-3}}$, \texttt{Sigma Aldrich}) is then poured slowly on top of the film (so that no oil invades the lower surface of the film, in contact with water). One of the PMMA plates is mounted on two perpendicular manual translation stages for alignment in the ($y$,$z$) direction and connected to a stepper motor (\texttt{Thorlabs}) of micrometer precision. The compression is quasi-static: the stepper motor displaces the plate in small increments at a constant speed. The motor stops for $5\mathrm{s}$ between each step allowing the system to relax to its equilibrium shape. The elastic film is imaged from the side and/or the top with two \texttt{Nikon} D800-E cameras mounted with macro objectives ($105$ mm). The images are then analysed using either \texttt{ImageJ} or \texttt{MATLAB}.

\noindent\textbf{Film analysis:}

We use two methods to determine the elastic films thicknesses. We weigh the films on a milligram scale: knowing the density, length and width we deduce an average thickness. We also illuminate the films with a tilted laser line. The laser deflection being a function of the height and tilting angle, height profiles are then extracted using the appropriate calibration. Both measurement techniques give similar results. In addition, the laser line allows us to check that our films are uniform: thickness variations are below the measurement standard deviation ($\pm\:10\mathrm{~\mu m}$) except at the edges of the film where a small ridge ($20\mathrm{~\mu m}$ thicker) develops during the spin coating/curing. The Young's modulus $E$ and Poisson ratio $\nu$ are determined by a tensile test on a dogbone shape using a \texttt{Shimadzu} tensile machine. For the pure VPS we find $E=1.0\mathrm{~MPa}$, $\nu=0.5$ while for the VPS/iron particle mixture the Young's modulus increases up  to $E=2.9\mathrm{~MPa}$ for a dogbone of density $\rho_s=2.5\mathrm{~g\,cm^{-3}}$. We checked that the film Young's modulus is the same as that obtained from the test samples by measuring the deflection of the film under its own weight. By letting a small portion of the film protrude out of a clamp we were able to vary the length easily and fit the bending stiffness $B=Et^3/\big[12(1-\nu^2)\bigr]$. The values found through this procedure are in good agreement with the previously measured $E$ and $t$. Finally the measured $E$ and $t$ are compatible with the direct measurement of the wrinkle wavelength $\lambda$ (see table, SI).

\section{Results}

We compress an elastic film floating at an interface between two fluids and study its behaviour as we increase the imposed displacement. At zero compression, $L=L_0$, the film lies flat at the interface. As soon as we start compressing the film ($\Delta=L_0-L>0$), it buckles out of plane with a characteristic wavelength $\lambda$ that develops along the film length (Fig.~\ref{FigEvolution} (a)). This is the wrinkled state. As $\Delta$ increases, the wrinkle amplitude  $A$ initially grows uniformly (Fig.~\ref{FigEvolution} (b)) until, suddenly, only one of the wrinkles continues to grow while the other wrinkles   progressively vanish. The excess  length of the vanished wrinkles is absorbed by the remaining  structure, the fold (Fig.~\ref{FigEvolution}(c)). We therefore observe a clear transition between the ``wrinkled state", in which the deformation is distributed throughout the film, and a ``fold state", in which all the deformation is localized in a narrow region of high curvature, \textit{i.e.}: the fold. As the compression $\Delta$ is increased beyond the wrinkle-to-fold transition, the fold continues to grow in amplitude and its curvature increases. Finally, the two opposing edges of the fold come  into contact  (self-contact). At this point a horizontal column of the upper fluid is encapsulated in this ``teardrop" shape (Fig.~\ref{FigEvolution} (d)). With further increases in compression this teardrop is forced down into the subphase and does not noticeably change shape; instead the length over which the film is in self-contact increases (Fig.~\ref{FigEvolution} (e)). In this paper, we  study the role of the density of the elastic film in the wrinkle-to-fold transition and its impact on the evolution of the fold as the uniaxial compression is increased to the point of self-contact and beyond.

\subsection{Wrinkle-to-fold transition and fold before self-contact}

 \begin{figure}
  \includegraphics[width=8.3cm]{Fig2.pdf}
  \caption{Side images of an elastic film of density $\rho_s=1.8\:g.cm^{-3}$ at an oil/water interface. Compression $\Delta$ increases from \textbf{(a)} to \textbf{(e)}, scale bar 5mm. Colours have been inverted to enhance contrast. \textbf{(a)-(b)} The film displays a quasi-periodic out of plane displacement of amplitude $A$ and wavelength $\lambda$, referred to as the wrinkled state. \textbf{(c)} The deformation localizes in a single fold. \textbf{(d)} The fold reaches self-contact. \textbf{(e)} A column of oil is encapsulated in the fold which grows deeper towards the bottom of the tank as compression is increased.
  }
  \label{FigEvolution}
\end{figure}

 In order for the fold to nucleate near the centre of the elastic film, at $L_0/2$, we ensure that the clamped edges are aligned very carefully. (The fold has to appear at least one wavelength $\lambda$ away from the clamped edges to avoid any boundary effects.) The results were compared to a fold generated by applying a small pressure at the film centre before starting to compress it. For three different densities $\rho_s=1.2,\:1.4,\:1.8\mathrm{~g\,cm^{-3}}$ and two different film lengths $L_0$, we used  images taken from the side to measure the fold amplitude $A$ during compression up to self-contact. The data are then normalized by the wavelength $\lambda$ measured experimentally (Fig.~\ref{fgr:Amppresc}). For small compressions, $\Delta / \lambda \lesssim 0.3$ we observe the wrinkled state with no influence of the film weight. In this wrinkled regime, Pocivavsek \textit{et al.}~\cite{Pocivavsek2008} found that the evolution of the wrinkle amplitude exhibits a weak dependence  on the film length, which we are not able to observe in our experiments. Moreover, Rivetti and Neukirch \cite{Rivetti2014} also predict that the buckling modes (symmetric/antisymmetric/non symmetric) of the film can evolve during the compression for different film lengths. We do observe two different buckling routes in our experiments and they are responsible for the two different trends in the amplitude that we observe in Fig.~\ref{fgr:Amppresc}: the experiments at an air--water interface remain symmetric throughout the compression while the experiments at an oil--water interface start antisymmetric but switch continuously to a symmetric mode. Since the antisymmetric mode has a smaller amplitude (for a given compression) than the symmetric one, when the fold switches to the symmetric mode its amplitude increases rapidly. This results in an inflection point at $\Delta/\lambda\sim 0.3$ for experiments at an oil--water interface. However, these details  are very sensitive to the experimental conditions (particularly to the clamped boundary conditions) and no quantitative comparison could be made with the predictions of Rivetti and Neukirch \cite{Rivetti2014} (see SI). For larger compressions,  $\Delta / \lambda > 0.4$, the behaviour does not depend on film length: the deformation always localizes in a downward symmetric fold with an amplitude that grows linearly with compression, as described previously \cite{Pocivavsek2008, Diamant2011}.

 \begin{figure}
  \centering
  \includegraphics[width=8.3cm]{Fig3.pdf}
  \caption{Dimensionless wrinkle/fold amplitude as a function of the dimensionless compression up to self-contact. (\protect\includegraphics[height=2mm]{op_purp_circ.pdf}, \protect\includegraphics[height=2mm]{op_orange_circ.pdf}) represent oil/water experiments, (\protect\includegraphics[height=2mm]{filled_blue_circ.pdf}, \protect\includegraphics[height=2mm]{filled_green_circ.pdf}) air/water experiments. Three film densities are presented (with different lengths/widths/thicknesses), giving rise to the four dimensionless weights presented in the legend. The black curve corresponds to the symmetric solution \cite{Diamant2011} ($M=0$), the orange and purple curves correspond to the symmetric numerical solution of equation \eqref{nlsystM} for $M=0.10$ and $M=0.18$. \textit{Inset}: Top view of the compressed film. The fold reaches self-contact at the edges of the film but is still open in the centre (see SI). Scale bar: 5mm.}
  \label{fgr:Amppresc}
\end{figure}

However, as we vary the film densities we find that the amplitude of the fold increases slightly with the film mass (Fig.~\ref{fgr:Amppresc}). To explain this dependence we first turn to the model developed by Diamant and Witten \cite{Diamant2011} for the nonlinear deflections of a floating elastic film. In the limit of an incompressible, weightless film of infinite length, this model has  an analytical solution\cite{Diamant2011}; here we shall have to develop a numerical solution but we first outline the derivation of the governing equation for a light film, so that the appropriate modifications can be made for the film weight. We introduce the intrinsic coordinates ($s$,$\theta$), where $\theta$ is the angle between the film and the horizontal axis and $s$ is the arc-length;  the film centreline is then parametrized in terms of arc-length, $[x(s),y(s)]$. The energy $U$ of a film of bending stiffness $B$ and width $W$ contains contributions from bending $U_b=\tfrac{1}{2}BW\int_{-\infty}^{\infty} \left(\partial_s \theta\right)^2~\mathrm{d} s$ and from the gravitational potential energy of the underlying fluid substrate $U_s$. In our case the gravitational energy of the upper and lower fluids  $U_{l_{up}}$ and $U_{l_{low}}$, respectively are given by: $U_s=U_{l_{up}}+U_{l_{low}}=\tfrac{1}{2}\Delta\rho gW\int_{-\infty}^{\infty} y(s)^2\cos{\theta}~\mathrm{d} s$ where the relevant density is the density difference between the two fluids $\Delta\rho=\rho_{low}-\rho_{up}$. The energy is to be minimized subject to the constraint of an imposed compression $\Delta=\int_{-\infty}^{\infty}(1-\cos\theta)~\mathrm{d} s$. To facilitate the calculation, lengths are  non-dimensionalized by $\ell_{eh}=\bigl[B/(\Delta\rho g)\bigr]^{1/4}=\lambda/(2\pi)$ and energies by $WB/\ell_{eh}$. The minimization itself is discussed in detail by Diamant and Witten\cite{Diamant2011} and yields a single equation for the intrinsic angle $\theta(s)$ with the boundary conditions $\theta (\pm \infty) = \partial_s \theta (\pm \infty) = y (\pm \infty) =0 $:
\begin{equation}
\partial_s^4 \theta + \bigl[\tfrac{3}{2}\left(\partial_s \theta\right)^2+P\bigr] \partial_s^2 \theta + \sin \theta=0
\label{Witteneq}
\end{equation}
Here $P$ is the Lagrange multiplier associated with the compression constraint, and  corresponds physically to the dimensionless horizontal force applied on the elastic film by the compression. Equation \eqref{Witteneq} has a family of analytical solutions \cite{Diamant2011, Diamant2013, Rivetti2013}:
\begin{equation}
\begin{split}
\theta & (s)=4\arctan \bigg( \frac{\kappa \sin\bigl[k(s+\phi)\bigr]}{k \cosh(\kappa s)} \bigg)\\
k=\tfrac{1}{2} & \sqrt{2+P} , \quad \kappa=\tfrac{1}{2}\sqrt{2-P}, \quad P=2-\frac{\Delta^2}{16}.
\end{split}
\label{Wittensol}
\end{equation} This family of solutions is parametrized by $\phi$, with $0 \le k\phi \le \pi$; the value of $k\phi$ selects the symmetry of the profile\cite{Rivetti2013}, with $k\phi=0$ corresponding to an even fold and $k\phi=\pi/2$ corresponding to an odd fold\cite{Diamant2011}.

For the experiments presented in this paper, the weight of the film is not negligible; we must therefore supplement  the bending and substrate energy with the  gravitational energy of the film itself, $U_g$. If $\rho_s$ is the density of the film and $\rho_{w}$ the density of the lower liquid (water) then the dimensionless gravitational energy of the film is
\[U_g = M \int_{-\infty}^{\infty} y(s)~\mathrm{d} s,\]
where
\begin{equation}
M=\frac{(\rho_s-\rho_w)t}{\Delta\rho \ell_{eh}}
\label{eqn:Mdefn}
\end{equation} is a dimensionless number that compares the weight of the film to the restoring force provided by buoyancy over the horizontal length $\ell_{eh}$. Applying the energy minimization procedure described by \cite{Diamant2011} to the total energy $U=U_s+U_b+U_g$  we obtain the system of equations (see SI~for details):
\begin{eqnarray}
&\partial_s^4 \theta + \left[\tfrac{3}{2}\left(\partial_s \theta\right)^2+P-My\right] \partial_s^2 \theta +(1-2M\partial_s \theta)\sin \theta=0 \nonumber \\
&\partial_s x=\cos \theta, \quad \partial_s y=\sin \theta \label{nlsystM}
\end{eqnarray}

using the same clamped boundary conditions as for the light film. We note that, as expected, the weightless equation  \eqref{Witteneq} is recovered by setting $M=0$.

For $M\neq0$, we are not able to find an analytical solution of the system of equations \eqref{nlsystM}; instead we obtain numerical solutions by using the \texttt{MATLAB} routine \texttt{bvp4c}, which uses a relaxation technique with the analytical solution \eqref{Wittensol} as an initial guess. We also use a simple  continuation algorithm to find the solutions as parameters are varied.

Our numerical solutions suggest that the antisymmetric solution no longer exists once $M\neq0$ (we are unable to find numerical solutions with this symmetry).  Even if an antisymmetric solution of equation \eqref{nlsystM} does exist, with $M>0$ (respectively $M<0$), the downward (respectively upward) symmetric solution has a lower energy due to the contribution of the gravitational term $U_g$. (In the limit $M=0$,  all of the solutions \eqref{Wittensol} are known to have the same energy \cite{Diamant2011, Diamant2013, Brau2013, Rivetti2013}.) The weight of the film therefore lifts the degeneracy that was inferred from the analytical solution \eqref{Wittensol} but this has not been observed experimentally\cite{Pocivavsek2008}; previously this breaking of symmetry was attributed to surface tension, but not quantified \cite{Pocivavsek2008}.

Experimentally, we do, in fact, observe antisymmetric configurations for sufficiently low compressions ($\Delta / \lambda \lesssim 0.25$); however, as the compression increases, the film always evolves to a downward symmetric configuration (see SI). This is in accordance with the theoretical prediction that  the energy difference between the two states increases  with $\Delta / \lambda$.

The evolution of the amplitude with increasing compression within the symmetric solution, as determined from the numerical solution of \eqref{nlsystM} for $M=0,\:0.10,\:0.18$, is shown in Fig.~\ref{fgr:Amppresc}. For $\Delta / \lambda \lesssim 0.2$, the system remains in the wrinkled state and there is very little influence of the film weight on the amplitude--compression curve. However, with $\Delta / \lambda > 0.25$, the fold amplitude increases with the film weight $M$, holding  compression, $\Delta$, fixed; furthermore, we recover the linear evolution of the fold amplitude with $\Delta$. As $M$ varies, the experimentally determined $(\Delta,A)$ curves are shifted as predicted by our theory (Fig.~\ref{fgr:Amppresc}). The  observed experimental fold is fully localized and downward symmetric for $\Delta / \lambda > 0.4$. Its amplitude can thus be quantitatively predicted by the numerical solution of \eqref{nlsystM} in this regime.

In fig.~\ref{fgr:Amppresc}, the last point of each curve corresponds to the point at which the downward symmetric fold first self-contacts, forming  a loop, \textit{i.e.}~ the teardrop shape. Experimentally, we find that this point is reached for $\Delta_{sc}/\lambda \approx 0.8$ whereas using equation \eqref{nlsystM}, this transition is found at a slightly higher value ($\Delta_{sc}/\lambda \approx 0.89$ for symmetric folds). The reason for this discrepancy is not immediately clear.
Some insight into this shift is obtained by examining a top view of the fold (inset of Fig.~\ref{fgr:Amppresc}): its shape is not the same at the edge of the film and at its centre . Capillary forces pull the elastic film on its sides closing the fold tighter at its edges than at its centre. Surface tension slightly distorts the fold's shape at the sheet's edges, which explains the small discrepancy between the theoretical and experimental $\Delta_{sc}$. However the fold's amplitude remains constant even when subjected to surface tension effects (see SI~for details).

\subsection{Evolution of the fold}

The analytical solutions \eqref{Wittensol} and the numerical solutions of \eqref{nlsystM} cease to be relevant beyond self-contact ($\Delta>\Delta_{sc}$) since, within the limitations of our model, the numerically generated profiles interpenetrate. Experimentally, we observe two different behaviours depending on the dimensionless mass $M$, defined in \eqref{eqn:Mdefn}, as shown in Fig.~\ref{fgr:amppostsc} (a)-(f): heavy films ($M\gtrsim0.14$) retain a symmetric configuration and encapsulate a column of the upper fluid in a teardrop shape fold (Fig.~\ref{fgr:amppostsc} (d)-(e)) that grows deeper as the compression increases (Fig.~\ref{fgr:amppostsc} (f)). Lighter films ($M\lesssim0.14$) start symmetric with the same teardrop shape but, as the compression increases, the loop starts to tilt and grows back up towards the interface (Fig.~\ref{fgr:amppostsc} (a)-(b)) where it finally reaches an obstacle such as the clamped edges or itself (Fig.~\ref{fgr:amppostsc} (c)). To quantify these observations, we first  define the amplitude after self-contact, $A^*$, as the depth of the centre of the loop, see the inset of Fig.~\ref{fgr:amppostsc} (g). This quantity is clearly defined beyond self-contact, while before self-contact,  $A^*\equiv A$. Fig.~\ref{fgr:amppostsc} (g) shows the experimentally measured amplitude as a function of the compression before and after self-contact. We vary the film density and length and present our results using the dimensionless mass $M$ described previously. At the scale of Fig.~\ref{fgr:amppostsc} the variation of amplitude due to the mass $M$ is not visible before self-contact; immediately after self-contact the amplitude initially keeps growing linearly ($A^*/\lambda \approx 0.5\Delta/\lambda$). However, as the compression is increased beyond self-contact, the influence of $M$ becomes more apparent, at a critical compression $\Delta_b$  the fold begins to tilt back up towards the interface with increasing compression \textit{i.e.}: $A^*$ decreases. (We emphasize that this tilting occurs quasi-statically.) Fig.~\ref{fgr:amppostsc} (g) shows that as $M$ increases, the transition from a straight to a tilted fold occurs at larger $\Delta_b$. Finally, when the elastic film is sufficiently dense ($M=0.18$ in Fig.~\ref{fgr:amppostsc} (g)), the fold  never bends back up and sinks indefinitely into the liquid subphase.

\begin{figure}
\centering
  \includegraphics[width=8.3cm]{Fig4.pdf}
  \caption{Images from the side of two folds as compression is increased. The two sets of pictures have the same dimensionless compressions, from left to right $\Delta/\lambda \approx \: 0.93,\: 1.19,\: 1.44$. Scalebars: $5\mathrm{~mm}$. \textbf{(a)-(c)} film with a low mass $M=0.03$, $\rho_s=1.20\mathrm{~g\,cm^{-3}}$. \textbf{(d)-(e)} film with a high mass  $M=0.18$, $\rho_s=1.8\mathrm{~g\,cm^{-3}}$. Colours have been inverted to enhance contrast. \textbf{(g)} \textit{Inset}: Schematic presenting the post buckling fold amplitude $A^*$ (fold centre in red). Dimensionless wrinkle/fold amplitude $A^*$ as a function of the dimensionless compression after self-contact. (\protect\includegraphics[height=2mm]{blue_open_circ.pdf}, \protect\includegraphics[height=2mm]{red_open_circ.pdf}, \protect\includegraphics[height=2mm]{op_purp_circ.pdf}) represents oil/water experiments, (\protect\includegraphics[height=2mm]{green_close_circ.pdf}, \protect\includegraphics[height=2mm]{brown_close_circ.pdf}) air/water experiments. \protect\includegraphics[height=2mm]{filled_black_circ.pdf} show the data corresponding to the pictures \textbf{(a)-(f)}. }
  \label{fgr:amppostsc}
\end{figure}

To understand the tilting behaviour of the loop, we first consider its properties at the point of self-contact. The key quantities of interest are the loop's  width $w$ and height $h$. Experiments show that these grow linearly with $\lambda$ (see Fig.~\ref{fgr:formebouclelarge} (b)) but do not vary with compression beyond self-contact  (see inset of Fig.~\ref{fgr:formebouclelarge} (b)). We believe that this lack of evolution of the teardrop shape is due mainly to the constant volume of the upper liquid that is trapped in the loop after the point of self-contact. The teardrop shape observed corresponds to that predicted by the solution of equations \eqref{Wittensol} and \eqref{nlsystM} \emph{at} self-contact (Fig.~\ref{fgr:formebouclelarge} (a)). However, the values of $w$ and $h$ are slightly overestimated by the model since we can only measure the shape at the film edges, where surface tension effects come into play;  at the centre of the film the loop is wider (see inset Fig.~\ref{fgr:Amppresc}). The numerical solution of equation \eqref{nlsystM} predicts that the loop formed at self-contact shrinks as $M$ increases. We were not able to observe this experimentally, and conclude that  this effect must be smaller than the role of surface tension at the film edges.

We compare our results to the numerical work by Demery \textit{et al.} \cite{Demery2014} where they study a compressed film  after self-contact; however, they do not incorporate the weight of the film itself. They find two possible configurations for the shape of the film in this case: a symmetric one (which we also observe experimentally) and an antisymmetric one (which we have never observed). In the weightless case the antisymmetric configuration has a lower energy --- its energy saturates with increasing $\Delta$, while the energy of the symmetric downward fold \emph{increases} linearly with $\Delta$. In our case the added weight of the film reduces the energy of the symmetric configuration by $U_g\sim-M \Delta^2$, explaining why this configuration is ultimately energetically favourable.

\begin{figure}
\centering
  \includegraphics[width=0.45\textwidth]{Fig5.pdf}
  \caption{\textbf{(a)} Comparison of symmetric profiles at self-contact given by the analytical solution \eqref{Wittensol} for $M=0$ (black) and the shape for $M=0.18$, determined from the numerical solution of equation \eqref{nlsystM} (red). The arrows define the height $h$ and the width $w$ of the teardrop. \textbf{(b)} The teardrop height $h$ (red circles) and width $w$ (blue squares) as a function of the wavelength. Closed symbols represent air/water experiments, open symbols oil/water experiments; here, solid lines correspond to the analytical prediction, based on equation \eqref{Wittensol}, while dashed lines correspond to the numerical solution of equation \eqref{nlsystM} for $M=0.18$. All other experimental parameters are found to be not relevant and vary across the data. \textit{Inset}: Experimental teardrop height and width in one experiment ($M=0.14$, $\rho_s=1.8\mathrm{~g\,cm^{-3}}$) as a function of the compression.}
  \label{fgr:formebouclelarge}
\end{figure}

To describe the evolution of the fold when $\Delta>\Delta_{sc}$ in more detail,  we need a new mathematical model. We continue to neglect surface tension  so that the film behaviour is invariant along its width and the problem can be treated as two-dimensional. Our model of the elastic film after self-contact is illustrated schematically in Fig.~\ref{fgr:heavyelexpdiagphase}. We break the system in two parts: a heavy beam corresponding to the part of the film in self-contact and a force $F$ acting at its tip (which is the teardrop shape encapsulating buoyant fluid). Here
\begin{equation}
F=\Delta\rho g \mathscr{A}-(\rho_s-\rho_{w}) g t \mathscr{L} ,
\label{fgr:force}
\end{equation}
$\mathscr{A}$ and  $\mathscr{L}$ are respectively the half-area and half the perimeter of  the upper fluid encapsulated in the teardrop (see schematic in the inset of Fig.~\ref{fgr:heavyelexpdiagphase}). Considering these half quantities means that we are considering each half of the two that are in contact separately, and assuming they do not exert force on one another. While this reduction suggests that the problem is equivalent to that of a beam subject to a constant load at one end, we emphasize that the self-weight of the beam is important and so instead we must consider a `heavy hanging column'\cite{Wang1981}, subject to a constant force pushing at the tip (since we have already found experimentally that the teardrop size does not evolve with further compression). We assume a clamped boundary condition at the top for simplicity (schematic Fig.~\ref{fgr:heavyelexpdiagphase}). At each compression step the portion in self-contact $L^*$ grows, increasing the length of the effective beam.

We introduce again  intrinsic coordinates ($s$,$\theta$), with $s$ the arc-length but now $\theta$ denotes the angle between the heavy beam and the vertical axis. The $(x, y)$ coordinates are thus rotated clockwise by $90$ degrees compared to the model presented in the previous section. In this system, the equation for the heavy hanging column is given by \cite{Wang1981, Wang1986}(see SI~for details):
 \[ B \partial_s^2 \theta = -[F-(\rho_s-\rho_{w})gts] \sin \theta\]
 $B$ is the beam bending modulus and $t$ its thickness. $\rho_s$ is the beam density and $\rho_{w}$ the lower liquid (water) density. The boundary conditions are that the end at  $s=0$ is free while that at  $s=L^*$ is clamped:
 \[ \partial_s \theta(s=0)=0, \quad \theta(s=L^*)=0\]
 The system is made dimensionless by dividing $s$ by $L^*$:
 \begin{equation}
\begin{split}
\partial_s^2 \theta + \left[\widetilde{F}-\bigg(\frac{L^*}{\ell_g}\bigg)^3 s\right]\sin \theta =0 \\
\end{split}
\label{heavyeleq}
\end{equation}
\[\partial_s \theta(s=0)=0, \quad \theta(s=1)=0.\]
Here $\widetilde{F}=F{L^{*}}^2/B$ is the dimensionless force due to the buoyancy of the teardrop and $\ell_g=(B/[(\rho_s-\rho_{w})gt])^{1/3}$ is the elasto-gravitational length, which compares elastic and gravitational forces; loosely speaking the elasto-gravitational length is the length above which the film will buckle under its own weight. The elasto-gravitational length is also closely related to the parameter $M$, which is defined in \eqref{eqn:Mdefn} and may be written $M=(\ell_{eh}/\ell_g)^3$.

For a fixed applied force $F$, the beam will buckle once the length reaches a threshold value. This threshold can be determined by linearizing \eqref{heavyeleq} and solving the resulting problem analytically  \cite{Wang1981} (see SI). This analysis yields the critical force $\widetilde{F}_c$ needed for the beam to buckle as a function of $(L^*/\ell_g)^3$, \textit{i.e.}~$\widetilde{F}_c=\frac{F_c{L^*}^2}{B}=\operatorname{f}\big((L^*/\ell_g)^3\big)$. In our experiments, the control parameter is the beam length $L^*$, which increases with increasing compression; the force is constant, since it originates from the buoyancy of the teardrop, which is fixed during the compression. We therefore write:
\begin{equation}
\frac{F\ell_g^2}{B}=\Big(\frac{\ell_g}{L^*_c}\Big)^2\operatorname{f}\bigg(\frac{L^*_c}{\ell_g}\bigg)
\label{heavyelthreshold2}
\end{equation}
where the function $\operatorname{f}(x)$ emerges from a solvability condition, see SI. For a given teardrop size, and hence buoyancy force, equation \eqref{heavyelthreshold2} may be solved to give the critical length at which the beam begins to bend.

\begin{figure}
\centering
  \includegraphics[width=8.3cm]{Fig6.pdf}
  \caption{Phase diagram of the post self-contact fold buckling. Green filled circles give the critical length $L^*$ before buckling in each experiment. Black open circles represent experiments in which the fold never buckles. The black dashed curve and background colours are guides for the eye to distinguish experimental phases. The red solid curve is the analytical result of the heavy hanging column (equation \eqref{heavyelthreshold2}), the blue (respectively orange) solid curve is the numerical result of the heavy hanging column theory with an initial angle of 10$^{\circ}$ (respectively 20$^{\circ}$). \textit{Inset}: Schematic of a perfectly symmetric fold post self-contact. The fold is split in two parts: the portion in self-contact is treated as a heavy hanging column of length $L^*$ and thickness $t$ and the teardrop is modeled as a constant point force (resulting from the buoyancy of the teardrop). The upper part is not taken into account here and a clamped boundary condition is assumed at the top of the heavy hanging column. The intrinsic coordinates ($s$,$\theta$) are drawn.}
  \label{fgr:heavyelexpdiagphase}
\end{figure}

The computed prediction for the length at which the fold should start to bend upwards is presented in Fig.~\ref{fgr:heavyelexpdiagphase}. Here, the buoyancy force of the teardrop is calculated using equation \eqref{fgr:force} with $\mathscr{A}$, $\mathscr{L}$ and $L^*$ measured from images taken from the side (and are hence subject to the inaccuracies discussed earlier) and $B$ is calculated based on the measurement of the film wrinkling wavelength ($B=\Delta\rho g \bigl[\lambda/(2 \pi)\bigr]^4$). This calculation, via equation \eqref{heavyelthreshold2}, captures qualitatively the experimentally observed fold behaviour. In particular, there is a critical buoyancy force below which the fold never buckles, regardless of the length of the film.

The experimental data points shown in Fig.~\ref{fgr:heavyelexpdiagphase} lie well below the prediction of equation \eqref{heavyelthreshold2}, although the trend is similar. We believe that in most experiments this discrepancy occurs because the fold is not perfectly aligned with the vertical axis when it reaches self-contact (Fig.~\ref{fgr:amppostsc}(a)). The average angle $\alpha$ between the fold and the vertical axis is $\left|\alpha\right|=10^{\circ}$. To account for this angle, we may include this effect by changing the boundary condition at the top of the beam: $\theta(s=1)=\alpha$. We now solve the problem numerically (using equation \eqref{heavyeleq}) to determine the critical length  $L^*_c/\ell_g$ for buckling (see SI~for details). The results are plotted in Fig.~\ref{fgr:heavyelexpdiagphase} (blue curve: $\alpha=10^{\circ}$, orange curve: $\alpha=20^{\circ}$) and are closer to the experimental values. An additional source of discrepancy between the experimental results and the numerical model may be, once again,  that we underestimate the size of the teardrop, since we measure its shape at the edge of the film (where it is affected by surface tension) and not at its centre. It may also be due to how we oversimplify the boundary condition at the top of the beam. Here, we have assumed that  a clamped boundary condition is appropriate but we have observed that the top of the fold can slide and rotate slightly as the compression is varied.

\section{Conclusion}
We have studied experimentally the behaviour of thin elastic films under uniaxial compression at a liquid-fluid interface. As previously observed, when confinement increases the film undergoes a transition from a uniform wrinkled state to a localized configuration in which a single fold is observed. We showed that the film weight affects the wrinkle-to-fold transition, breaking the up-down symmetry of the light film problem considered previously. More importantly, however, the film weight plays a large role in the evolution of the fold beyond the transition: it may bend back towards the fluid interface or grow deeper as compression is increased. A simple model predicts the two possible behaviours in a phase diagram.

The experimental results presented here allow us to confirm some  previous theoretical assumptions but also highlight the role of the film weight and surface tension in the formation and evolution of a fold. The teardrop shape observed after self-contact seems similar to the so-called ``self-encapsulation" of an elastic rod which is observed in a solely elastic situation when a rod covering a fixed span is loaded at the middle with a transverse force such that two points of the rod come into contact with each other\cite{Bosi2015}. Finally, this study could give insight into the elastic properties of compressed particle-laden interfaces, which exhibit a wrinkling instability similar to that described for elastic films in this Letter. For example, our results show that such films may encapsulate a teardrop of the upper (lower) fluid if the dimensionless parameter $M>0$ ($M<0$); this may provide a novel route for the creation of cylindrical, particle-coated droplets \cite{Subramaniam2005} by compressing a particle-laden interface. In such settings, the granular aspect of the interface allows the interface's mechanical properties to be tuned by probing new parameters such as the particle size, density or packing fraction \cite{Petit2016} as well as the solid-liquid contact angles. However, the granular aspect also affects the stress state within the layer at the onset of wrinkling \cite{Cicuta2009} meaning that the transition from wrinkles to localized folds upon further compression remains to be carefully characterized in such systems.

\section{Acknowledgement}
We are grateful to S\'ebastien Neukirch, Vincent D\'emery and Claude Perdigou for useful discussions as well as Aur\'elie Fargette for her help with the preparation of thin films and Julien Chopin for his expertise with the tensile testing.