Scholarly article on topic 'Designing Dual-functionalized Gels for Self-reconfiguration and Autonomous Motion'

Designing Dual-functionalized Gels for Self-reconfiguration and Autonomous Motion Academic research paper on "Nano-technology"

Share paper
Academic journal
Sci. Rep.
OECD Field of science

Academic research paper on topic "Designing Dual-functionalized Gels for Self-reconfiguration and Autonomous Motion"

Designing Dual-functionalized Gels for Self-reconfiguration and Autonomous Motion

Olga Kuksenok & Anna C. Balazs

Chemical Engineering Department, University of Pittsburgh, Pittsburgh, PA 15261, USA.

Human motion is enabled by the concerted expansion and contraction of interconnected muscles that are powered by inherent biochemical reactions. One of the challenges in the field of biomimicry is eliciting this form of motion from purely synthetic materials, which typically do not generate internalized reactions to drive mechanical action. Moreover, for practical applications, this bio-inspired motion must be readily controllable. Herein, we develop a computational model to design a new class of polymer gels where structural reconfigurations and internalized reactions are intimately linked to produce autonomous motion, which can be directed with light. These gels contain both spirobenzopyran (SP) chromophores and the ruthenium catalysts that drive the oscillatory Belousov-Zhabotinsky (BZ) reaction. Importantly, both the SP moieties and the BZ reaction are photosensitive. When these dual-functionalized gels are exposed to non-uniform illumination, the localized contraction of the gel (due to the SP moieties) in the presence of traveling chemical waves (due to the BZ reaction) leads to new forms of spontaneous, self-sustained movement, which cannot be achieved by either of the mono-functionalized networks.




Received 5 December 2014

Accepted 11 March 2015

Published 30 April 2015

Correspondence and requests for materials should be addressed to A.C.B. (balazs@pitt.

Polymers containing both the SP and BZ functionalities have just recently been synthesized1. To the best of our knowledge, there have been no systematic studies of these dual-functionalized "SP-BZ" gels. Here, we uncover the properties of these active materials and provide guidelines for controlling the complex synergy between the structural reconfiguration and motion.

Before detailing the behavior of the SP-BZ gels, we first describe the distinctive properties of the mono-functio-nalized polymer networks. Notably, the SP-functionalized poly(N-isopropylacrylamide) (poly(NIPAAm)) gels2-5 can be molded remotely and reversibly into a variety of shapes via non-uniform illumination6. In acidic aqueous solutions and the absence of light, the spirobenzopyran chromophores remain primarily in the open-ring form (the protonated merocyanine form, or McH) and are hydrophilic2,3,7. Illumination with blue light causes the isomerization of these chromophores into the closed-ring form (the spiro form, or SP), which is hydrophobic2,3. Consequently, when this gel is illuminated through a photomask, the areas exposed to light shrink, while the non-exposed areas remain swollen. The resulting relief patterns can be "erased" simply by removing the light source2,6.

The remarkable feature of the BZ gels8,9 is that they are the only known polymer networks that undergo periodic swelling and deswelling in the absence of external stimuli and an imposed flow. (Other remarkable examples of oscillating gels include pH-responsive systems based on the bromate-sulfite (BS) reaction that undergo oscillations in a continuously-stirred tank reactor10 or gel-based oscillators that occur due to a chemo-mechanical feedback mechanisms in multi-component systems such as hydrogel-enzyme oscillators11-13 or self-regulating, homeostatic materials14.) When Ru(bpy)3 is grafted to a poly(NIPAAm) gel and the material is immersed in a solution of BZ reagents, the ensuing chemical reaction provides the fuel for this BZ gel's rhythmic mechanical oscillations. In particular, the gel swells as the Ru catalyst is oxidized and collapses when Ru is reduced by the oscillating reaction. These mechanical oscillations are accompanied by traveling chemicals that propagate throughout the sample (if the gel size exceeds the wavelength of the chemical wave). Millimeter-sized BZ gels can undergo self-oscillations for a few hours, i.e., until the reagents in the host solution are consumed9,15; the system, however, can be readily "re-fueled" by replenishing these solutes16.

The BZ functionality forms a perfect complement to the SP-functionalization because the BZ gels are also responsive to blue light. Non-uniform illumination with blue light can be harnessed to control the pattern ofwave propagation and direct the motion of the BZ gels17-19. With both functionalities in one sample, the illuminated portion of the SP-BZ gel will shrink, altering the gel's overall shape. This, in turn, can drastically affect the

dynamics of the wave propagation and the sample's motion, and thereby lead to completely new dynamic behavior.

The recently synthesized poly(NIPAAm-co-Ru(bpy)3-co-Sp) chains that encompass both the BZ and SP functionality are referred to as "PNRS"1. The SP-BZ gels can be formed by cross-linking PNRS chains1. Published experiments were, however, solely focused on probing the solution properties of the dissolved PNRS chains (i.e., no cross-linking)1. Importantly, these studies revealed that: 1) the chains exhibit an oscillatory coil-globule transition (due to the BZ reaction) and 2) the lower critical solution temperature (LCST) of the polymer solution can be shifted by light (due to the light-induced conversion of the chromophore to the hydrophobic form). These experiments indicate that the SP-BZ gels will integrate the distinctive properties provided by both the chromophores and the catalyst.

