Abstract
Hygroscopic soft actuators offer an attractive means to convert environmental energy to mechanical motions as they use water vapor, a ubiquitous substance in the atmosphere. To overcome the limits of existing hygroactuators, such as simplistic actuation mode, slow response, and low efficiency, here we present three kinds of humidity-powered soft machines adopting directionally electrospun hygroresponsive nanofibrous sheets. The wheels, seesaws, and vehicles developed in this work utilize spatial humidity gradient naturally established near moist surfaces such as human skin, so that they operate spontaneously, realizing energy scavenging or harvesting. We also constructed a theoretical framework to mechanically analyze their dynamics, which allowed us to optimize their design to obtain the highest motion speed physically possible.
Introduction
Transport of water in soft media can induce volume change and consequently power movements in artificial (e.g., hydrogels) and natural (e.g., plants) actuators. In particular, water vapor, which is ubiquitous in the atmosphere, is utilized by many plants to drive self-propulsion,1–3 self-sealing, 4 and seed dispersal.5,6 Those natural actuators commonly adopt a bilayer configuration, where one layer is hygroscopically active, while the other layer is inactive. For example, seed awns of wild wheat, Erodium and Pelargonium species, exploit hygroresponsivity of polymers such as lignin and hemicellulose in the walls of cells that constitute the active layer. On the contrary, the walls of cells serving as an inactive layer are tightly bound by inextensible cellulose microfibrils.7–10
Although the hygroscopic actuators in nature have served their functions for survival and reproduction over millions of years,11,12 several issues remain to be solved before the biological principle can be translated to practical applications including soft robotics and energy harvesters. First, most of the botanical actuators rely on changes of humidity over a long timescale as days or even seasons. Thus, the spatial, as well as temporal, gradient of humidity needs to be exploited as a practically reliable energy source, which should drive actuators at a feasible level of frequency (of the order of 0.1–1 Hz). Second, various modes of actuation beyond what are exhibited by plants should be devised using the spatial gradient of humidity. Third, the design of actuation system should be optimized to maximize the efficiency.
Spatial gradient of humidity is naturally established near moist surfaces, including human skin, wet soil, shores, and kitchen and laboratory benches, where the high humidity level drops down to the atmospheric dry level over a distance on the order of 1 cm. 13 To exploit the spatial gradient, we need to either make hygroscopic actuators continually move between the regions of high and low humidity or devise a way to continually change the humidity level around the actuators. The former scheme was used by a so-called hygrobot, whose bilayer actuator crosses the humidity boundary layer by bending–unbending cycles. 14
The latter scheme was realized by deploying a great number of actuators in parallel to scale up the actuation force to drive the mechanical system. 15 Scaling up through increasing the number of actuators, however, grows the size and mass of the machines, which may prevent them from utilizing humidity gradient established only within ∼1 cm from natural moist surfaces. Despite those early attempts, diverse design options to exploit the naturally established spatial humidity gradient are wanted to cope with a range of practical situations.
Just as most of the mechanical systems that harvest or scavenge environmental energy, the power conversion efficiency and energy density play a critical role for the success of hygroresponsive actuators. Although many previous researches on soft actuation with environmental humidity focused on the development of novel materials16–21 and elucidation of mechanical principles, 7 design optimization through rational analysis has been rare.22–25 Seeking the optimized machine design is important because it enables not only a dramatic improvement in energy conversion performance, but also exploitation of a naturally occurring humidity gradient on wet surfaces without unnecessary enlargement of the system. 14
To address these issues, here we introduce three types of machines—hygrowheels, hygroseesaws, and hygrovehicles—that are spontaneously driven by spatial humidity gradient without controlled external stimuli. We mechanically analyze their dynamics, and optimize their design to ensure sufficient force to drive centimetric machines. We start with explaining our bilayer actuator, an engine that drives our machines. We term all of these machines hygromachines, for they commonly convert the chemical potential energy of environmental humidity to mechanical work.
Materials and Methods
Fabrication of hygrowheels
The rim and struts of hygrowheel were made of 100-μm-thick polyvinyl chloride (PVC) film, which was cut by a film cutting machine (Silhouette portrait). To prevent bilayers from directly contacting the substrate, we attached polyethylene terephthalate films to the free end of bilayers. The wet substrate was realized by damped filter paper (GB005 Blotting Paper; Whatman). The width of bilayer was 8 mm and the inactive layer comprises a 35-μm-thick polyimide (PI) film and a 25-μm-thick silicone adhesive.
Fabrication of hygroseesaws
The transparent rod of hygroseeaw was made of 200-μm-thick PVC film. The central axis and fixed rods were made of acrylic shaft of 2-mm diameter. The sliders were made of white V4 resin using a stereolithography (SLA)-type 3D printer (Formlabs; Form 2). To reduce friction of sliders, we lubricated the contact surface between machine parts with a lubricant (WD-40 Multi-Use Product). The electromagnetic generator consisted of a circular magnet plate coiled with copper wire. The magnet used was neodymium magnet 10 mm in diameter and 3 mm in height, which has a shaft at the center. Copper coils were formed by winding a total of 300 m of 40-gauge magnet wire.
To connect the axes of seesaw and generator, we used plastic gear and pinion, which have 40 and 8 teeth, respectively. The width of bilayer was 24 mm and the inactive layer comprises a 50-μm-thick PI film and a 30-μm-thick silicone adhesive.
Fabrication of hygrovehicles
The transparent chamber, shutter, and wheel of hygrovehicle were made of 200-μm-thick PVC films. The gears, axles, and pins were made of white V4 resin using an SLA-type 3D printer (Formlabs; Form 2). To reduce friction of shutter and cylindrical gears, we lubricated the contact surface with a lubricant (WD-40 Multi-Use Product). Two pairs of magnets were made of magnet paper (Canon magnetic photographic paper MG-101). The wet floor was realized by damped filter paper (GB005 Blotting Paper; Whatman). The radius of wheel was 15 mm and the distance between the bottom window and wet floor was maintained at 2 mm. The width of bilayer was 32 mm and the inactive layer comprises a 50-μm-thick PI film and a 30-μm-thick silicone adhesive.
Results
Hygroengines: Hygroresponsive bilayer actuators
Machines powered by environmental humidity variation need an engine to convert the surrounding chemical potential energy into the mechanical work, and a mechanical structure to translate the cyclic motion of engine into a desired motion. We adopt a phytomimetic bilayer actuator as an engine of our machines. The active layer that expands via absorbing environmental humidity is fabricated by electrospinning nanofibers of polyethylene oxide (PEO). PEO was selected as the humidity-responsive material due to its superior hygroexpansivity, low toxicity, and ease of electrospinning. 26 As shown in Figure 1A, a 10 wt% aqueous PEO solution is pumped through multiple metallic nozzles, at whose tips nanojets are ejected onto a rotating drum collector under a high electrical field (1.5 kV/mm). By tuning the spooling rate to the linear speed of the jets, we get a sheet of directionally aligned nanofibers, whose scanning electron microscopy image is shown in Figure 1B.

