Self-oscillating chemoelectrical interface of solution-gated ion-sensitive field-effect transistor based on Belousov–Zhabotinsky reaction

The Belousov–Zhabotinsky (BZ) self-oscillation reaction is an important chemical model to elucidate nonequilibrium chemistry in an open system. However, there are only a few studies on the electrical behavior of pH oscillation induced by the BZ reaction, although numerous studies have been carried out to investigate the mechanisms by which the BZ reaction interacts with redox reactions, which results in potential changes. Needless to say, the electrical characteristic of a self-oscillating polymer gel driven by the BZ reaction has not been clarified. On the other hand, a solution-gated ion-sensitive field-effect transistor (ISFET) has a superior ability to detect ionic charges and includes capacitive membranes on the gate electrode. In this study, we carried out the electrical monitoring of self-oscillation behaviors at the chemoelectrical interface based on the BZ reaction using ISFET sensors, focusing on the pH oscillation and the electrical dynamics of the self-oscillating polymer brush. The pH oscillation induced by the BZ reaction is not only electrically observed using the ISFET sensor, the electrical signals of which results from the interfacial potential between the solution and the gate insulator, but also visualized using a large-scale and high-density ISFET sensor. Moreover, the N-isopropylacrylamide (NIPAAm)-based self-oscillating polymer brush with Ru(bpy)3 as a catalyst clearly shows a periodic electrical response based on the swelling–deswelling behavior caused by the BZ reaction on the gate insulator of the ISFET sensor. Thus, the elucidation of the electrical self-oscillation behaviors induced by the BZ reaction using the ISFET sensor provides a solution to the problems of nonequilibrium chemistry.

www.nature.com/scientificreports/ ion sensitive electrodes 7,8 . On the other hand, the changes in the concentrations of catalysts and intermediates such as Br − seemed to have been analyzed without considering the change in [H + ] (pH) in many studies [9][10][11][12][13][14] . However, the cyclic behavior of pH based on the BZ reaction was not necessarily observed using an ion sensitive array 7 , nor was it imaged using an ion-sensitive field-effect transistor (ISFET) 8 , which is utilized in this study. Therefore, the monitoring and imaging of the pH self-oscillation based on the BZ reaction is valuable for the elucidation of nonequilibrium reactions. The BZ reaction proceeds in a nonlinear chemical reaction system under open and batch conditions 15,16 . A polymer gel is a cross-linked three-dimensional polymer network that swells upon immersion in a solvent, providing an open system for heat and mass transfer. So far, we have found a self-oscillating system in a stimulusresponsive polymer gel based on N-isopropylacrylamide (NIPAAm) that is driven by the BZ reaction. In particular, the poly[NIPAAm-co-Ru(bpy) 3 ] gel clearly exhibited periodical volume changes (swelling and deswelling) in an aqueous solution containing the reactants of the BZ reaction, except for the catalyst 5,17 . This indicates that the chemomechanical self-oscillation of the polymer gel is realized in a closed solution without external stimuli. Therefore, the self-oscillating polymer gel is expected to be applied to the observation of peristatic motions in biological systems and actuators in bioanalytical systems [18][19][20] . In addition, the chemomechanical energy derived from the self-oscillating polymer gel may be converted into electrical energy for autonomous sensors when its electrical properties are clarified.
Considering the above, we examined the self-oscillating chemoelectrical interface based on the BZ reaction using a solution-gated ISFET sensor. This sensor was proposed for detecting ions and biomolecules in biological environments 21,22 . In this device, an electrolyte solution induces the interfacial potential between the solution and the gate insulator instead of a metal gate in a metal-oxide-semiconductor (MOS) transistor; this device, however, requires a reference electrode in the solution. The gate insulator is often composed of oxide or nitride membranes such as Ta 2 O 5 , Al 2 O 3 , Si 3 N 4 , and SiO 2 ; therefore, hydroxyl groups at the oxide or nitride surface in the solution undergo the equilibrium reaction with hydrogen ions through protonation (-OH + H + ⇄ -OH 2 + ) and deprotonation (-OH ⇄ -O − + H + ). This is why the change in the surface charge is detected from the change in pH without being mostly affected by other cations in accordance with the principle of the field effect 23,24 . Such ISFET sensors mostly have a silicon substrate, but a variety of semiconducting materials have recently been applied to pH-sensitive ISFET sensors 25,26 . Moreover, a large-scale and high-density ISFET (arrayed-gate ISFET) sensor has recently been developed with the progress of complementary MOS (CMOS) integrated circuit technologies [27][28][29][30][31][32][33][34] . Using the arrayed-gate ISFET sensor, one can image the change in pH in the solution by detecting it on each gate insulator, which is called the pH image sensor 27 . In a previous paper, an ideal pH sensitivity of approximately 55 mV/pH, which was near the Nernstian response at the electrolyte solution/gate insulator interface, and a high frame speed of 1933 fps were reported 34 . Notably, this sensor had the smallest pitch of 2 mm, which resulted in a high spatial resolution.
In addition, the gate insulator surface of the ISFET sensor was chemically modified to specifically and selectively detect biomolecules and to reduce electrical noise from interfering species 35,36 . Also, such chemical modifications on the gate insulator surface contributed to the electrical detection of the phase transition (swelling and deswelling) in the poly(NIPAAm) gel 37 . Therefore, the self-oscillating polymer gel-coated gate ISFET sensor is expected to be used for the electrical monitoring of self-oscillation based on the phase transition on the modified surface, which is driven by the BZ reaction.
In this paper, we report the electrical monitoring of self-oscillation at the chemoelectrical interface based on the BZ reaction using the ISFET sensors, focusing on the pH oscillation and the electrical dynamics of the self-oscillating polymer brush.