To explore potentially interesting behaviors of these new SP-BZ gels, we developed the computational model described in the Methods. (Unless noted otherwise, the simulation parameters and their correspondence to physical values are given in the SI.) In particular, we modified the gel lattice spring model, or gLSM20-22, which combines a finite element approach for solving the elastodynamic equations that characterize the propagation of chemo-mechanical waves in the gel, and a finite difference approximation for solving the terms that characterize the chemical reactions and diffusion of reactive species in the system. The gLSM was originally developed to model the dynamic behavior of BZ gels8,23-26, allowing us to capture the self-oscillation of the gel induced by the internalized BZ reaction. We build on this model to describe the SP-BZ gels; specifically, we now account for the kinetics of the photo-induced isomerization of the SP chromophores, and the corresponding decrease in the hydration of the polymer network caused by the isomerization of the chromophores to their spiro form. Moreover, we account for the effect of light on both the BZ reaction and SP moieties within the SP-BZ gels. The latter effect is introduced through two parameters: the reaction rate constant kL(I(r)), which describes the light-induced ring-closure of the spirobenzopyran moiety, and the coefficient W(I(r)), which accounts for the additional production of bromide ions in BZ reaction in the presence of light27. Both parameters are assumed to be proportional to the light intensity at a given point, I(r). In accordance with experimental findings2, we also assume that the temperature of the SP-functionalized gels remains constant when they are illuminated with blue light. Here, we investigate the dynamic behavior of the dual-functionalized gels in the presence of non-uniform light, focusing on the illumination patterns described below. In addition, we elucidate the effects of varying the total concentration of the spirobenzopyran chromophores, CSt P , on the directed motion of these systems.

The output from the simulations in Fig. 1 illustrates the dynamic behavior in these dually responsive gels. Here, an initially flat sample roughly 6.5 mm X 1.7 mm X 0.2 mm in size is illuminated at both edges, leaving a central non-illuminated region of radius Rd = 20 units (approximately 0.8 mm). The sample morphs into a bent structure that promotes its self-sustained downward motion. Neither the SP nor the BZ gel alone would yield such net translational motion (see Supplementary Information).

The SP-BZ gel encompasses two important properties arising from the BZ functionalization that contribute to the above behavior. First, due to the inter-diffusion of the polymer and solvent in the system, BZ gels autonomously move in a direction that is opposite to the direction of the traveling chemical wave within the material17-19. Second, light of a sufficiently high intensity suppresses oscillations within an illuminated region, causing the wave propagation to originate in the non-illuminated area of the BZ gel17-19.

In Fig. 1, the intensity of the light illuminating the ends is higher than the critical intensity needed to suppress the oscillations in a uniformly illuminated sample of the BZ gel. (In eq. (11) below, we set W = 4 X 1024 > Wc within the illuminated regions.) Hence, the