Hygroengine made of a hygroscopically responsive bilayer.
The aligned microstructure in the active layer, as inspired by seed awns of Erodium 1 and Pelargonium 2 species, is known to enhance both the actuation speed and magnitude. 14 By moving the nozzles parallel to the drum axis using a linear stage, as shown in Figure 1A, we can obtain a wide fibrous sheet (5 cm in maximum width). The active layer thickness can be varied by adjusting the jetting duration. The empirically measured properties of active PEO layer are listed in Supplementary Table S1. A humidity-insensitive PI film is attached to the hygroscopically responsive PEO layer to result in a bilayer structure (Fig. 1B). The bonding of the active and inactive layers is achieved via a 25-μm-thick silicone adhesive film having an adhesive shear strength of 0.5 MPa.
When a bilayer is exposed to humid air, the water vapor diffuses into the active layer and the water concentration ϕ(z,t) follows the one-dimensional unsteady diffusion equation: ∂ϕ/∂t = D∂ 2 ϕ/∂z 2 with D being the vapor diffusivity in the porous fibrous layer. In our previous research, we empirically determined the water vapor diffusivity in the directionally electrospun PEO sheet. 14 With the initial uniform water concentration ϕ0, the initial condition is given by ϕ(z,0) = ϕ0. Under the atmospheric humidity of ϕ∞, the boundary conditions are given by ϕ(0,t) = ϕ∞ and ∂ϕ(h,t)/∂z = 0, where we considered no moisture flux at the interface of active and inactive layers (z = h). As the bilayers situated in our hygromachines do not cross the humidity boundary layer, we assume that the humidity surrounding the active layer is invariant (Supplementary Text S6 and Supplementary Fig. S5).
Solving the equation gives the spatiotemporal evolution of the humidity in the active layer, which allows us to find the strain tensor:
The equilibrium equations for force and moment lead us to predict the bilayer curvature as a function of time t and external loads, the detailed procedure explained in Supplementary Text S1. Figure 1C shows the experimental images of a bilayer actuator clamped at one end, which is initially uncurved as fabricated at ϕ∞ = 0.4, bends downward when exposed to humid air (ϕ∞ = 0.8), and upward when dry (ϕ∞ = 0.2).
Hygrowheels
A hygrowheel is a circular rim to which multiple bilayer actuators and their guiding struts are attached, as shown in Figure 2A. The wheel is placed on a filter paper under which is a bath of 40°C warm water. As shown in the upper right panel of Figure 2A, the high humidity near the substrate drops down to the atmospheric value over ∼2 cm, as measured by a humidity sensor (Sensirion SHT85). The bilayers are curved close to the rim when dry, but bend away from the rim when wet. The rolling sequence is depicted in Figure 2B, which starts when a strut A entering the humidity boundary layer touches the substrate. The bilayer actuator at B bends to lift and rotate the entire wheel with point A as a pivot.