Results and discussion
Real-time monitoring of pH oscillation induced by BZ reaction with single-gate ISFET sensor. The Ta 2 O 5 film was used as the gate insulator of the ISFET sensor in this study. The Ta 2 O 5 surface has hydroxy groups in an electrolyte solution, which shows the equilibrium reaction with hydrogen ions depending on pH. Therefore, the change in potential at the solution/Ta 2 O 5 gate insulator interface is electrically output with varying pHs, depending on the Nernstian response (59.2 mV/pH at 25 ℃). The Nernst equation for pH is shown as Eq. 1 and strictly reflected by the coefficient [β/(β + 1)] with the parameter β (Eq. 2): where φ 0 is the interfacial potential at an electrolyte solution/gate insulator interface, pH pzc is the point of zero charge, q is the elementary charge, N S is the site density of hydroxyl groups at the oxide membrane, K a is the equilibrium constant between hydrogen ions and hydroxyl groups, k is the Boltzmann constant, T is the absolute temperature, and C EDL is the capacitance of the electric double layer 21 . The Ta 2 O 5 -gate ISFET shows higher pH sensitivity near the ideal sensitivity owing to the higher density of hydroxy groups at the Ta 2 O 5 surface (ca. 10 15 / cm 2 ) 38 , which is sufficiently high for β/(β + 1) to be assumed as 1. On the other hand, the BZ reaction is induced under a strongly acidic condition. This is why the change in pH during the BZ oscillation may have never been expected. However, the pH oscillation was induced upon adding malonic acid to the BZ presolution (malonicacid-free) and was clearly detected using the single-gate ISFET sensor, as shown in Fig. 1. The actual pH was estimated from the calibration curve (ΔV out with ΔpH) of the single-gate ISFET sensor used in this study (Fig. S1). www.nature.com/scientificreports/ The average of the pH sensitivity was 55.2 ± 2.1 mV/pH (n = 10) and then the pH of the BZ presolution was determined to be around 1.0 (Fig. S1). In the BZ solution with the concentration of χ, the peak-to-peak ΔV out V peak out was approximately 43 mV, which corresponded to ΔpH ≈ 0.8, and the period of the pH oscillation was approximately 67 s ( Fig. 1). Moreover, the peak-to-peak ΔpH based on V peak out and the oscillation period depended on the concentration of the BZ solution, as shown in Fig. 2. The initial concentration ratio for each component in the BZ reaction was constant, regardless of the concentration of the BZ solution (1/2 × χ, χ, and 2 × χ). At the lower concentration (1/2 × χ), the pH oscillation was observed upon adding malonic acid, but it was gradually attenuated over time owing to the consumption of the substrates (Fig. S2A). On the other hand, the pH oscillated more rapidly and lasted longer at the higher concentration (2 × χ). This is because there were sufficient amounts of reactants in the BZ solution. Thus, a linear relationship between the peak-to-peak ΔpH and the period of pH oscillation was found, depending on the concentration of the BZ solution. That is, the pH oscillation induced by the BZ reaction was clearly demonstrated even under the strongly acidic condition because the ISFET sensor can selectively detect the change in pH. However, we should investigate further to more convincingly show how the pH oscillation is involved in other reactions such as a cyclic redox reaction involving metallic catalysts.
The redox reaction of the ruthenium ion complex on the Pt electrode is From the Nernst equation, the electrode potential E , which is induced between the Pt electrode and the reference electrode, is expressed by  Change in interfacial potential (ΔV out ) for BZ reaction detected using single-gate ISFET sensor. Malonic acid was added to the BZ presolution without malonic acid (pH 1.0). ΔV out was translated to pH in accordance with the pH sensitivity of the ISFET sensor (i.e., the Nernstian response) used in this study (Fig. S1). www.nature.com/scientificreports/ where E 0 is the standard electrode potential, R is the gas constant, T is the absolute temperature, and F is the Faraday constant. Therefore, the potential based on the BZ reaction is shown by Indeed, the redox potential of the ruthenium ion complex oscillated with a period consistent with that of [H + ] and [Br − ] oscillations, as shown in Fig. 3. [Br -] was simultaneously measured using a bromide-selective electrode (8005-10C, HORIBA). The redox reaction of the ruthenium ion complex was also confirmed by the change in the color of the BZ solution (Fig. S3). That is, the G component in the RGB color analysis was oscillated, similarly to the redox potential oscillation (orange Ru bpy 2+ 3 from/to light-green Ru bpy 3+ 3 ). Considering the Field, Koros, and Noyes (FKN) mechanism 9 , the BZ reaction takes place via three processes (A, B, and C), as shown in Table S2. The inhibitor Bris removed in Process A; the activator species HBrO 2 is produced autocatalytically and the catalyst is oxidized in Process B; and the inhibitor Bris reformed upon reduction of the catalyst in Process C. From Fig. 3 Fig. 3), although each period was consistent. Considering each reaction in the three processes and the sequence (e.g., Processes A → B → C) shown in Table S2, the sequence for the change in [H + ] also seemed to shift one by one (i.e., Processes B → C → A). That is, the pH oscillation could not be simply explained from the reaction equations shown for the FKN mechanism. In the pH oscillation obtained in the FET measurements, malonic acid was added to the BZ presolution (pH 1.0). Then, dissociated malonic acids may have resulted in the change in [H + ] in the BZ solution with lower pH, showing its equilibrium reaction CH 2 (COOH) 2 ⇄ CH 2 (COOH) COO − + H + ⇄ CH 2 COO − 2 + 2H + , pK a 2.83 and 5.69 . In fact, ΔV out ([H + ]) gradually decreased after the addition of malonic acid, which resulted in pHs 2-3 with the generation of BrCH(COOH) 2 before ΔV out ([H + ]) first increased (the pH oscillation was observed) (Figs. 1 and S2). Therefore, the equilibrium reaction of malonic acid may be considered for the pH oscillation induced by the BZ reaction, on the basis of our consideration. However, the driving force for the increase in [H + ] cannot be explained from this point. Basically, ruthenium ion complexes showed trivalent Ru bpy 3+ 3 Ru bpy 2+ 3 .  This assumption also provides the generation of Ru bpy 2+ 3 as well as H + , which may explain the slight deviation of peak potentials in [H + ] (dotted circle in Fig. 3). From the above considerations, malonic acid with higher pH was added to the BZ presolution with lower pH and then the dissociated malonic acid underwent H + on the basis of its equilibrium reaction, resulting in the decrease in [H + ] (Process X). Subsequently, undissociated malonic acid reacted with Br 2 , which generate BrCH(COOH) 2 , in accordance with the usual process. BrCH(COOH) 2 was incorporated in the ruthenium ion complex-bromomalonic acid subsystem, resulting in [H + ] increase (Process Y). After the increase in [H + ] based on this assumption, malonic acid spontaneously maintained its equilibrium reaction to reduce [H + ] again (Process X). This is why Processes X and Y may have periodically occurred in turn; that is, the pH oscillation may have been electrically observed using the ISFET sensor. However, the peak pH and the period of pH oscillation simulated in a previous work gradually increased over time 8 , whereas in this study they were almost constant (Figs. 1 and S2). Therefore, such a pH oscillation induced by the BZ reaction should be further simulated in the future to elucidate the reaction mechanism. Although there have already been numerous studies in which the reaction rates of catalysts and intermediates were analyzed, changes in [H + ] were not considered. pH oscillation imaging with arrayed-gate ISFET sensor. The large-scale and high-density ISFET sensor (arrayed-gate ISFET sensor) chip, which was developed on the basis of CMOS integrated circuit technologies 34 , enabled the two-dimensional imaging of pH oscillation induced by the BZ reaction. Figure S4 shows ΔV out with the change in pH from 4.01 to 9.18, detected using the arrayed-gate ISFET sensor. The data shown in Fig. S4A were obtained from three ISFETs, which were randomly chosen from among the 256 × 256 ISFETs. ΔV out changed immediately after changing the pH buffer, indicating the good responsivity to the change in pH. ΔV out shifted in the negative direction when pH was changed from 1.68 to 10.01, and then ΔV out returned to the original value in a stepwise manner for the repeated pH measurements, similarly to the single-gate ISFET sensor. This result demonstrates that the Ta 2 O 5 membrane as the gate insulator of the arrayed-gate ISFET sensor developed in this study remained electrochemically stable when immersed in solutions of various pHs and that each ISFET operated stably. Considering the data shown in Fig. S4A, a linear plot of ΔV out vs pH was obtained from the 256 × 256 arrayed ISFETs, 0.5% of which were excluded from the evaluation owing to electrical failures, as shown in Fig. S4B. This calibration curve indicated a pH sensitivity of 55.6 mV/pH near the Nernstian response. Thus, the electrical responsivity of the arrayed-gate ISFET sensor to changes in pH was in good agreement with that of the single-gate ISFET sensor shown in Fig. S1. In addition, the arrayed-gate ISFET sensor measures two-dimensionally the change in pH, which is in contrast to the single-gate ISFET that measures the change in pH at a local point, and can be used as a pH image sensor.  www.nature.com/scientificreports/ Using the arrayed-gate ISFET sensor, we monitored the pH oscillation induced by the BZ reaction in real time, as shown in Figs. 4 and S5. From ΔV out for the arrayed-gate ISFET sensor, the pH (−log[H + ]) oscillation was clearly observed upon adding malonic acid to the BZ presolution at the concentration of χ (Fig. 4), similarly to ΔV out monitored by the single-gate ISFET (Fig. 1). In addition, the arrayed-gate ISFET sensor, as the pH image sensor, visually demonstrated the pH oscillation induced by the BZ reaction in real time via the change in the color, which corresponded to ΔV out , as shown in the movie (Fig. S5). The oscillation period was approximately 30-40 s, whereas V peak out for the arrayed-gate ISFET sensor was approximately 10 mV (ΔpH ≈ 0.2) smaller than that for the single-gate ISFET, despite the similar trend of pH sensitivity. This may be due to the fact that the volume of the measurement solution in the arrayed-gate ISFET sensor was 1 mL smaller than that in the singlegate ISFET (10 mL); moreover, the effect of stirring of the BZ solution was different. Optimal agitation would be necessary to observe the periodic pH oscillation using the ISFET sensors. Thus, the chemical self-oscillation was imaged in real time on the basis of the change in pH using the CMOS integrated circuit device. In addition, to observe the patterns similar to a series of expanding concentric rings of the BZ reaction observed in petri dishes (target pattern) 39 , several arrayed-gate ISFET sensor chips may be arranged on a wider area in the future.
Chemoelectrical self-oscillation of polymer brush with single-gate ISFET sensor. As shown in the above sections, the solution-gated ISFET sensor can directly detect the pH oscillation induced by the BZ reaction on the Ta 2 O 5 gate insulator with hydroxy groups, which show the equilibrium reaction with hydrogen ions. On the other hand, such a sensing surface is chemically modified and functionalized to detect target samples by ion-sensitive membranes and biomolecular recognition sites such as DNA probes and enzymes 22,40 . The chemical modification on the Ta 2 O 5 gate insulator contributes to the specific and selective detection of target ions and biomolecules at the expense of losing the original pH sensitivity owing to the blocking of hydroxy groups. In particular, the thermoresponsive poly(N-isopropylacrylamide) (PNIPAAm) brushes grafted on the Ta 2 O 5 gate insulator were evaluated as a bioelectrical interface for biosensing on the basis of the electrical property of the ISFET sensor. In this case, the difference between the capacitances of the PNIPAAm brushes at the swelling state and the deswelling state was found in previous studies. 37,41 Moreover, the PNIPAAm gel copolymerized with ruthenium complexes (one of the catalysts) shows the self-oscillating behavior on the basis of the BZ reaction with the redox reaction of complexes 5,18 . That is, the self-oscillating behavior based on the BZ reaction can also be observed as a change in the capacitance of the PNIPAAm brush with Ru bpy 3+/2+ 3 grafted on the Ta 2 O 5 gate insulator, using the ISFET sensor. The self-oscillating polymer [designed to be a random copolymer containing NIPAAm, Ru(bpy) 3 , and NAPMAm as binding sites of Ru(bpy) 3 ], that is, the poly(NIPAAmr-NAPMAm-r-Ru(bpy) 3 NAPMAm) brush, was successfully grafted on the Ta 2 O 5 gate insulator of the ISFET sensor by SI-ARGET ATRP, followed by the conjugation of Ru(bpy) 3 -NHS ester (Fig. 5). After cleaning the Ta 2 O 5 gate surface, ATRP initiators were tethered at the surface with hydroxy groups by silane coupling. Then, NIPAAm and NAPMAm monomers were copolymerized by grafting from the surface with the initiators for 20 h to control the thickness of poly(NIPAAm-r-NAPMAm) in the aqueous solution with the mixture of reagents for ARGET ATRP. After the ARGET ATRP reaction, Ru(bpy) 3 -NHS ester was covalently reacted with amino groups of the poly(NIPAAm-r-NAPMAm) brush grafted on the Ta 2 O 5 gate insulator. The thickness was measured to www.nature.com/scientificreports/ be approximately 20 nm in atmosphere by AFM (Fig. S6). Using the prepared device, we electrically monitored the self-oscillation of the polymer brush, as shown in Fig. 6. In the measurement, ruthenium ion complexes and malonic acid were not included in the BZ presolution at the concentration of χ, and then malonic acid only was added to the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm) brush-grafted Ta 2 O 5 gate surface. In particular, the glass substrate grafted with the same polymer brush was used to cover over the modified Ta 2 O 5 gate surface to prevent the BZ components such as Brfrom diffusing away from the active surface, that is, to maintain the BZ reaction at the surface (Fig. S7). In fact, the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm) brush-grafted Ta 2 O 5 gate ISFET demonstrated electrical self-oscillation, as detected on the basis of the swelling-deswelling behavior of the polymer brush. On the other hand, the poly(NIPAAm-r-NAPMAm) without Ru(bpy) 3 brush-grafted Ta 2 O 5 gate ISFET (control sensor) showed no electrical responses, although the signals of both sensors gradually and slightly decreased owing to the change in pH upon the addition of malonic acid, as shown in Fig. 6A. Moreover, the difference between the signals of the two sensors ( V diff out ) clearly showed the self-oscillation with an amplitude of approximately 3 mV and a period of 40 s (Fig. 6B). This self-oscillation signal was not caused by the change in pH, judging from the electrical behavior of the control sensor. Therefore, the electrical self-oscillation of the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm) brush-grafted Ta 2 O 5 gate ISFET should have resulted from the change in the capacitance of the polymer brush attributable to its swelling and deswelling behavior, judging from the results reported in previous papers 37,41 . From the view point of ionic charges, however, the modified ISFET sensor may have detected simultaneously the change in the charges of ruthenium ion complexes owing to the redox reaction ( Ru bpy 2+ 3 ⇄ Ru bpy 3+ 3 ) on the Ta 2 O 5 gate surface, similarly to the principle of pH sensing (see "Conclusion" section). In this case, the effect of these charges on ΔV out would have depended on the density of ruthenium ion complexes in the polymer brush grafted on the Ta 2 O 5 gate surface. Also, the increase in the density of positive charges based on the oxidation reaction should have induced the increase in ΔV out , in accordance with the source-follower circuit employed in this study 41 . On the other hand, the effect of the change in the capacitance of the polymer brush on ΔV out should be additionally considered for the self-oscillating polymer (i.e., swelling ⇄ deswelling). In previous works, the thermoresponsive PNIPAAm-brush-grafted Ta 2 O 5 gate ISFET clearly showed ΔV out for the change in the temperature, which induced the swelling-deswelling behavior of the PNIPAAm brush past the lower critical solution temperature (LCST), resulting in the change in the capacitance of the PNIPAAm brush 37,41 . Considering the discussion in previous works, the capacitance of the polymer brush in the swelling state would have been larger than that in the deswelling state owing to the hydration of the polymer brush, resulting in the increase in ΔV out . That is, ΔV out of the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm)-brush-grafted Ta 2 O 5 gate ISFET should have increased on the basis of the phase transition caused by the oxidation of Ru bpy 2+ 3 to Ru bpy 3+ 3 (i.e., swelling). This indicates that the positive shift in ΔV out for the self-oscillating polymer brush may have been synergistically induced by the increases in the density of ionic charges of ruthenium ion complexes (bivalent to trivalent ions) and the capacitance of the polymer brush with swelling, which were primarily induced by the oxidation reaction. Moreover, ruthenium ion complexes were covalently introduced into the polymer brush, which then enhanced the ionization of the polymer brush; therefore, the transition between hydration and dehydration appeared tosmoothly occur at the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm)-brush-grafted Ta 2 O 5 gate surface (Fig. 6). The amount of output signals (i.e., the amplitude) and periods would have been mainly derived from the length and density of the polymer brush chain and the concentration of BZ solution. From the slope of V diff out for its increase or decrease, the reaction rate for the oxidation appeared to be higher than that for reduction 5 . This indicates that the swelling behavior of the poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm) brush would have occurred more smoothly than its deswelling behavior at the Ta 2 O 5 gate surface. Thus, the chemical self-oscillation behavior that resulted in the structural change of the polymer brush with the swelling-deswelling transition was successfully monitored as the electrical self-oscillating behavior using the ISFET www.nature.com/scientificreports/ sensor. This suggests that such chemical oscillations can generate electrical energies for future self-oscillating sensors and actuators.