chemical waves originate only in the dark central region and then propagate to the illuminated ends. Due to the spirobenzopyran chro-mophores, the illuminated ends of the gel shrink. Compared to these collapsed ends, the central region is relatively swollen. This uneven distribution of solvent within the gel causes the central region to bulge out of the plane (in the negative z-direction for the case in Fig. 1). With the bending of the gel, the traveling chemical waves move not only in the lateral direction (the x-direction), but also upward from the depressed center to the ends of the sample (along the positive z-direction). (The pattern development during a single oscillation cycle is shown in Figure S1 a-c.) Due to the inter-diffusion of the polymer and solvent20,21,28, the movement of the chemical wave in the positive z-direction causes the gel to move in the opposite direction17-19; i.e., the negative z-direction, and hence, migrate to the bottom of the simulation box for the example in Fig. 1. (The differences in the chemical patterns among Figs. 1a-c are due to the fact that these images are taken at different phases ofthe chemical oscillations.)

Since we neglect the attenuation of light in our relatively thin samples, the "upward" and "downward" directions of motion occur with equal probability in our simulations (Supplementary Fig. S1, d). If the sample is illuminated from the top, even a small degree of attenuation would break this symmetry and favor the downward bending (due to the higher degree of shrinking on the illuminated top of the sample) and would result in the downward motion shown in Fig. 1.

To characterize the observed motion, we monitor the temporal evolution of the z-coordinate of the central node on the bottom face of the gel, zc for different values of C'SP (Fig. 2). In all cases, the system's dynamics can be divided into two stages. During the first stage, the sample remains flat and localized in one plane (as indicated by the flat portion of the curve). During the second stage, the sample moves downward (along the negative z-direction) with approximately constant velocity. (The letters within the curve for CStP = 0.2 in Fig. 2a correspond to the respective images in Fig. 1.) The seemingly monotonic movement of the sample involves small-scale oscillatory motion, as can be seen from the enlargement of a portion of the curve (Fig. 2a).

Figure 2b reveals two key effects of varying CtSP, indicating means of controlling the autonomous motion. First, the duration of the first stage is the shortest for the highest concentration of CtSP (see inset). During this first stage, the out-of-plane bending occurs due to the mismatch in the degree of swelling of the central region and illuminated, shrunken ends. The ends shrink significantly more at higher C'SP (Supplementary Fig. S2 showing the equilibrium degree of swelling of as a function of C'SP); this results in a higher mismatch between the dark and illuminated regions, and correspondingly, a faster loss in the stability of the in-plane configuration. This behavior can be seen by comparing the early-time images for samples with low and high values of ClP (Fig. 3a and c, respectively).

Second, the velocity of the downward motion is fastest for the lowest CSt P and decreases with an increase in CSt P. Notably, the degree of bending also depends on ClP and is lower for the higher ClP due to the lower degree of swelling (Supplementary Fig. S2) and higher rigidity of the sample; namely, the presence of additional bulky chro-mophores makes the sample less flexible at higher C'SP. Hence, the vertical components of the chemical wave propagating from the center to the ends of the gel are significantly more pronounced for the lower CStP sample, which exhibits a higher degree ofbending (see late-time images in Fig. 3b and d). The latter behavior leads to the faster downward motion at lower CSt P.

We emphasize that both types of functionalization are necessary for a flat, thin gel to undergo structural reconfiguration (bending) that leads to autonomous, directed motion. A pure SP gel would simply bend, but not move under a stationary light6. If the pure BZ

Figure 1 | Dynamics of a SP-BZ gel with total concentration of chromophores C'SP = 0.2. Simulation times are: (a) t = 29088, (b) t = 6 X 104, (c) t = 1.4

X 105 and (d) t = 1.8 X 105. Both ends of the sample are illuminated; radius of the masked (non-illuminated) region in the center is Rd = 20 units (approximately 0.8 mm). The detailed correspondence between the dimensionless simulation values and physical values are given in the SI; the dimensionless units of time and length in our simulations correspond to ~1 sec and ~40 mm, respectively. Color bar indicates the concentration of oxidized ruthenium catalyst, v.