Hygrowheels.
When the center of mass of the wheel crosses the vertical pivot line (blue dashed line), a unit step of rolling is completed. Because dry bilayers keep entering the humid region with wheel rolling, our hygrowheels can be continually driven by the spatial humidity gradient (Supplementary Video S1). The wheel can even run uphill as long as the bilayer bending can produce sufficient force to rotate the wheel beyond the pivot line, as shown in Figure 2C.
We can calculate the bending rate of bilayer within the humidity boundary layer and the time taken by the wheel's center of mass to cross the pivot line, τw, as discussed in Supplementary Text S2 and Supplementary Figure S1. Since the wheel rotates an angle of αw = 2π/n, when the pivot line is crossed, the linear translation speed of the wheel is given by Uw = 2lwsin(αw/2)/τw, where n is the number of struts and lw the distance of pivot from the center of wheel. We can find the maximum tilt angle of substrate the wheel can climb (Supplementary Text S2 and Supplementary Fig. S1). The wheel speed is a function of the substrate tilt angle and the geometry of bilayer and wheel. We plot the dependence of wheel translational speed on the thickness of active layer and the tilt angle in Figure 2D.
We compare our theoretical predictions with the experimental data at a fixed active layer thickness and at a zero tilt angle in Figure 2E and F, respectively, to verify our model. Figure 2E shows that the speed decreases with the increase of tilt angle because the horizontal distance from the wheel center to the pivot line increases. Figure 2F shows that there exists an optimal thickness of active layer that maximizes the wheel speed. If the active layer is too thin (thinner than 17.4 μm as calculated in this plot), it cannot generate enough force to push the wheel to rotate. Too thick an active layer needs too long a time for vapor to diffuse into the soft medium, slowing its response.
The theoretical maximum speed (26 mm/s) at 28 μm thickness is 4.4 times higher than at 55 μm thickness, clearly demonstrating the importance of mechanical analysis for efficient design of hygromachines.
Hygroseesaws
A hygroseesaw is a rod with the axis of rotation at the center and two hygropistons loaded at two ends, as depicted in Figure 3A. The hygropiston consists of a bilayer actuator, one end of which is fixed to the rod tip, while the other end is attached to a cylinder of mass ms capable of sliding along a slot in the rod. The hygropiston pushes the cylindrical slider as the bilayer is straightened out under high humidity, but pulls the slider as the bilayer bends with drying. Figure 3B schematically shows the working principle of a hygroseesaw along with the corresponding experimental images.