Conclusion
In this study, we demonstrated the electrical monitoring of self-oscillation behaviors at the chemoelectrical interface based on the BZ reaction between the electrolyte solution and the Ta 2 O 5 gate insulator of the ISFET sensor. First, the pH oscillation induced by the BZ reaction was electrically monitored using the nonmodified single-gate and arrayed-gate ISFET sensors. This is because the ISFET sensors had the superior ability to detect ionic charges; otherwise, the oxide gate surface with hydroxy groups in the solution showed the equilibrium reaction with hydrogen ions with positive charges, depending on the pH at such a chemoelectrical interface. In particular, the arrayed-gate ISFET sensor provided images of the pH oscillation owing to the large-scale and highdensity format. To clarify the mechanism of the BZ reaction with pH oscillation, the ruthenium ion complexbromomalonic acid subsystem may be added to it, considering the simulation of the pH oscillation shown in a previous paper 8 . In addition, the equilibrium reaction of malonic acid should be considered to elucidate the BZ reaction mechanism including the pH oscillation owing to the relatively large difference between the pH of BZ presolution ( ≈ 1) and the pKa of the malonic acid ( ≈ 3) added. Second, the self-oscillating polymer brush, poly(NIPAAm-r-NAPMAm-r-Ru(bpy) 3 NAPMAm), grafted on the Ta 2 O 5 gate insulator by SI-ARGET ATRP, clearly exhibited the periodical electrical response of the ISFET sensor. This self-oscillation behavior was mainly based on the change in the capacitance of the polymer brush derived from its swelling-deswelling behavior based on the BZ reaction on the surface of Ta 2 O 5 gate insulator, that is, at the chemoelectrical interface, although the effect of ionic charges may be considered as well.

