Previous Article |
Table of Contents
| Next Article
CHEMISTRY
Magnetic self-assembly of three-dimensional surfaces from planar sheets



*Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, MA 02138;
Division of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, MA 02138;
Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027;
Department of Physics, Harvard University, 17 Oxford Street, Cambridge, MA 02138; and ¶School of Computer Science, McGill University, 3480 University Street, Montreal, QB, Canada H3A 2A7
Contributed by George M. Whitesides, February 1, 2005
| Abstract |
|---|
|
|
|---|
This report describes the spontaneous folding of flat elastomeric sheets, patterned with magnetic dipoles, into free-standing, 3D objects that are the topological equivalents of spherical shells. The path of the self-assembly is determined by a competition between mechanical and magnetic interactions. The potential of this strategy for the fabrication of 3D electronic devices is demonstrated by generating a simple electrical circuit surrounding a spherical cavity.
folding | microfabrication | 3D structure | soft lithography | soft electronics
Folding of connected, 2D plates [using robotics (4) or spontaneous folding (59)] can yield 3D microelectromechanical systems (MEMS) and microelectronic devices. We (10, 11) and others (4, 12) have explored a number of routes to small 3D shapes based on self-assembly. These strategies are still early in their development.
Here, we explore a new strategy for formation of 3D objects that combines the advantages of planar microfabrication with those of 3D self-assembly. Our approach comprises four steps (Fig. 1a): (i) cutting the 3D surface of interest into connected sections that "almost" unfold into a plane (unpeeling a sphere as one unpeels an orange is an example); (ii) flattening this surface and projecting it onto a plane; (iii) fabricating the planar projection in the form of an elastomeric membrane patterned with magnetic dipoles; and (iv) allowing this patterned membrane to fold into an "almost-correct" 3D shape by self-assembly. This strategy offers the potential to transform easily patterned, functionalized planar sheets into 3D structures and devices. It also raises the problem of designing and generating stable 3D structures by decomposing and projecting these structures into 2D shapes and then balancing the shapes of 2D cuts, the placement of magnetic dipoles, and the mechanical characteristics of the membrane.
|
Folding a planar sheet into a surface that is curved locally along two orthogonal axes causes the middle surface of the sheet to be stretched because of changes in the Gaussian curvature (22). Because stretching a sheet is energetically more expensive than bending it, we minimized the geometric incompatibilities between the planar and the 3D forms by using appropriately shaped, 2D, soft, elastomeric membranes. We chose magnetic forces to guide the self-assembly for three reasons: (i) magnetic interactions are insensitive to the surrounding medium and to the details of surface chemistry (8); (ii) the distances over which they act can be engineered to cover a range of sizes (nanometers to meters); and (iii) magnetic dipoles tend to form stable closed loops, and these loops are features easily translated into design rules (23).
| Design of 2D Elastomeric Sheets Patterned with Magnetic Dipoles |
|---|
|
|
|---|
We explored elastomeric sheets having three different 2D shapes: an equatorial cut (Fig. 1bi), an orange-peel cut (Fig. 1bii), and a flower-petal cut (Fig. 1biii). We generated the shapes of the flat sheets by using a modified version of the software package UNFOLD POLYTOPE for MATHEMATICA (www.cs.mcgill.ca/~fukuda/download/mathematica). In the original package, a polytope (specified as a collection of faces and edges) is unfolded by cutting a number of edges and rotating the faces around the remaining edges until they all lie in the same plane. Fig. 2 illustrates the unfolding process for one vertex, D. The new location, D', of the vertex, D, is calculated as
![]() |
|
|
|
| Results and Discussion |
|---|
|
|
|---|
Finite-Element Simulation of Magnetic Field Profiles. We performed a finite-element simulation of the 2D magnetic field profiles of the unfolded flat precursors and of cross sections of the folded 3D structures by using the software package FINITE ELEMENT METHOD MAGNETICS (http://femm.foster-miller.com). We set the magnetic field intensity at the surface of each magnet to 0.04 T (the experimentally measured upper limit) and treated the poly(dimethylsiloxane) (PDMS) membrane as a diamagnetic material with the magnetic characteristics of air (a preset feature of the software). From an input of the shape of the flat sheet (or a cross section of the 3D structure) patterned with magnetic features of given shape, polarity, and magnetic permeability, this software calculates and maps the magnetic field intensity and the magnetic field lines. In the output of the simulation, the intensity of the magnetic field is depicted by color intensity: Higher field intensity corresponds to a darker shade of gray. The strength of the interactions between the patterned magnets can be estimated from the strength of the magnetic field and the number of field lines in the space surrounding the magnets. Thus, a magnetic energy minimum corresponds to a map in which the field intensity and the field lines are confined only to the regions of high magnetic permeability (i.e., within the magnets) and do not extend (or extend only to a negligible degree) into the regions of low magnetic permeability (i.e., in the plane of the elastomeric membrane, outside of the magnets).
Fig. 5 shows the results of the simulation. In the unfolded equatorial cut (Fig. 5a), there were significant interactions between adjacent magnets positioned in the middle of the segments. No interactions were observed between the magnets positioned at the tips of the segments and between magnets positioned in the middle and at the tips of each segment. In the unfolded flower-petal cut (Fig. 5b), there was appreciable interaction between the magnets positioned in the middle of the segments. The central region containing six triangular magnets formed a completely closed ring of magnetic dipoles and showed no stray field lines extending away from the magnets. There was no interaction between magnets positioned at the tips of adjacent sections of the membrane; there was also no interaction between those magnets and the magnets positioned in the middle of the segments. Upon folding of both the equatorial cut and the flower-petal cut, the magnets positioned at the tips of the segments (Fig. 5e) and the magnets positioned in the middle of the segments (Fig. 5f) formed continuous, closed rings of magnetic dipoles; the field intensity and the field lines in these areas were confined exclusively within the rings, indicating minimized magnetic energy of the folded structures.
|
The flower-petal cut patterned with magnets arranged anti-parallel to one another (Fig. 5d) did not form a sphere; it folded, instead, into stable structures comprising groups of two or four petals (see Fig. 4h). Fig. 5 h and i shows simulations of an aggregate comprising three pairs of magnetized features and an aggregate comprising one pair and a group of four magnetized features, respectively. These aggregates correspond to cross sections of the polar regions of a sphere formed from the flower-petal cut shown in Fig. 5d, in which the groups of two or four magnetized features are separated by a gap of 400 µm. The simulation shows that the field intensity in these gaps is very low or equal to zero; the dipoles aggregated in pairs or in a group of four have formed closed, (magnetically) stable loops. The energy gain of bringing these small loops into a single loop comprising all six segments would be insignificant compared with the energy gain in the case of a flower-petal cut patterned with equatorially magnetized features (as the one shown in Fig. 5e).
Estimation of the Elastic and Magnetic Free Energy of the Folded Structures. The process of folding minimizes the sum of the elastic and magnetic energies of the system: The increase in mechanical energy that occurs on folding the elastic sheet is balanced by the decrease in magnetic energy as the magnetic dipoles approach one another. We used the balance of the two energies to estimate how changing the size of the components and the thickness of the elastic membrane affects the stability of folded structures. We compared the magnetic energy gained in bringing non-interacting magnetic dipoles into the equatorial loop of the structure shown in Fig. 4 a and b with the elastic energy cost for bending the sheet around that loop.
To estimate the elastic energy of the system, as a first approximation we considered only the bending energy, Ub, needed to bend a sheet of thickness h made of a material with Young's modulus E into a cylinder by using magnets of length R placed in the middle of the segments. Thin elastic sheets are most easily folded by isometric bending of the middle surface, without stretching it (24); here, we ignored the contribution of the stretching energy that arises because of the presence of boundary layers at ridges (25, 26) and peaks (21, 27). Thus, the elastic energy scales as
![]() | [1] |
where
1/R is the local mean curvature at the center of the sheet and A
R2 is its area. The elastic energy due to bending alone is largely independent of the radius of the sphere.
To estimate the magnetic energy of the system, we approximated the magnetic field generated by permanent magnets of length L (proportional to the radius of the sphere) and square cross section b with the field generated by a finite, current-carrying solenoid (28). The magnetic field, Bisolated, inside each isolated magnet in the unfolded configuration can be described as the magnetic field along the axis of the finite solenoid:
![]() | [2] |
where M0 is the magnetic strength of the material, z is the distance from the point at which we measure the field to the face of the solenoid, L is the length of the solenoid, and b is its radius (proportional to the width of the magnets). To estimate the magnetic field, Bring, of the folded configuration, we made the approximation that the magnetic field inside the ring of magnets is constant (similar to the field inside a solenoid torus); thus, it can be expressed (28) as
![]() | [3] |
Using Eqs. 2 and 3, we can express the change of the magnetic energy, Um, between the unfolded and the folded configurations as
![]() | [4] |
where the integration is done over the volume, V, of the magnets. We calculated the magnetic energy difference due to folding numerically, for the limiting case of L » b, as
![]() | [5] |
This rough estimation of the magnetic energy ignores the shape of the closed ring of magnets (in the structures considered in this work, the ring is shaped as a polygon and not as a circle), the variation of the magnetic field away from the main axis of the solenoid, and the magnetic energy outside the volume of the magnets. Nevertheless, the scaling of the magnetic energy with the dimensions of the magnets is largely independent of these effects.
Using Eqs. 1 and 5, we can then express the balance of the magnetic and elastic energies as
![]() | [6] |
It is readily apparent that the ratio Um/Ub depends on the ratios
, b/h, and L/b. The ratio
describes properties of the materials used (magnets and elastic sheets) and does not depend on the dimensions of the folded structure. Using the ratios b/h and L/b, we can estimate how changing the dimensions of the components and the thickness of the elastic membrane will affect the stability of a folded structure: If the width and the length of the magnets are rescaled by a factor
, the thickness of the elastic sheet must be rescaled by the same factor to preserve the balance of the two energies and, thereby, the stability of the folded structure.
An important factor that this simple calculation does not take into account is the possibility of coexistence of several stable, folded shapes for a given set of parameters. The next step in refining the scaling estimate will include studies of the energy landscape of the folded structures and the kinetics of the folding process.
Three-Dimensional Electrical Circuit. We used this strategy to generate elementary 3D electrical circuit (Fig. 6). We fabricated each section of a sheet having the design shown in Fig. 6a to include an electrically isolated wire (see supporting information for fabrication details). This wire connected the cathode and the anode of a LED to two pairs of solder pads placed at the tips of each section, near the edge of the membrane; these solder pads were intended to provide electrical connectivity between LEDs placed on adjacent sections of the sheet. The solder features (wires and contact pads) were embedded within the PDMS membrane, with the wires close to the middle surface of the sheet that is not stretched when rolled into a cylindrical surface (29, 30) (Fig. 6b). When cutting the sheet to shape, we cut through the outside edges of the embedded solder pads and, thus, produced four edges exposing bare solder adjacent to the magnets at the tips of each section. For self-assembly, we suspended the sheet in water at 60°C (the melting temperature of the solder was 47°C) and agitated gently. Within 13 min, the planar sheet folded into a sphere; simultaneously, electrical connections formed between the six sections by fusion of the drops of molten solder (Fig. 6c). The heat also restored the electrical continuity of wires that had broken during bending of the flat sheet. After self-assembly was complete, and the resulting structure had cooled to room temperature and dried, we connected one top and one bottom contact pad to a battery. Fig. 6d shows that all six LEDs illuminate and demonstrates the continuity of an electrical circuit that traces a path through all six sections of the structure.
|
| Conclusions |
|---|
|
|
|---|
This work also identifies a remarkably interesting problem in applied mathematics: the analytical design of planar sheets that can fold into 3D structures with minimized global (magnetic and mechanical) free energy. Answering the core questionhow to translate a global goal (designing and generating a stable 3D structure) into local interactions defined by the shape of the 2D cut, the placement of the magnetic dipoles, and the mechanical characteristics of the membranewill require further development in both fabrication and applied mathematics.
This approach combines the efficient methods of photo- and soft lithography for fabrication of flat, patterned surfaces with self-assembly to form 3D from 2D structures. It can be generalized to other interactions (e.g., interactions between patterned electrical monopoles or dipoles, and capillarity) and to a broad range of materials. This strategy is relevant to the fabrication of "soft" electronics (31) (e.g., displays, thin film transistors, and sensor skins) and deployable structures (32). The most immediate problem that should be solved in moving toward applications is to improve the method of fabrication of the precursor membranes. In principle, both fabricating the sheet (with embedded electrical components) and patterning the magnetic dipoles can be accomplished by using planar microfabrication (33, 34).
| Acknowledgements |
|---|
| Footnotes |
|---|
Abbreviations: LED, light-emitting diode; PDMS, poly(dimethylsiloxane).
|| To whom correspondence should be addressed. E-mail: gwhitesides{at}gmwgroup.harvard.edu.
© 2005 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
A. W. Feinberg, A. Feigel, S. S. Shevkoplyas, S. Sheehy, G. M. Whitesides, and K. K. Parker Muscular Thin Films for Building Actuators and Powering Devices Science, September 7, 2007; 317(5843): 1366 - 1370. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||