gels were placed in the non-uniform illumination considered here, the sample would undergo spontaneous movement that depends on initial random fluctuations in the concentration of BZ reagents (Supplementary Figs. S3 and S4). In the dual-functionalized gel, the SP chromophores give rise to the controllable out-of plane bending, which modulates the propagation of the traveling BZ chemical waves. These synergistic interactions enable the self-sustained motion shown in Fig. 1.

Our simulations reveal that the light-directed shape-changes and resultant motion are independent of initial fluctuations (Supplementary Fig. S1). (We ran four independent simulations for all the values of in Fig. 2b and in all cases observed the robust bending and directed motion of the sample displayed in Figs. 1-3.). Additional simulations show that such behavior remains robust for a range of radii of the dark region, 8 < Rd # 25. If, however, we increase Rd beyond 25 units, and thereby decrease the size of the ends that undergo light-induced shrinking, we find that the SP-BZ gel undergoes complex motion (Fig. 4), which depends on initial random fluctuations and is similar to that of pure BZ gels (see

Figs. S3 and S4). In Fig. 4, where Rd = 30, the sample effectively reorients to ''escape'' from the illuminated areas and moves towards the dark region in the center. The latter behavior is similar to the response of pure BZ gels, which exhibit negative photo-taxis17-19,29,30. These results indicate that the illuminated region must be above a critical size to drive the gel to undergo the robust out-of-plane bending (and subsequent directed motion).

To demonstrate the generality of behavior seen in Fig. 1, we considered a relatively large, square sample that is 90 X 90 X 5 nodes in size and illuminated on all four edges, leaving the central radius of Rd = 20 in the dark (Supplementary Fig. S5). Again, the out-of-plane bulging of the swollen center is accompanied by the bending of the collapsed edges, resulting in the directed motion of the gel. For the case in Supplementary Fig. S5, all four corners are seen to fold upward and the gel moves downwards. Decreasing the size of the illuminated area so that Rd = 40 eliminates this effect and the sample displays behavior (Supplementary Fig. S6) analogous to that in Fig. 4.

Besides the size of the illuminated region, the actual pattern of illumination also affects the systems' dynamics. When the light and

path ofthe traveling chemical waves in the material. Hence, with the combination of light and the dual-functionalization, the gel's shape, movement and directionality can be altered "on the fly''. In essence, this light-responsive material offers a blank canvass for molding self-propelled, reconfigurable soft robots that could operate in an autonomous manner.


To perform these studies, we modified our gel lattice spring model, or gLSM20 22, which was originally developed to capture the dynamic behavior of BZ gels8,23 26. The total energy of these responsive gels is the sum of the interaction energy, UFH, and the elastic energy associated with the deformation of the gel, Uel. The interaction energy is given as:

Ufh = Vh [(1 -w)ln(1 -W)+Xfh (W,T)W(1 -4 + a(1-W)csp -x' v(1 -4)]

The first two terms in eq. (1) describe the mixing energy of the system. Here, Xfh(4,T) is the polymer-solvent interaction parameter31, which depends on the polymer volume fraction, 4 , and temperature, T. The last two terms in eq. (1) describe changes in the hydration of the gel due to interactions involving the chemically anchored moieties. The parameters a and x* in eq. (1) are the coupling parameters that characterize the strengths of the respective interactions. By setting a > 0, we can account for the photo-induced shrinking of the gel6 caused by the photo-conversion of the chromophores to the spiro (SP) form2,3. Here, cSP is the concentration of chromo-phores in the SP state normalized by the total chromophore concentration, CSP. By setting x* > 0, we can account for relative swelling of the gel due to the hydrating effect of the oxidized ruthenium catalyst on the network21. The concentration of oxidized ruthenium catalyst in the system is given by v. The factor I3 in front ofeq. (1) is given by I3 = det B and is an invariant of the left Cauchy-Green (Finger) strain tensor B; it characterizes the volumetric changes in the deformed gel32. The local volume fractions of polymer in the deformed state, 4 , and undeformed state, 4 0, are related as Ref. 21 4 = 4qI3 1=2.

The elastic energy contribution, Uei, to the total energy describes the rubber elasticity of the cross-linked network33,34. It is proportional to the crosslink density, cq, as follows:

Ul = ^ (I1-3- ln If),

Figure 2 | (a) Time evolution of the z-coordinate of the bottom face of the center of the gel, zc, for the simulation in Figure 1. Points marked (a)-(d) correspond to the respective images in Fig. 1 a-d. (b) Time evolution of zc for gels at four different values of CtSP. Inset shows evolution at early times.

dark regions for both the rectangular (Fig. 1) and square-shaped samples (Supplementary Fig. S5) are reversed (so that the center is illuminated and the ends are held in the dark), the more swollen regions are located at the edges of the gels (Supplementary Fig. S7). Hence, the samples lack a central "focal" point (bent region) that leads to the directed motion; hence, the resulting motion in these cases strongly depends on initial fluctuations (Supplementary Fig. S7).

While certain modes of non-uniform illumination can be harnessed to drive the gel's directed motion, uniform illumination over the entire sample can be exploited to arrest the motion. At this light intensity (W > Wc), the oscillations in the BZ gel would be suppressed. Moreover, under the uniform illumination, the entire sample would exhibit the same equilibrium degree of swelling, and hence, it would straighten out to its initial, flat state. The prior motion can, however, be reinstated by illuminating just the ends of the gel.

With this newly derived model, we could probe the properties of a new material, which permits remarkable forms of controllable, bio-mimetic movement. Namely, the photo-induced, coordinated contraction and expansion ofthe gel is coupled to the chemo-mechanical transduction (enabled by the internalized reaction) to produce self-sustained motion. The process of structural reconfiguration leading to net movement resembles the action in muscle-driven motion. Notably, the BZ functionality allows these materials to move autonomously, without the need for external stimuli. In effect, the system generates the energy to power its motion. The SP-functionality allows one to not only reshape and bend the material, but also capitalize on the photosensitivity of the BZ gel and thus, modulate the

where v0 is the volume of a monomeric unit and I1 = trB 32.

Equations (1) and (2) yield the following constitutive equation for SP-BZ gels20,21:

S=-P(4,V,Csp )I + co vo -J4- B, (3)

where I is the unit tensor, S is the dimensionless stress tensor measured in units of v0{1 kT. The isotropic pressure in eq. (3) depends on both the concentration of oxidized ruthenium catalyst, v, and concentration ofchromophores in the spiro form, cSP, and can be written as:

P(4,v,csp)=-(4 + ln(1-<4) + x(<4,T)<42) +covo4 (2Jo)-1-acspW + x'vJ. (4)

Here, x(4,T) = x0(T) 1 x14, where x0(T) = [Sh — TSs]/kT, with Sh and 5s being the respective changes in the enthalpy and entropy per monomeric unit of the gel31.

We describe the dynamics of these gels within the framework of the two-fluid model34, where both the respective polymer and solvent velocities, v(p) and v(s), contribute to the total velocity of the system as v = 4v(p) 1 (1 —4)v(s). We further assume that only the polymer-solvent inter-diffusion contributes to the gel dynamics20,21,28; therefore, we set v= 0 20,21.

The above equations for the elastodynamics of the gel must be supplemented by kinetic equations for the chromophores, which undergo the following inter-conversion:

McH u SP

where kL and kD are the reaction rate constants for the forward and backward reaction, respectively. Typically, spontaneous conversion back to the McH form is significantly slower than the photo-induced isomerization to the SP form3,7. This inter-conversion reaction can be described as Ref. 6:

^ =kL(I(r))(1-csP)-

kdcsp ■

The reaction rate constant kL is proportional to the local light intensity, I(r). We consider only thin samples, and thus, neglect the attenuation of the light across the thickness of the sample6. The photo-stationary concentration of chromophores in the SP form yields6 csP = (1 + kD/kL(I))-1-

We must also account for the reaction kinetics in BZ gels. This can be captured by a modified version20,21 of the two-variable Oregonator model35,36 that explicitly accounts for the polymer volume fraction, 420,21.

Figure 3 | Dynamics of a SP-BZ gel at low concentration of chromophores, C£p = 0.1 (left column) and high concentration of chromophores, C£p = 0.4 (right column). Simulation times are: (a) t 5 29056 and (b) t 5 1.3 x 105 for C£p = 0.1, and (c) t 5 29044 and (d) t 5 1.3 x 105 for Cp = 0.4.

As noted above, the polymer and solvent move with the respective velocities of v(p) and v(s). Hence, the chemical reactions within the SP-BZ gels occur simultaneously with the following three dynamic processes: 1) the movement of the grafted oxidized catalyst, v, and chromophores, csp, with the polymer at the velocity v(p), 2) the transport of the activator for the BZ reaction, u, along with the solvent at the velocity vs) = — W/(1_ W)v(p), and 3) the diffusion of the dissolved activator u throughout the polymer network with a diffusion flux given by21 j(u) = — (1 — W)V(u(1 — W)—1). Hence, the dynamics of the SP-BZ gels can be described by the following set of equations:

= -V.(Csp v<p') + ki(1(r))(1-Csp)-kDCSp,

= -V.(v v(p))+ G(u,v,w

L« = -V.(u v(s))-V.j(")+F(«,i

G(u,v,W) = (1 —W)2u-(1 -

F(u,v,W) = (1 -W)2«-«2 -(1 -f)[/v + ®(1(r))

u-q(1 -W M + q(1 - W

The dimensionless parameters q, /, and e in the above equations have the same meaning as in the original Oregonator model35. The dimensionless variable Ф(/(г)) in eq. (11) accounts for the additional production of bromide ions in the presence of light27 and is dependent on the light intensity /. In our simulations of BZ gels, the above approach allowed us to reproduce17,18 the experimentally observed suppression of oscillations within BZ gels by visible light above a critical intensity25.

Within the framework of the gLSM model20,21,37, the dynamics of the polymer network is assumed to be purely relaxational38, so that the forces acting on the deformed gel are balanced by the frictional drag due to the motion of the solvent, resulting in the following force balance equation21:

V._= V0T-1D« f(W)(v^)-v(s)).

The terms G(u,v,W) and F(u,v,W) in eqs. (8)-(9) are the reactive terms that characterize the BZ reaction and are written as:

Here, f (W) is the friction coefficient, Dw is the diffusion coefficient of the activator, and S is defined above (eq. (3)). We assume38 that the gel/solvent system satisfies the incompressibility condition, V-v = 0; in addition, we set the total velocity, v ; (1 — W)v(s) = 0 21. In other words, it is solely the polymer-solvent inter-diffusion that contributes to the gel dynamics21,38 and, correspondingly, there is no net momentum exchange between the gel and the external fluid. Hence, by setting v = 0 21,38, we

Figure 4 | Dynamics ofa SP-BZ gel with CSp =0.2. Simulation times are: (a) t = 4 X 104,(b) t = 9 X 104, (c) t = 1.3 X 105and(d) t = 1.5 X 105. Both ends of the sample are illuminated; radius of the non-illuminated region in the center is Rd = 40.

neglect the hydrodynamic interactions within the gels. In fact, neutral, non-responsive polymer gels were utilized as a medium for the BZ reaction to specifically suppress the hydrodynamics effects39. In a previous study, we showed that we can neglect hydrodynamic effects even for the case of multiple, interacting gels due to the slow dynamics of the gels and the low viscous forces in the system40.

The effect of light in our SP-BZ gels is introduced through the reaction rate constant kL(J(r)) in eq. (7) and the coefficient W(J(r)) in eq. (11); both values are assumed to be proportional to the light intensity at a given point, J(r). Notably, the temperature of the SP-functionalized gels remains constant when it is illuminated with blue light2, and hence, it is only the light-induced ring closure that causes the dehydration of the polymer matrix. In other words, the physical origin of the photo-induced volume change of the SP-functionalized gels is distinctly different from gel collapse arising from direct light-induced heating.

1. Yamamoto, T. & Yoshida, R. Self-oscillation of polymer and photo-regulation by introducing photochromic site to induce LCST changes. React. Funct. Polym. 73, 945-950 (2013).

2. Szilagyi, A. et al. Rewritable microrelief formation on photoresponsive hydrogel layers. Chem. Mater. 19, 2730-2732 (2007).

3. Satoh, T., Sumaru, K., Takagi, T. & Kanamori, T. Fast-reversible light-driven hydrogels consisting of spirobenzopyran-functionalized poly(N-isopropylacrylamide). Soft Matter 7, 8030-8034 (2011).

4. Sugiura, S. et al. Photoresponsive polymer gel microvalves controlled by local light irradiation. Sensor Actuat. A-Phys 140, 176-184 (2007).

5. Sumaru, K., Ohi, K., Takagi, T., Kanamori, T. & Shinbo, T. Photoresponsive properties of poly(N-isopropylacrylamide) hydrogel partly modified with spirobenzopyran. Langmuir 22, 4353-4356 (2006).

6. Kuksenok, O. & Balazs, A. C. Modeling the Photoinduced Reconfiguration and Directed Motion of Polymer Gels. Adv. Funct. Mater. 23, 4601-4610 (2013).

7. Satoh, T., Sumaru, K., Takagi, T., Takai, K. & Kanamori, T. Isomerization of spirobenzopyrans bearing electron-donating and electron-withdrawing groups in acidic aqueous solutions. Phys. Chem. Chem. Phys. 13, 7322-7329 (2011).

8. Yoshida, R., Kokufuta, E. & Yamaguchi, T. Beating polymer gels coupled with a nonlinear chemical reaction. Chaos 9, 260-266 (1999).

9. Sasaki, S., Koga, S., Yoshida, R. & Yamaguchi, T. Mechanical oscillation coupled with the Belousov-Zhabotinsky reaction in gel. Langmuir 19, 5595-5600 (2003).

10. Horvath, J., Szalai, I., Boissonade, J. & De Kepper, P. Oscillatory dynamics induced in a responsive gel by a non-oscillatory chemical reaction: experimental evidence. Soft Matter 7, 8462-8472 (2011).

11. Dhanarajan, A. P., Misra, G. P. & Siegel, R. A. Autonomous chemomechanical oscillations in a hydrogel/enzyme system driven by glucose. J. Phys. Chem. A 106, 8835-8838 (2002).

12. Mujumdar, S. K., Bhalla, A. S. & Siegel, R. A. Novel Hydrogels for Rhythmic Pulsatile Drug Delivery. Macromol. Symp. 254, 338-344 (2007).

13. Siegel, R. A. Autonomous rhythmic drug delivery systems based on chemical and biochemomechanical oscillators. in Chemomechanical Instabilities in Responsive Materials Ch. 7, 175-201 (Springer, 2009).

14. He, X. etal. Synthetic homeostatic materials with chemo-mechano-chemical self-regulation. Nature 487, 214-218 (2012).

15. Miyakawa, K., Sakamoto, F., Yoshida, R., Kokufuta, E. & Yamaguchi, T. Chemical waves in self-oscillating gels. Phys Rev E 62, 793-798 (2000).

16. Chen, I. C. et al. Shape- and size-dependent patterns in self-oscillating polymer gels. Soft Matter 7, 3141-3146 (2011).

17. Dayal, P., Kuksenok, O. & Balazs, A. C. Using Light to Guide the Self-Sustained Motion of Active Gels. Langmuir 25, 4298-4301 (2009).

18. Dayal, P., Kuksenok, O. & Balazs, A. C. Designing autonomously motile gels that follow complex paths. Soft Matter 6, 768-773 (2010).

19. Lu, X. et al. Photophobic and phototropic movement of a self-oscillating gel. Chem. Commun. (Camb) 49, 7690-7692 (2013).

20. Yashin, V. V. & Balazs, A. C. Pattern formation and shape changes in self-oscillating polymer gels. Science 314, 798-801 (2006).

21. Yashin, V. V. & Balazs, A. C. Theoretical and computational modeling of self-oscillating polymer gels. J. Chem. Phys. 126, 124707 (2007).

22. Kuksenok, O., Yashin, V. V. & Balazs, A. C. Three-dimensional model for chemoresponsive polymer gels undergoing the Belousov-Zhabotinsky reaction. Phys. Rev. E 78, 041406.1-041406.16 (2008).

23. Yoshida, R., Takahashi, T., Yamaguchi, T. & Ichijo, H. Self-oscillating gel. J. Am. Chem. Soc. 118, 5134-5135 (1996).

24. Yoshida, R. Self-oscillating polymer and gels as novel biomimetic materials. B. Chem. Soc. Jpn. 81, 676-688 (2008).

25. Shinohara, S., Seki, T., Sakai, T., Yoshida, R. & Takeoka, Y. Photoregulated Wormlike Motion of a Gel. Angew. Chem. Int. Edit. 47, 9039-9043 (2008).

26. Kuksenok, O. et al. Exploiting gradients in cross-link density to control the bending and self-propelled motion of active gels. J Mater. Chem. 21, 8360-8371 (2011).

27. Krug, H. J., Pohlmann, L. & Kuhnert, L. Analysis of the modified complete Oregonator accounting for oxygen sensitivity and photosensitivity of Belousov-Zhabotinsky systems. J. Phys. Chem. 94, 4862-4866 (1990).

28. Boissonade, J. Self-oscillations in chemoresponsive gels: A theoretical approach. Chaos 15, 023703 (2005).

29. Kuksenok, O. et al. Chemo-responsive, self-oscillating gels that undergo biomimetic communication. Chem. Soc. Rev. 42, 7257-7277 (2013).

30. Dayal, P., Kuksenok, O. & Balazs, A. C. Directing the Behavior of Active, Self-Oscillating Gels with Light. Macromolecules 47, 3231-3242 (2014).

31. Hirotsu, S. Softening of bulk modulus and negative Poisson's ratio near the volume phase transition of polymer gels. J. Chem. Phys. 94, 3949-3957 (1991).

32. Atkin, R. J. & Fox, N. An Introduction to the Theory of Elasticity. Longman: New York, 1980.

33. Hill, T. L. An Introduction to Statistical Thermodynamics. Addison-Weley: Reading, MA, 1960.

34. Onuki, A. Theory of Phase-Transition in Polymer Gels. Adv. Polym. Sci. 109, 63-121 (1993).

35. Tyson, J. J. & Fife, P. C. Target patterns in a realistic model of the Belousov-Zhabotinskii reaction. J. Chem. Phys. 73, 2224-2237 (1980).

36. Tyson, J. J. A Quantitative Account of Oscillations, Bistability, and Traveling Waves in the Belousov-Zhabotinskii Reaction. in Oscillations and Traveling Waves in Chemical Systems. Ch. 3, 93-144 (Wiley, 1985).

37. Yashin, V. V., Kuksenok, O. & Balazs, A. C. Modeling autonomously oscillating chemo-responsive gels. Prog. Polym. Sci. 35, 155-173 (2010).

38. Barriere, B. & Leibler, L. Kinetics of solvent absorption and permeation through a highly swellable elastomeric network. J. Polym. Sci. Polym. Phys. 41, 166-182 (2003).

39. Yamaguchi, T., Kuhnert, L., Nagy-Ungvarai, Z., Mueller, S. C. & Hess, B. Gel systems for the Belousov-Zhabotinskii reaction. J. Phys. Chem. 95, 5831-5837 (1991).

40. Dayal, P., Kuksenok, O., Bhattacharya, A. & Balazs, A. C. Chemically-mediated communication in self-oscillating, biomimetic cilia. J. Mater. Chem. 22, 241-250 (2012).


The authors gratefully acknowledge financial support from the United States Department of Energy under grant number DE-FG02 90ER45438 (for studies on the dual-functionalized SP-BZ gels) and the United States Army Research Office under grant number W911NF-13-1-0489 (for studies on the pure BZ gels).

Author contributions

O.K. and A.C.B. designed the study; O.K. developed the model and numerical code and performed the simulations. Both authors wrote the manuscript.

Additional information

Supplementary information accompanies this paper at scientificreports

Competing financial interests: The authors declare no competing financial interests. How to cite this article: Kuksenok, O. & Balazs, A.C. Designing Dual-functionalized Gels for Self-reconfiguration and Autonomous Motion. Sci. Rep. 5, 9569; DOI:10.1038/ srep09569 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder in order to reproduce the material. To view a copy of this license, visit