Hygroseesaws.
In the experiment, a transparent chamber houses vapor supplied from the underlying water bath. At stage 1, the hygropiston at the left end, which has just entered the humid chamber, pushes the slider closer to the center of rod. At the right end in the low humidity region, the hygropiston pulls the slider away from the center of rod. When the mass of a bilayer is mb and the chord length difference between two bilayers is Δl, the resulting torque imbalance is calculated as ΔT = (ms + mb/2)gcos θΔl with θ being the tilt angle of the seesaw, and g the gravitational acceleration. Because of intrinsic friction at the rotation axis, ΔT should exceed a critical value ΔTc determined by the friction at the rotation axis and the seesaw weight, to rotate the seesaw.
Through an intermediate rotation stage 2, which takes much shorter than stage 1, the seesaw enters stage 3 where the hygropiston at the right end is under high humidity while it is dry at the left end. Now the right piston pushes and the left piston pulls their own slider, causing the seesaw to rotate in the opposite sense in stage 4 and to return to stage 1. In this manner, the seesaw can continually rotate reciprocally as driven by the spatial gradient of humidity (Supplementary Video S2).
By computing the sliding distances of both the hygropistons as a function of time, we can find the time taken for ΔT to exceed ΔTc, which corresponds to a half period of seesawing cycle (Supplementary Text S3 and Supplementary Fig. S2). The period is a function of the active layer thickness ha and the mass of slider ms, once the bilayer length, its inactive layer thickness, and the rod length are given. We plot the seesawing frequency as a function of ha and ms in Figure 3C. We compare our theoretical predictions with the experimental data at fixed ha and ms in Figure 3D and E, respectively, to verify our model. Figure 3D shows that the frequency increases with the slider mass because larger ms facilitates the tipping of seesaw by breaking the torque balance more easily.
However, ms cannot be increased indefinitely because too heavy a slider overstretches or overcurls the hygropiston due to gravitational effects, which was observed when ms exceeded 1 g in our experiments. Figure 3E shows that there exists an optimal thickness of active layer that maximizes the seesawing frequency. Just as in the hygrowheel, too thin an active layer cannot generate enough force to move the hygropiston overcoming the slider weight, and too thick a layer needs excessively long time for vapor to diffuse into the medium. The theoretical maximum frequency (0.076 Hz) at 56 μm thickness is 1.8 times higher than at 100 μm thickness.
Our hygroseesaw can continually convert the chemical potential energy associated with spatially varying environmental humidity into electrical energy by connecting a simple electromagnetic generator, as shown in Figure 3F and Supplementary Video S3. The axes of seesaw and generator are linked with gears of 5:1 ratio to increase the rotating angle and angular velocity of a magnet in the generator. We used different electric circuits to measure the instantaneous voltage induced on a load resistor, to rectify the alternating current to direct current, and to power two light emitting diodes (LEDs), with details given in Supplementary Figure S3.
The peak voltage across a load resistance over 10 kΩ was measured to be 3.0 V as shown in Figure 3G, but the maximum power was measured to be 2.5 mW with a load resistance of 500 Ω. The rectified output had a duration of 2 s, Figure 3H. The LEDs blinking according to the seesaw rotation are indicated in Figure 3F.
Hygrovehicles
A hygrovehicle consists of two pairs of wheels that carry an engine bay between them, as schematically shown in Figure 4A. The linear reciprocal motion of hygropiston within the bay is converted to rotation of an axle via a one-way rotating gear. Two pins moving with the hygropiston slide back and forth on spiral grooves of cylindrical gears, causing the gears to rotate in one way (Supplementary Video S4). To harness the spatial gradient of humidity while keeping the hygroresponsive bilayer's relative vertical distance from the underlying moist surface, the humidity around the hygropiston is periodically varied by using a hygrochamber, which is operated as depicted in Figure 4B and Supplementary Video S5. The chamber has a shutter that can simultaneously open the top window and close the bottom window, or vice versa, when a critical force from the bilayer is applied.