Methods
Chemicals. The following chemicals and materials were used in this study. N-isopropylacrylamide Devices. We used two types of ISFET sensor for monitoring the pH oscillation induced by the BZ reaction in this study. One was the single-gate ISFET sensor, which was composed of a silicon-based n-channel depletionmode FET with a Ta 2 O 5 /SiO 2 (100 nm/50 nm) layer as the gate insulator with a width (W) of 340 μm and a length (L) of 10 μm, respectively (ISFETCOM Co., Ltd.). The other type was the arrayed-gate ISFET sensor, which was composed of 256 × 256 pixels with a Ta 2 O 5 (300 nm) layer as the gate insulator. Also, the sensing area of one pixel was 1.42 × 1.42 mm 2 and the pixel pitch was 2 × 2 mm 2 . Other sensor performance parameters are shown in Table S1 and described in a previous paper 34 .
The extended-gate ISFET sensor was used for monitoring a self-oscillation behavior of a polymer brush grafted to the Ta 2 O 5 gate surface. In this case, Ta 2 O 5 film (100 nm) with the surface area of 12 × 12 mm 2 was sputtered on a Au (100 nm)/Cr (20 nm) electrode, which was connected to the gate of a silicon-based n-channel junction-type FET (K246-Y9A, Toshiba), as an extended-gate electrode, considering ease of surface treatments 42,43 .
The Ta 2 O 5 thin film was used as the passivation layer to prevent the leakage of currents as well as the pHresponsive layer in an electrolyte solution. A polycarbonate ring of 10 or 20 mm inner diameter was encapsulated to pour measurement solutions onto the Ta 2 O 5 gate surface (0.5 or 1 mL) as a measurement well. Also, the singlegate ISFET sensor was immersed in a 5 mL measurement solution in a beaker for monitoring the pH oscillation.
Electrical measurement of pH oscillation induced by BZ reaction using ISFET sensors. To confirm the fundamental electrical properties of the single-gate ISFET sensor, its gate voltage (V G )-drain current (I D ) electrical characteristics were measured using a semiconductor parameter analyzer (B1500A, Agilent). A change in V G in V G -I D electrical characteristics was estimated as a threshold voltage (V T ) shift, which was evaluated at a constant I D of 700 μA and a constant drain voltage (V D ) of 2 V. A Ag/AgCl reference electrode with a KCl solution was connected to the measurement solution through a salt bridge. As the pH measurement solution, standard buffer solutions with pHs of 1.68, 4.01, 6.86, 9.18, and 10.01 (Wako Pure Chemical Industries, Ltd.) were prepared. The time course of the change in surface potential at the gate surface (ΔV out ) was monitored using a source follower circuit 35 , with which the potential change at the interface between an aqueous solution and a gate insulator can be read out directly at a constant I D (RadianceWare Inc.). In this study, V D and I D were set to 1 V and 700 μA, respectively. On the other hand, the pH responsivity of the arrayed-gate ISFET sensor was monitored from pH 1.68 to pH 10.01 at the frame speed of 3.4 fps in this study, in accordance with the measurement setup shown in a previous work 34 . The measurement data (ΔV out ) was automatically transformed into output images using a measuring software installed on a PC. For both ISFET sensors, ΔV out was calibrated on the basis of the change in pH to analyze the pH sensitivity.
The BZ solution was composed of 1 M HNO 3 , 1 M NaBrO 3 , 1 M CH 2 (COOH) 2 , and 1 mM Ru(bpy) 3 2+ 2Clin deionized water, the concentration of which is expressed as χ. That is, the BZ solution with 1/2 × χ was adjusted to 0.5 M HNO 3 , 0.5 M NaBrO 3 , 0.5 M CH 2 (COOH) 2 , and 0.5 mM Ru(bpy) 3 2+ 2Clfor each component in deionized water, for example. In the electrical measurement, malonic acid was added to the BZ presolution without www.nature.com/scientificreports/ malonic acid and the BZ solutions were stirred under batch conditions. The total volume of the BZ solution was controlled to be 0.5, 1, or 5 mL, depending on the size of measurement wells. The pH oscillation induced by the BZ reaction was observed on the basis of ΔV out of both ISFET sensors without chemical modifications. Ta  www.nature.com/scientificreports/