Benny Davidovitch

Prototypical model for tensional wrinkling in thin sheets

The buckling and wrinkling of thin films has recently seen a surge of interest among physicists, biologists, mathematicians, and engineers. This activity has been triggered by the growing interest in developing technologies at ever-decreasing scales and the resulting necessity to control the mechanics of tiny structures, as well as by the realization that morphogenetic processes, such as the tissue-shaping instabilities occurring in animal epithelia or plant leaves, often emerge from mechanical instabilities of cell sheets. Although the most basic buckling instability of uniaxially compressed plates was understood by Euler more than two centuries ago, recent experiments on nanometrically thin (ultrathin) films have shown significant deviations from predictions of standard buckling theory. Motivated by this puzzle, we introduce here a theoretical model that allows for a systematic analysis of wrinkling in sheets far from their instability threshold. We focus on the simplest extension of Euler buckling that exhibits wrinkles of finite length—a sheet under axisymmetric tensile loads. The first study of this geometry, which is attributed to Lamé, allows us to construct a phase diagram that demonstrates the dramatic variation of wrinkling patterns from near-threshold to far-from-threshold conditions. Theoretical arguments and comparison to experiments show that the thinner the sheet is, the smaller is the compressive load above which the far-from-threshold regime emerges. This observation emphasizes the relevance of our analysis for nanomechanics applications.

Curvature-induced stiffness and the spatial variation of wavelength in wrinkled sheets

Significance Thin elastic sheets buckle and wrinkle to relax compressive stresses. Wrinkling metrologies have recently been developed as noninvasive probes of mechanical environment or film properties, for instance in biological tissues or textiles. This work proposes and experimentally tests a prediction for the local wavelength of wrinkles in nonuniform curved topographies. Wrinkle patterns in compressed thin sheets are ubiquitous in nature and technology, from the furrows on our foreheads to crinkly plant leaves, from ripples on plastic-wrapped objects to the protein film on milk. The current understanding of an elementary descriptor of wrinkles—their wavelength—is restricted to deformations that are parallel, spatially uniform, and nearly planar. However, most naturally occurring wrinkles do not satisfy these stipulations. Here we present a scheme that quantitatively explains the wrinkle wavelength beyond such idealized situations. We propose a local law that incorporates both mechanical and geometrical effects on the spatial variation of wrinkle wavelength. Our experiments on thin polymer films provide strong evidence for its validity. Understanding how wavelength depends on the properties of the sheet and the underlying liquid or elastic subphase is crucial for applications where wrinkles are used to sculpt surface topography, to measure properties of the sheet, or to infer forces applied to a film.

Indentation of ultrathin elastic films and the emergence of asymptotic isometry.

We study the indentation of a thin elastic film floating at the surface of a liquid. We focus on the onset of radial wrinkles at a threshold indentation depth and the evolution of the wrinkle pattern as indentation progresses far beyond this threshold. Comparison between experiments on thin polymer films and theoretical calculations shows that the system very quickly reaches the far from threshold regime, in which wrinkles lead to the relaxation of azimuthal compression. Furthermore, when the indentation depth is sufficiently large that the wrinkles cover most of the film, we recognize a novel mechanical response in which the work of indentation is transmitted almost solely to the liquid, rather than to the floating film. We attribute this unique response to a nontrivial isometry attained by the deformed film, and we discuss the scaling laws and the relevance of similar isometries to other systems in which a confined sheet is subjected to weak tensile loads.

Elastic sheet on a liquid drop reveals wrinkling and crumpling as distinct symmetry-breaking instabilities.

Smooth wrinkles and sharply crumpled regions are familiar motifs in biological or synthetic sheets, such as rapidly growing plant leaves and crushed foils. Previous studies have addressed both morphological types, but the generic route whereby a featureless sheet develops a complex shape remains elusive. Here we show that this route proceeds through an unusual sequence of distinct symmetry-breaking instabilities. The object of our study is an ultrathin circular sheet stretched over a liquid drop. As the curvature is gradually increased, the surface tension stretching the sheet over the drop causes compression along circles of latitude. The compression is relieved first by a transition into a wrinkle pattern, and then into a crumpled state via a continuous transition. Our data provide conclusive evidence that wrinkle patterns in highly bendable sheets are not described by classical buckling methods, but rather by a theory which assumes that wrinkles completely relax the compressive stress. With this understanding we recognize the observed sequence of transitions as distinct symmetry breakings of the shape and the stress field. The axial symmetry of the shape is broken upon wrinkling but the underlying stress field preserves this symmetry. Thus, the wrinkle-to-crumple transition marks symmetry-breaking of the stress in highly bendable sheets. By contrast, other instabilities of sheets, such as blistering and cracking, break the homogeneity of shape and stress simultaneously. The onset of crumpling occurs when the wrinkle pattern grows to half the sheet’s radius, suggesting a geometric, material-independent origin for this transition.

On the Stabilization of Ion Sputtered Surfaces