Hygrovehicles.
Stage 1 starts as soon as the shutter has opened the bottom window and closed the top window. Then the humid air enters the chamber through the bottom window, increasing the interior humidity. Although the active layer in the bent bilayer absorbs water vapor in this stage, the bilayer cannot unbend until its force on the shutter gets strong enough to overcome the resisting force provided by a pair of magnet 1 fixed to the chamber bottom and magnet 2 on the shutter. Upon the force from the bilayer exceeding the magnetic force Fm, the released shutter quickly opens the top window and closes the bottom window, which corresponds to stage 2.
At stage 3, the dry air above enters the chamber through the open top window, drying the active layer while the shutter is locked by a pair of magnet 3 on the chamber bottom and magnet 4 on the shutter. When the force produced by the dried bilayer exceeds the magnetic force, identical to Fm above, the shutter is released to close the top window and open the bottom window, corresponding to stage 4. We plot the diagram of RH in chamber versus the hygropiston stroke, which completes a cycle, in Figure 4B.
In stages 1 and 3, we can calculate the time taken for the bilayer to overcome the magnetic force Fm and the curvature change at the end of each stage by combining the vapor diffusion and elasticity models, as given in Supplementary Text S1. In the model, we neglect the time taken for the vapor to fill the initially dry chamber or escape from the humid chamber because it is much shorter than the time for diffusion in the active layer (Supplementary Text S5 and Supplementary Fig. S4). Furthermore, stages 2 and 4 are much shorter than the other stages, and thus, a cycle approximately takes a sum of durations of stages 1 and 3.
The distance the hygrovehicle travels through a hygrochamber cycle, as shown in Figure 4C, is given by L = Rβ, where R is the radius of the wheel and β is the rotation angle of the axle per cycle. The average speed of the vehicle is
Because the vehicle's linear speed can be varied by altering the gear size and wheel radius, we rather plot the hygrochamber frequency as a function of the active layer thickness and Fm in Figure 4D. We compare our theoretical predictions with the experimental data at fixed ha and Fm in Figure 4E and F, respectively, to verify our model. Figure 4E shows that the frequency decreases with the magnetic force because a stronger magnetic resisting force requires a bilayer to absorb vapor longer to build up enough force to overcome Fm. Figure 4F implies that a thicker active layer helps to move the shutter by easily overcoming the magnetic forces. However, too thick an active layer tended to limit the hygropiston stroke by reducing the amount of curvature change.
In our model, active layers thicker than 50 μm failed to fully open and close the shutters although the magnetic forces were overcome.
Discussion
We have designed, fabricated, and demonstrated three machines that can harness the environmental humidity energy to generate mechanical motions. Furthermore, we have analyzed the machine dynamics, to optimize their design and obtain the best performance within the given parameter space. The designs have commonly revolved around the idea that the spatial gradient of humidity should be used to continually drive the machines without human intervention. Using bilayers consisting of hygroscopically active and inactive layers as an engine, our machines exhibit diverse motions, that is, rolling of the entire machine (hygrowheels), reciprocal rotation (hygroseesaws), and linear locomotion by wheels (hygrovehicles). The hygrowheels and hygroseesaws periodically change the distance of the bilayers from the moist substrate to realize repeated cycles of bending and unbending.
The hygrovehicles periodically change the humidity level surrounding the bilayer whose vertical position is kept constant, using a hygrochamber. The hygroseesaws and hygrochambers commonly use the bilayers as a piston that pushes a slider when wet, but pulls when dry. The mechanical analyses of the machine dynamics allow us to optimize the performance of those machines to obtain speeds or rates of reciprocation several times higher than those poorly designed.
Unlike other soft actuators that respond to external stimuli controlled by artificial means, for example, electric27–29 and magnetic fields,30,31 light intensity,32,33 temperature, 34 and pH, 35 our hygromachines work spontaneously under the spatial humidity gradient. From the thermodynamic point of view, our machines produce work as the chemical potential energy flows from the humid down to dry regions, just as heat engines produce work as the thermal energy is transported from a hot to cold reservoir. If the machines work in the spatial gradient of humidity naturally established around human skin, wet soil, shores, and kitchen or laboratory benches, they can be regarded as energy harvesters or scavengers, as exemplified by our hygroseesaws lightening up LEDs. Environmentally sustainable and scalable energy generation can be achieved through the utilization of multiple hygroseesaws; batteries can be used to store the generated power.
A variety of machines are possible by connecting shafts, gears, or cargos to our hygromachines for a range of purposes including those attempted in this work. In particular, these machines can be utilized to explore small inaccessible spaces without the risk of energy depletion.
In Figure 5, we compare the speed of the hygromachines with those of the autonomous and spontaneous robots,14,15,36 as well as robots that are powered by controlled external stimuli.14,23,24,37–43 The hygrowheel demonstrates the highest speed (in body length/min) among the autonomous robots, and the hygrovehicle registers a speed that is comparable with the other robots. When the hygroseesaw is used for electric power generation, it produces a voltage peak of 3 V and power peak of 2.5 mW. Compared with the power generation device of Chen et al., 15 the seesaw's power peak is 50 times higher and the voltage peaks are comparable.

