Open Access Article
Kyoko Shibata†
a,
Hana Kanda†b,
Yoshimi Tanakac and
Yutaka Sumino
*ade
aDepartment of Applied Physics, Tokyo University of Science, 6-3-1 Nijuku, Katsushika-ku, Tokyo 125-8585, Japan. E-mail: ysumino@rs.tus.ac.jp
bGraduate School of Innovative and Practical Studies, Yokohama National University, Yokohama 240-8501, Japan
cResearch Collaborator, Institute of Science and Engineering, Kanazawa University, Kanazawa 920-1192, Japan
dWaTUS and DCIS, Research Institute for Science & Technology, Tokyo University of Science, 6-3-1 Nijuku, Katsushika-ku, Tokyo 125-8585, Japan
eFaculty of Engineering and Physical Sciences, University of Surrey, Guildford, Surrey GU2 7XH, UK
First published on 5th June 2026
We investigate the encapsulation of water by a thin elastic film as a minimal model of elastocapillary self-folding with fluid transport. An equilateral triangular polydimethylsiloxane film is lifted quasi-statically from a water surface, while its side length and thickness are systematically varied. Depending on these parameters, the film exhibits three distinct morphologies: folding, recoiling, and liquid encapsulation. We show that the morphology is governed by the interplay of surface, gravitational, and bending energies, and that encapsulation occurs only within a narrow parameter region where the elastocapillary, elastogravity, and capillary length scales become comparable. This provides a simple physical criterion for liquid encapsulation by elastic films.
. Under such conditions, capillary forces can drive counterintuitive fluid motions, such as capillary rise and the spontaneous climbing of droplets against gravity.1–4 These phenomena originate from the tendency of interfacial tension to minimize surface area, leading to cohesive interactions between fluids and surrounding structures. When the characteristic length scale is sufficiently small, this balance between surface energy and gravity provides a simple physical framework for understanding a wide range of capillarity-driven behaviors.
When solid boundaries interacting with a liquid are soft and deformable, interfacial tension can induce significant elastic deformation, leading to elastocapillary phenomena. A droplet on an elastomer creates ridges at its contact line, thereby modifying its wetting dynamics.5 Similarly, thin elastic filaments can undergo spontaneous twisting driven by the draining of liquids,6 and micro to nanoscale resist structures may collapse during drying due to capillary forces.7 These examples demonstrate that elastocapillary effects not only alter local wetting behavior but can trigger large-scale deformations of soft solids. Such deformation of soft boundaries, in turn, can strongly influence the dynamics of interacting fluids. Beyond engineered systems, similar elastocapillary and capillarity-driven mechanisms are widely exploited in nature, for example, in the water-collecting structures of plants,8–10 the locomotion of water-walking insects,11–13 and emerging soft robotic systems.14,15
Exploiting the coupling between elasticity and capillarity, capillary origami enables the spontaneous formation of three-dimensional structures from two-dimensional elastic films.16,17 By appropriately designing the film shape, such self-folding can be harnessed to capture and lift liquid droplets.18 In their study, the folding behavior and resulting functions were primarily controlled by modifying the film geometry, for example, through the number of petals. Related approaches have also demonstrated elastomer-based devices capable of liquid pipetting driven by capillary forces.19
Despite these advances, the physical criteria governing the onset of liquid encapsulation and self-folding remain unclear, particularly in relation to the competing roles of elasticity, surface tension, and gravity. In this study, we address this issue by systematically varying the side length and thickness of an equilateral triangular film. This minimal geometry allows us to construct a phase diagram for elastocapillary self-folding and liquid encapsulation, which we interpret using characteristic length scales associated with bending, capillarity, and gravity.
In the work by Reis et al.,18 the film geometry was optimized into a petal-like shape to achieve controlled deformation. Nakamura et al.19 employed fixed flower-like geometries and treated the lifting speed and gravity direction as control parameters. In the present study, by contrast, we focused on a minimal equilateral triangular film and observed three deformation morphologies corresponding to folding, recoiling, and liquid encapsulation. Whereas previous studies achieved encapsulation using specifically designed geometries such as petal-shaped films,18 the present study demonstrates that encapsulation can emerge even in a simple triangular geometry through the competition between elastocapillary, elastogravity, and capillary length scales. This provides a more general physical criterion that is not tied to a particular geometric design. We further tested this interpretation by varying the surface tension of the aqueous phase using surfactant solutions.
In the following sections, we describe an experimental system in which an equilateral triangular elastic film is lifted quasi-statically from a water surface while its side length and thickness are systematically varied. Depending on these parameters, the film undergoes folding, recoiling, or liquid encapsulation. We analyze these behaviors in terms of the competition between bending, capillary, and gravitational effects.
:
1. The mixture was spin-coated onto a PTFE substrate at rotation speeds ranging from 500 to 1500 rpm to form thin films. The films were then cured at 100 °C for 1 h to obtain elastic PDMS sheets. After curing, the films were cut into equilateral triangular shapes with a side length w, as shown in Fig. 1(a). The Young's modulus of the PDMS film was taken as E = 0.7 × 106 Pa.20
The liquids used in this study were pure water and a 30 mM aqueous solution of polyethylene glycol mono-4-octylphenyl ether (Triton X). Pure water was prepared using a Millipore Milli-Q system. Triton X was purchased from Tokyo Chemical Industry Co., Ltd (P0873). The density ρ of both liquids was assumed to be 1.0 × 103 kg m−3. The surface tension of water, γw, was taken as 72 mN m−1, whereas that of the Triton X solution, γT, was taken as 30 mN m−1, since the Triton X concentration (30 mM) is well above the critical micelle concentration (CMC ≈ 0.2–0.3 mM). Here, γi denotes the surface tension of the liquid phase, where i = w or T.
The film was initially positioned at the water surface and then lifted vertically at a constant speed of 1 mm s−1. The deformation of the film during lifting was recorded using a CMOS camera (DMK37BUX273, The Imaging Source) at 20 Hz. The force measured by the load cell was recorded using a data logger (GL7000, Graphtec Co.) and used to evaluate the weight of the lifted water, denoted by mg. The value before lifting the film from the water surface was defined as zero, and the increase from this value was plotted as the apparent lifted weight in Fig. 2(b and d). Although the measured force includes the contribution from the film itself, the film weight is at most 0.1 mN (for w = 20 mm and h = 150 µm), which is negligible compared with the measured force. Therefore, the measured force mainly reflects the weight of the water lifted by the film.
![]() | ||
| Fig. 2 (a) and (c) Phase diagrams of the deformation morphologies of the lifted film as functions of film thickness h and side length w for (a) pure water and (c) Triton X solution. Orange triangles, green squares, and blue circles correspond to the folding, recoiling, and encapsulation regimes, respectively [see Fig. 1(c)]. The morphologies were classified based on the final film configuration observed during lifting. The dashed line indicates the parameter set used in panels (b) and (d). All experimental data points are shown in the phase diagrams. (b) and (d) Apparent weight of the liquid lifted by the film, mg, as a function of film thickness h for a fixed side length of w = 10 mm, measured for (b) pure water and (d) Triton X solution. | ||
The liquid encapsulation regime is further characterized by the weight of the lifted liquid, mg, as shown in Fig. 2(b). These data were obtained for films with a fixed side length w = 10 mm. Larger values of mg were observed when the film exhibited the encapsulation morphology. As the system transitions from folding to encapsulation, the lifted liquid weight increases gradually. In contrast, a sharp decrease in mg was observed when the system transitions from encapsulation to recoiling at larger h.
In addition to these dominant behaviors, marginal cases were observed, including intermediate encapsulation–recoiling states, where the deformed film retained a finite liquid droplet at its bottom, as well as asymmetric folding without threefold symmetry. These marginal events are indicated by black crosses in Fig. 2(a) and (b).
We have also conducted additional experiments using a Triton X solution. The corresponding results are shown in Fig. 2(c) and (d). The phase diagram is shifted toward the lower-left corner, and the parameter region for the encapsulation regime becomes smaller. The lifted liquid weight shown in Fig. 2(d) also exhibits a strong dependence on the observed morphology. We also examined the effect of lifting speed and liquid viscosity on the final morphology. The lifting speed was varied from 0.1 to 30 mm s−1. The liquid viscosity was varied from approximately 1 mPa s (pure water) to 50 mPa s using a 30 wt% PEG6000 aqueous solution. However, no substantial change in the final morphology was observed. Under these conditions, the capillary number, Ca = μV/γi, ranged from approximately 10−6 in the standard experimental condition to 10−1, but remained below unity. The Reynolds number Re = ρLV/μ remained of order unity or smaller in the present system. The Weber number We = ρV2L/γi ranged from 10−4 to 10−1, indicating that inertial effects are negligible. These observations indicate that viscous effects associated with the lifting motion do not dominate the morphology selection within the accessible experimental range. We therefore focus on the quasi-static balance between bending, capillary, and gravitational effects.
![]() | ||
| Fig. 3 Characteristic deformation modes of the elastic film, classified as (a) mode I, (b) mode II, and (c) mode III. The dashed lines indicate regions of large mean curvature. Bottom images show the corresponding regions of high mean curvature in the experimental snapshots shown in Fig. 1(c). | ||
In the folding regime, the deformation is initially dominated by mode I [Fig. 3(a)], where ridges form along lines connecting the lifting point to the midpoints of the film edges. Geometrical constraints at the vertices due to capillary force at the air–water interface prevent vertical deformation, leading to the coexistence of mode I and mode III [Fig. 3(c)] at later stages. This combined deformation ultimately results in expulsion of the liquid.
In contrast, recoiling behavior is associated with mode II [Fig. 3(b)], in which the regions near the vertices bend downward. As the film is lifted, capillary breaking occurs, and the elastic energy stored in the film drives the system back to a nearly flat configuration.
Liquid encapsulation emerges in an intermediate regime, where the deformation evolves from a mixed mode I–II configuration toward a mode II–dominated state. In this case, the approach of the three vertices facilitates the formation of a transient liquid bridge, which subsequently undergoes capillary breaking to form a closed capsule.
These deformation modes are not universally traversed in all cases; rather, each observed behavior follows a distinct deformation pathway. The central attachment point enforces a symmetric boundary condition, which stabilizes the threefold deformation patterns but does not qualitatively alter the underlying energy balance. This classification allows us to interpret the phase diagram in terms of transitions between dominant deformation modes.
, where ℓridge ∼ Lec, following the estimation of loop structures in ref. 16 and 21. For the encapsulated state, we estimate Uc ∼ αρgw4 + εB, where the first term represents the gravitational potential energy of the encapsulated liquid volume and the second term represents the bending cost of the capsule. Using simple geometric estimates for these energies, we constructed an energy-minimum phase diagram shown in Fig. 5. Despite the simplified assumptions, the obtained phase structure reproduces the experimentally observed morphology sequence and its shift upon reducing surface tension. These results support the interpretation that liquid encapsulation emerges in the crossover regime where bending, capillary, and gravitational effects become comparable.
![]() | (1) |
When L ≫ Leg, the bending energy becomes subdominant compared to the gravitational energy, and the film tends to undergo folding with significant out-of-plane deformation when the film is lifted. Indeed, Leg is known to correspond to the typical wavelength of the elastic film under the effect of bottom liquids.18 A similar argument applies to the present system. On the other hand, when L ≪ Leg, bending rigidity suppresses large out-of-plane deformation, and the film tends to remain flat and recoil back.
Similarly, when L ≫ Lec, capillary forces are sufficiently strong, and the film undergoes strong capillary-induced bending deformation as reported for elastocapillary loop structures.16 When L ≪ Lec, capillary forces are insufficient to bend the film, resulting in a recoiling state rather than water encapsulation.
When L ≃ Lec and L ≃ Leg, the capillary length Lcg = (γi/ρg)1/2 also becomes comparable to L. Therefore, the crossover regime between these characteristic length scales provides a natural condition for liquid encapsulation. This suggests that encapsulation emerges in the crossover regime where bending, capillary, and gravitational effects compete. We also note the relation
, showing that only two independent characteristic scales exist in the present system.
We note that B = Eh3/{12(1 − ν2)}, where B is proportional to h3. Thus, Leg ∼ h3/4 and Lec ∼ h3/2. Using the parameters listed in Table 1, these characteristic lengths are plotted in Fig. 4(a) and (b). As a function of the film thickness h, Leg increases as h3/4 while Lec increases as h3/2. As a result, the two characteristic lengths become comparable over a finite range of h, irrespective of the precise numerical prefactors.
| Variable | Description | Values |
|---|---|---|
| E | Young's modulus | 0.7 × 106 Pa |
| ν | Poisson's ratio | 0.5 |
| ρ | Density of water | 1.0 × 103 kg m−3 |
| g | Gravitational acceleration | 10 m s−2 |
![]() | ||
| Fig. 4 The characteristic length of the system, elastogravity, Leg (black solid line), elastocapillary, Lec (red dash-dotted line), and capillary length Lcg (green dotted line). Each line is drawn based on eqn (1) with (a) water (γw = 72 mN m−1) and (b) Triton X solution (γT = 30 mN m−1). | ||
![]() | ||
| Fig. 5 Energy-minimum phase diagrams calculated from the simplified energy estimates for (a) water (γw = 72 mN m−1) and (b) Triton X solution (γT = 30 mN m−1). | ||
This crossover regime corresponds to the parameter range where water encapsulation was observed in the experiments (Fig. 4(a) water and (b) Triton X solution), as seen in the phase diagram shown in Fig. 2(a) and (c). We confirmed the overlapping region shifts toward the lower-left corner when Triton X solution (γT = 30 mN m−1) was used, which is in agreement with Fig. 2(a) and (c).
Here, we note that L should be interpreted as an effective lateral deformation scale of the triangular film, rather than the geometric side length w itself. The deformation is localized along curved ridges as illustrated in Fig. 3, and each energy term depends on the film geometry. As a result, the relevant length scale governing bending, capillary, and gravitational effects depends not only on w but also on the detailed deformation geometry and boundary conditions. Therefore, the plot in Fig. 4 captures the ordering and relative location of the experimentally observed morphology regimes, although the comparison should be understood as semi-quantitative because geometry-dependent numerical prefactors are neglected in the present model.
Based on simplified geometric assumptions, we confirmed that the film's encapsulation morphology appears in the crossover regime where the elastogravity, elastocapillary, and capillary length scales become comparable. Despite the simplified estimation of these energies, the obtained phase diagram showed semi-quantitative agreement with the experiment.
In this study, we kept the film's lifting speed constant at 1 mm s−1. Preliminary experiments suggested that dynamical effects become important when the depth of the liquid reservoir is as small as 1 mm and the liquid viscosity is 1 Pa s so that the effective capillary number for the experiment is of the order of 1. Under such conditions, the dynamic effect should be interesting to consider. We leave this regime for future study.
We note that our system is quasi-static in the sense that the observed morphology is essentially independent of lifting speed. This suggests that the system evolves toward a minimum-energy configuration. However, unlike the capillary origami case,16 the lifting process using a rod attached to the center of the triangular film constrains the deformation pathway. Accordingly, a precise understanding of the phase diagram will require not only length-scale and energy-landscape arguments but also consideration of force balance and dynamical effects, even though most of the observed morphology selection is governed by quasi-static energetics. A full elastodynamic treatment would require appropriate handling of film-edge effects and contact-line dynamics. Such an analysis remains an important direction for future study.
Our results show that simple elastic films can be used to manipulate liquids by appropriately choosing the thickness and side length. For potential application, it would be interesting to consider the effect of handling non-Newtonian liquids with a film for future study. Such approaches may enable the transport and manipulation of liquids and soft materials using simple elastic films.
For the recoiling and folding states, the triangular configuration gives
and
. For the encapsulated state, assuming a tetrahedron-like liquid volume with ridge structures [Fig. 3(b): mode II], we estimate the enclosed volume and lifted height as
for the gravitational energy term. For the bending contribution, using the characteristic area and curvature
.
Using these numerical prefactors, we constructed the energy-minimum phase diagrams shown in Fig. 5 for (a) water and (b) Triton X solution.
Footnote |
| † These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2026 |