The classical theory of ion beam sputtering predicts the instability of a flat surface to uniform ion irradiation at any incidence angle. We relax the assumption of the classical theory that the average surface erosion rate is determined by a Gaussian response function representing the effect of the collision cascade, and consider the surface dynamics for other physically motivated response functions. We show that although instability of flat surfaces at any beam angle results from all Gaussian and a wide class of non-Gaussian erosive response functions, there exist classes of modifications to the response that can have a dramatic effect. In contrast to the classical theory, these types of response render the flat surface linearly stable, while imperceptibly modifying the predicted sputter yield vs incidence angle. We discuss the possibility that such corrections underlie recent reports of a “window of stability” of ion-bombarded surfaces at a range of beam angles for certain ion and surface types, and describe some characteristic aspects of pattern evolution near the transition from unstable to stable dynamics. We point out that careful analysis of the transition regime may provide valuable tests for the consistency of any theory of pattern formation on ion sputtered surfaces.

Spreading of viscous fluid drops on a solid substrate assisted by thermal fluctuations

We study the spreading of viscous drops on a solid substrate, taking into account the effects of thermal fluctuations in the fluid momentum. A nonlinear stochastic lubrication equation is derived, and studied using numerical simulations and scaling analysis. We show that asymptotically spreading drops admit self-similar shapes, whose average radii can increase at rates much faster than these predicted by Tanners law. We discuss the physical realizability of our results for thin molecular and complex fluid films, and predict that such phenomenon can in principal be observed in various flow geometries.

Capillary interactions among spherical particles at curved liquid interfaces

We study the effect of interfacial curvature on the binding energy and forces exerted on small spherical particles that adsorb on an interface between two immiscible liquids. When the interface has anisotropic curvature, the constant-contact-angle condition at the particle-fluid boundary requires a deformation of the interface. Focusing on the case of an initially cylindrical interface, we predict the shape after a spherical particle binds. We then calculate the energy of adsorption and find that it depends on the shape of the interface very far from the binding site. Turning to the problem of two adsorbed spherical particles, we predict a capillary interaction that arises purely from the deformations caused by the contact-angle condition. An analogy is made between these curvature-induced capillary forces and electrostatic forces between quadrupoles in two dimensions. We conclude with a conjectured general form for the interaction of a single spherical particle with the Gaussian curvature of the underlying fluid interface, which we compare to previous work.

Elastic building blocks for confined sheets

We study the behavior of thin elastic sheets that are bent and strained under a weak, smooth confinement. We show that the emerging shapes exhibit the coexistence of two types of domains. A focused-stress patch is subject to a geometric, piecewise-inextensibility constraint, whereas a diffuse-stress region is characterized by a mechanical constraint-the dominance of a single component of the stress tensor. We discuss the implications of our findings for the analysis of elastic sheets under various types of forcing.

Dynamics of conformal maps for a class of non-Laplacian growth phenomena.

Time-dependent conformal maps are used to model a class of growth phenomena limited by coupled non-Laplacian transport processes, such as nonlinear diffusion, advection, and electromigration. Both continuous and stochastic dynamics are described by generalizing conformal-mapping techniques for viscous fingering and diffusion-limited aggregation, respectively. The theory is applied to simulations of advection-diffusion-limited aggregation in a background potential flow. A universal crossover in morphology is observed from diffusion-limited to advection-limited fractal patterns with an associated crossover in the growth rate, controlled by a time-dependent effective Péclet number. Remarkably, the fractal dimension is not affected by advection, in spite of dramatic increases in anisotropy and growth rate, due to the persistence of diffusion limitation at small scales.

Universal collapse of stress and wrinkle-to-scar transition in spherically confined crystalline sheets

Imposing curvature on crystalline sheets, such as 2D packings of colloids or proteins, or covalently bonded graphene leads to distinct types of structural instabilities. The first type involves the proliferation of localized defects that disrupt the crystalline order without affecting the imposed shape, whereas the second type consists of elastic modes, such as wrinkles and crumples, which deform the shape and also are common in amorphous polymer sheets. Here, we propose a profound link between these types of patterns, encapsulated in a universal, compression-free stress field, which is determined solely by the macroscale confining conditions. This “stress universality” principle and a few of its immediate consequences are borne out by studying a circular crystalline patch bound to a deformable spherical substrate, in which the two distinct patterns become, respectively, radial chains of dislocations (called “scars”) and radial wrinkles. The simplicity of this set-up allows us to characterize the morphologies and evaluate the energies of both patterns, from which we construct a phase diagram that predicts a wrinkle–scar transition in confined crystalline sheets at a critical value of the substrate stiffness. The construction of a unified theoretical framework that bridges inelastic crystalline defects and elastic deformations opens unique research directions. Beyond the potential use of this concept for finding energy-optimizing packings in curved topographies, the possibility of transforming defects into shape deformations that retain the crystalline structure may be valuable for a broad range of material applications, such as manipulations of graphene’s electronic structure.