Velocity of various soft robots measured in BL per minute versus BL. The red data points represent autonomous robots; robots powered by external stimuli are represented in other colors. BL, body length.
Despite the advantages in exploiting the ubiquitous energy source of humidity in the atmosphere, there are limitations to be considered in designing hygromachines. First, the actuation speed and the mechanical stiffness of hygroresponsive materials can hardly go together because higher porosity, desirable to expedite vapor diffusion, tends to lower the stiffness. Increasing the hygroscopic response speed of recently reported strong hygroresponsive materials16–19 through porosity control44,45 may resolve this problem. Second, the spatial gradient of humidity is established over a distance on the order of centimeter, which may restrict the size of the actuators. We have shown how to make larger machines using centimetric actuators with hygrowheels and hygroseesaws, promoting more ideas for other purposes. Third, the durability of polymeric actuators can be an issue. It was reported that the bilayer actuators made of PEO can endure 104 cycles of RH variation between 0.2 and 0.8 at 25°C. 14
One way to address this duration issue is by using cheap, disposable, and environmentally benign polymeric materials and fabrication processes. Overall, the future studies to widen the application area of hygroscopic actuation would commonly benefit from our theoretical framework to analyze and optimize the machines driven by hygroresponsive engines.
Footnotes
Authors' Contributions
H.-Y.K. conceived and supervised the research. M.C. and B.S. designed and fabricated the hygromachines, analyzed the experimental data, constructed theoretical models, and performed numerical computations. All the authors discussed the results and prepared the article and figures.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This work was supported by the National Research Foundation of Korea (NRF) grants funded by the Korean Government (MSIT) (Grant Nos. NRF2018R1A3B1052541 and NRF2021M3F7A1017476). H.-Y.K. acknowledges administrative support from SOFT Foundry Institute.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
