DOI:
10.1039/D2RA01422J
(Paper)
RSC Adv., 2022,
12, 15145-15149
Modulation of spatiotemporal dynamics in the bromate–sulfite–ferrocyanide reaction system by visible light
Received
3rd March 2022
, Accepted 13th May 2022
First published on 18th May 2022
Abstract
We have carried out the first systematic study of the effects of visible light on the homogenous dynamics in the bromate–sulfite–ferrocyanide (BSF) reaction. Under flow conditions, the reaction system displayed photoinduction and photoinhibition behavior, and the oscillatory period decreased with the increase of light intensity, which is due to the fact that light irradiation mainly enhanced the negative process and affected the positive feedback. The light effect on positive and negative feedback is studied by analyzing the period length of pH increasing and decreasing in detail. With the increase of light intensity, the period length of pH increasing decreases monotonically, while the period length of pH decreasing changes nonmonotonically. These results suggest that light could be used as a powerful tool to control homogenous dynamics. Results obtained from numerical simulations are in good agreement with experimental data.
1. Introduction
Photo-sensitive chemical reactions play important roles in many physiological processes,1 e.g., photosynthesis,2 circadian clock,3,4 bioluminescence,5,6 and phototaxis.7,8 In addition, light as an external control parameter is particularly appealing in the study of nonlinear chemical reactions due to the easy tunability of intensity and spectral composition.9–11
Many chemical oscillators are known to be sensitive to light irradiation, showing photoinduction and photoinhibition behavior in homogeneous and inhomogeneous systems. In the Ru(bpy)32+-catalyzed Belousov–Zhabotinsky (BZ) oscillator, the oscillatory frequency was influenced by the illumination intensity,12 while the rotating spiral waves transform to a labyrinthine standing-wave pattern by periodic light-perturbation.13 By virtue of the process, photophobic and phototropic movements of a self-oscillating gel hosting the photosensitive BZ reaction were designed to approach its favorable environments spontaneously.14,15 In another case, the photosensitive dissociation of molecular iodine was utilized in the chlorine dioxide-iodine-malonic acid (CDIMA) reaction, Briggs–Rauscher reaction, and Bray–Liebhafsky reaction, leading to the cessation or promotion of oscillations,16–18 and pattern structures modulation11,19,20 under the forcing of light illumination.
In addition, light effects on oscillatory dynamics have recently been studied in other pH chemical oscillations, especially the hydrogen autocatalytic processes.21–26 Photo-induced oscillations and pulse waves have been studied in the hydrogen peroxide–sulfite–ferrocyanide (HPSF) reaction.27 The number of oscillatory-peaks was conveniently controlled by changing the illumination period and duration. Original of the phenomena was accounted for the reversible photolysis of ferricyanide, which was followed by the protonation of newly produced cyanide. In the HPSF reaction, the negative feedbacks of the oscillation contained two loops (the consumption of hydrogen ions by ferrocyanide22 and the autocatalytic formation of hydroxide radical),26 both of which were promoted under the light illumination. Hence, the oscillation frequency was influenced by multiple factors and was not easy to control.
The bromate–sulfite–ferrocyanide (BSF) reaction also shows excellent photosensitivity under visible light. Moreover, the negative feedback of the BSF reaction is mainly the oxidation of ferrocyanide by bromate, which consumes considerable hydrogen ions produced in the positive feedback. Kaminaga and coworkers have proposed that light illumination can be used to directly tune the dynamics of the BSF reaction through accelerating the negative feedback (BrO3−–Fe(CN)64−) and have no effect on the positive feedback (BrO3−–SO32−).23 However, in the combined system (BrO3−–SO32−–Fe(CN)64−), the enhancement of negative feedback must indirectly affect the positive feedback process. Therefore, in this work, we systematically investigated the light effect on the frequency of oscillations in the homogenous medium of BSF reaction system by experimental and numerical methods for revealing the light effect on both negative and positive feedback.
2. Experimental content
2.1 Material
Analytical-grade reagent KBrO3, Na2SO3, H2SO4 and K4[Fe(CN)6]·3H2O (Sinopharm Chemical Regents) were used without purification. Throughout the experiment, the four reaction solutions were prepared immediately before each run by the ultrapure water supplied by a water purification system (Millipore, Milli-Q Jr) and were continuously bubbled by N2 in order to avoid the effect of air oxidation. It has been known to all that Fe(CN)64− are affected by light. For this reason, Fe(CN)64− aqueous solution was stored in the brown bottle wrapped up in shading cloth. In addition, all experiments were done in the dark box to avoid the influence of daylight.
2.2 Method
2.2.1 Experimental section. The flow reactor with a total volume of 13.0 mL was made of quartz glass allowing the light to pass more easily. The following input feed concentrations were kept fixed in all CSTR experiments: [BrO3−]0 = 35 mM, [SO32−]0 = 25 mM, [Fe(CN)64−]0 = 15 mM, [H2SO4]0 = 4.2 mM. Four stock solutions were separately introduced into the reactor by a four-channel peristaltic pump (ISMATIC, Switzerland) through glass capillary tubes.To avoid local acidification, the SO32− and H2SO4 tubes were premixed before their entering into the reactor. The reactor mixture was vigorously stirred by a Teflon-coated magnetic stirrer bar. The thermostat (PolyScience, USA) was used to maintain the reaction temperature of T = 30.0 °C in the reactor during the experiment. The reactor was equipped with a combined pH compound electrode (Cole Parmer, USA) to monitor the pH value. The pH-time data was recorded by a computer through a pH meter and an E-coder (Edaq 201, Australia).
An LED light source (λ = 385 ± 10 nm) was used as an accessory light source for controlling the homogenous dynamics by modulating the light intensity with a digital control unit and calibrating with a photometer (Model 1 L1400A, International Light). The optical path length was 2 cm.
2.2.2 Simulation section. The simulations of homogenous kinetics were carried out with a commercial software package (Berkeley Madonna, error control parameter set at 10−10) for stiff differential equations. The same results were obtained with the error control parameter set at 10−6 and 10−12 and with another stiff algorithm.
3. Results and discussion
3.1 Effect of irradiation intensity on the oscillatory dynamics of BSF reaction systems in a CSTR
In the BSF reaction system, Fe(CN)64− acts as a negative feedback agent, which mainly exists in the form of protonation (HFe(CN)63−) through (R1).22 In the dark, the negative feedback loop consumes six protons combined with (R1) and (R2). Under the light irradiation, the more reactive Fe(CN)5(H2O)3− is formed by the key photosensitive step (R3). Meanwhile, the photo-oxidation process (R5) gives rise to the hydrated electron if the aqueous solution of Fe(CN)64− is exposed in UV light.22 In the present work, the aqueous solution of Fe(CN)64− was is strictly shielded from light to avoid (R5) before entering to the CSTR. Actually, the photo-oxidation process increases with decreasing excitation wavelength, as for excitation wavelength below 315 nm, the primary process is photo-oxidation with a secondary aquation process.28,29 Reinhard et al. also reported that the main product at the excitation wavelength of 266 nm are the products of the photo-oxidation process, while at 355 nm are the product of the photo-aquation process.30 Apart from the influence of wavelength, pH also affects the yield of photo-oxidation, the photo-oxidation process is more important on the base condition than on the acid condition.31 The wavelength of the light source in our experiment is 385 ± 10 nm and the pH range is between 4 and 7, which is an acid environment. Therefore, we propose that the photo-aquation process plays a dominant role in our experiment. |
![image file: d2ra01422j-t1.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t1.gif) | (R1) |
|
![image file: d2ra01422j-t2.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t2.gif) | (R2) |
|
![image file: d2ra01422j-t3.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t3.gif) | (R3) |
|
![image file: d2ra01422j-t4.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t4.gif) | (R4) |
|
![image file: d2ra01422j-t5.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t5.gif) | (R5) |
Fig. 1 shows the photo-sensitive dynamics in the CSTR. Under dark, the system shows a steady low pH state which illustrated that autocatalytic reaction was dominant (Fig. 1a). When the system was illuminated with the light, the low pH state immediately entered into an oscillatory state if irradiation intensity increased to 1.154 mW cm−2 (Fig. 1b). As the increasing of the irradiation intensity, the shape of the oscillatory waveform varied which mainly resulted from enhancing the negative feedback process by (R4). As the oscillatory amplitude ranged from 4 to 7, Fe(CN)64− mainly exists in the protonation form, so the photooxidation process was inhibited, while the main photo-response reaction was the photo-aquation HFe(CN)63−. However, when light is applied at I = 29.3 mW cm−2, the reaction system displayed photoinhibition and reached to low pH steady-state (Fig. 1f).
![image file: d2ra01422j-f1.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-f1.gif) |
| Fig. 1 The kinetic curves under different light intensities in a CSTR. k0 = 2.78 × 10−3 s−1; I = (a) 0 mW cm−2; (b) 1.154 mW cm−2; (c) 3.56 mW cm−2; (d) 11.24 mW cm−2; (e) 26.9 mW cm−2; (f) 29.3 mW cm−2. | |
To establish the relationship between the light intensity and the oscillatory period, we have measured the periods under different light intensities, and the results show in Fig. 2. With increasing light intensity, the whole oscillatory periods change monotonically until the system goes into a steady state. The period length of the pH increasing and decreasing are also analyzed shown in Fig. 3. The pH increasing period represents the time required for the negative feedback process, while pH decreasing represents the time required for the positive feedback process. Fig. 3 shows the influence of light intensity on the positive and negative feedback process of the reaction system. As the increasing of I, the pH rising period decreased firstly when I ≤ 10.26 mW cm−2. At I > 10.26 mW cm−2, there is no obvious effect on the pH rising period. On the contrary, the pH descending length changed nonmonotonically with the increase of light intensity. In other words, when I ≤ 10.26 mW cm−2, visible light irradiation not only affected the negative feedback but also modulated the positive feedback. This is easy to understand, under this range of light intensity, due to the production of more reactive Fe(CN)5H2O3−, the negative feedback loop increased gradually, which resulted in the weaken of positive feedback. However, when I > 10.26 mW cm−2, the pH descending period changed slowly as the increasing the light intensity, this is probably because much more HFe(CN)63− covert to Fe(CN)5H2O3− which produced more HCN and dissociated proton by (R6), resulting in strengthening positive feedback processes. So, the oscillatory period decreased as the increase of the light intensity when I > 10.26 mW cm−2.
![image file: d2ra01422j-f2.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-f2.gif) |
| Fig. 2 The oscillatory period as the function of light intensity in the CSTR at k0 = 2.78 × 10−3 s−1. | |
![image file: d2ra01422j-f3.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-f3.gif) |
| Fig. 3 Measured pH rising period and descending period as a function of light intensity in a CSTR. | |
3.2 Computational results
For explaining the experimental results in the CSTR, the mechanism of BSF reaction is given in Table 1. In the dark, the BSF reaction mechanism can be described by (M1)–(M7) reactions which has been successfully applied to simulate the homogenous dynamics of the BSF reaction.23,32–35 With light illumination, the composition reactions (M8)–(M10) could be considered. In addition, the reaction of Fe(CN)63− (M11) could be considered, which also lead to the decreasing of pH. The rate laws and rate constants used in the simulations are given in Table 2. During simulation, the rate constants of (M1), (M2), (M4), (M5) was obtained from ref. 32. And others were fitted according to our experimental results.
Table 1 Mechanism for the photosensitive BSF reaction
No. |
Reactions |
M1 |
SO32− + H+ ↔ HSO3− |
M2 |
HSO3− + H+ ↔ H2SO3 |
M3 |
3HSO3− + BrO3− → 3SO42− + Br− + 3H+ |
M4 |
3H2SO3 + BrO3− → 3SO42− + Br− + 6H+ |
M5 |
6H2SO3 + BrO3− → 3S2O62− + Br− + 6H+ + 3H2O |
M6 |
H+ + Fe(CN)64− ↔ HFe(CN)63− |
M7 |
6HFe(CN)63− + BrO3− → 6Fe(CN)63− + 3H2O + Br− |
M8 |
CN− + H+ ↔ HCN |
M9 |
![image file: d2ra01422j-t6.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t6.gif) |
M10 |
6Fe(CN)5H2O3− + BrO3− + 6H+ → 6Fe(CN)5H2O2− + 3H2O + Br− |
M11 |
SO32− + 2Fe(CN)63− + H2O → 2Fe(CN)64− + SO42− + 2H+ |
Table 2 Rate laws and constants for the photosensitive BSF reaction
Rate law |
Rate constants |
v1 = k1[SO32−][H+] − k−1[HSO3−] |
k1 = 1010 M−1 s−1, k−1 = 103 s−1 |
v2 = k2[HSO3−][H+] − k−2[H2SO3] |
k2 = 6 × 109 M−1 s−1, k−2 = 108 s−1 |
v3 = k3[BrO3−][HSO3−] |
k3 = 0.13 M−1 s−1 |
v4 = k4[BrO3−][H2SO3] |
k4 = 18 M−1 s−1 |
v5 = k5[BrO3−][H2SO3] |
k5 = 0.7 M−1 s−1 |
v6 = k6[Fe(CN)64−][H+] − k−6[HFe(CN)63−] |
k6 = 1010 M−1 s−1, k−6 = 6 × 1010 s−1 |
v7 = k7[HFe(CN)63−][BrO3−] |
k7 = 0.085 M−1 s−1 |
v8 = k8[CN−][H+] − k−8[HCN] |
k8 = 1010 M−1 s−1, k−8 = 6 × 1010 s−1 |
v9 = αk9[HFe(CN)63−] |
k9 = 0−100 s−1 |
v10 = k10[BrO3−][Fe(CN)5H2O3−][H+] |
k10 = 2 M−2 s−1 |
v11 = k11[Fe(CN)63−][SO32−] |
k11 = 2 M−1 s−1 |
The change of chemical species in time in the CSTR can be described by the following eqn (E1):
|
![image file: d2ra01422j-t7.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-t7.gif) | (E1) |
here
ci,cstr,
ci,0 are the concentrations of
ith chemical species in the CSTR.
si,j is the stoichiometric number of the
ith species in the
jth reaction step,
ri is the reaction rate of the
jth reaction, and
n is the number of reaction steps,
k0 is the flow rate in the reactor.
During the simulation process, we define the photokinetic rate constant K9, which is proportional to the light intensity I (K9 = k9I). The photokinetic factor α = (1 − exp(−2.3A))/A involves the absorption of light by all the components of the system, where A is the total absorbance at the irradiation wavelength of 385 nm, defined as
, where ε is the molar extinction coefficients and l is the optical path length in the CSTR, and
is the molar extinction coefficients of HFe(CN)63−, l is the optimal length of the reactor.
Based on the above mechanism in Table 1, the relationship between the oscillatory periods and photokinetic rate k9 in the CSTR was obtained, as shown in Fig. 4. If k9 equals zero, the system has no oscillations. As increasing k9 to 10 s−1, sustained pH oscillations started. The oscillatory periods gradually decreased with a further increase of k9. When the k9 is greater than 90 s−1, the oscillation is suppressed. The pH rising period and descending period under different k9 were also calculated, as shown in Fig. 5. The variation trend of oscillatory period obtained by simulation was consistent with the experimental results. The monotonic change of pH rising period and non-monotonic change of pH descending period obtained by simulation were also consistent with the experimental results. However, the oscillation period of the simulation results is inconsistent with the experiment, because the rate constants of many reactions cannot be obtained. Our next work will further study photoreaction kinetic of the sub-reaction system for detecting the relevant kinetic parameters.
![image file: d2ra01422j-f4.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-f4.gif) |
| Fig. 4 Calculated oscillation period with the increasing of photokinetic rate k9 in a CSTR. Input concentrations are the same as Fig. 2. | |
![image file: d2ra01422j-f5.tif](/image/article/2022/RA/d2ra01422j/d2ra01422j-f5.gif) |
| Fig. 5 Calculated pH rising period and descending period as a function of photokinetic rate k9 in a CSTR. Input concentrations are the same as Fig. 3. | |
4. Conclusion
To summarize, in this work we observed the light effect on the oscillations in the bromate–sulfite–ferrocyanide reaction system, which provides that light can be used as a convenient control parameter to regulate the reaction kinetics. As the increasing of light intensity, the negative feedback was enhanced, while the positive feedback was modulated. The rate law model involved eleven elemental reactions that can simulate the homogenous phenomenon well with the experimental results. The relationship between the oscillatory dynamics and illumination intensity could be obtained to facilitate further multiple application of the BSF, such as driven for self-assembly of nanoparticles,36 switching the conformation of a DNA molecule,37 and aggregation of supra-amphiphile molecules.38,39
Conflicts of interest
There are no conflicts to declare.
Note added after first publication
This article replaces the version published on 18th May 2022, which contained an error in (R6).
Acknowledgements
This work was supported by “the Fundamental Research Funds for the Central Universities” (grant number: 2021QN1047).
References
- P. L. Gentili and J. C. Micheau, J. Photochem. Photobiol. C Photochem. Rev., 2020, 43, 100321 CrossRef CAS.
- M. J. West-Eberhard, J. A. Smith and K. Winter, Science, 2011, 332, 311–312 CrossRef CAS PubMed.
- Q. Wang, Z. Zuo, X. Wang, L. Gu, T. Yoshizumi, Z. Yang, L. Yang, Q. Liu, W. Liu, Y. Han, J. I. Kim, B. Liu, J. A. Wohlschlegel, M. Matsui, Y. Oka and C. Lin, Science, 2016, 354, 343–347 CrossRef CAS PubMed.
- T. A. LeGates, D. C. Fernandez and S. Hattar, Nat. Rev. Neurosci., 2014, 15, 443–454 CrossRef CAS PubMed.
- S. T. Adams Jr, D. M. Mofford, G. S. Reddy and S. C. Miller, Angew. Chem., Int. Ed., 2016, 55, 4943–4946 CrossRef PubMed.
- M. G. Lagorio, G. B. Cordon and A. Iriel, Photochem. Photobiol. Sci., 2015, 14, 1538–1559 CrossRef CAS PubMed.
- K. Wakabayashi, Y. Misawa, S. Mochiji and R. Kamiya, Proc. Natl. Acad. Sci. U. S. A., 2011, 108, 11280–11284 CrossRef CAS PubMed.
- I. Radu, C. Bamann, M. Nack, G. Nagel, E. Bamberg and J. Heberle, J. Am. Chem. Soc., 2009, 131, 7313–7319 CrossRef CAS PubMed.
- H. Luo, C. Wang, L. Ren, Q. Gao, C. Pan and I. R. Epstein, Angew. Chem., Int. Ed., 2016, 55, 4988–4991 CrossRef CAS PubMed.
- A. Kaminaga, V. K. Vanag and I. R. Epstein, Angew. Chem., Int. Ed., 2006, 45, 3087–3089 CrossRef CAS PubMed.
- R. Nagao, R. C. C. de Miranda, I. R. Epstein and M. Dolnik, J. Phys. Chem. A, 2019, 123, 992–998 CrossRef CAS PubMed.
- L. Ren, B. Fan, Q. Gao, Y. Zhao, H. Luo, Y. Xia, X. Lu and I. R. Epstein, Chaos, 2015, 25, 064607 CrossRef PubMed.
- V. Petrov, Q. Ouyang and H. Swinney, Nature, 1997, 388, 655–657 CrossRef CAS.
- X. Lu, L. Ren, Q. Gao, Y. Zhao, S. Wang, J. Yang and I. R. Epstein, Chem. Commun., 2013, 49, 7690–7692 RSC.
- L. Ren, L. Yuan, Q. Gao, R. Teng, J. Wang and I. R. Epstein, Sci. Adv., 2020, 6, eaaz9125 CrossRef PubMed.
- A. P. Muñuzuri, M. Dolnik, A. M. Zhabotinsky and I. R. Epstein, J. Am. Chem. Soc., 1999, 121, 8065–8069 CrossRef.
- S. Kéki, G. Székely and M. T. Beck, J. Phys. Chem. A, 2003, 107, 73–75 CrossRef.
- B. Z. Chowdhry, A. P. Mendham, J. Tetteh and R. Withnall, ChemPhysChem, 2002, 3, 443–446 CrossRef CAS PubMed.
- D. Feldman, R. Nagao, T. Bansagi Jr, I. R. Epstein and M. Dolnik, Phys. Chem. Chem. Phys., 2012, 14, 6577–6583 RSC.
- M. Dolnik, T. Bansagi Jr, S. Ansari, I. Valent and I. R. Epstein, Phys. Chem. Chem. Phys., 2011, 13, 12578–12583 RSC.
- Y. Mori, G. Rabai and I. Hanazaki, J. Phys. Chem., 1994, 98, 12968–12972 CrossRef CAS.
- V. K. Vanag, Y. Mori and I. Hanazaki, J. Phys. Chem., 1994, 98, 8392–8395 CrossRef CAS.
- A. Kaminaga, G. Rábai, Y. Mori and I. Hanazaki, J. Phys. Chem., 1996, 100, 9389–9394 CrossRef CAS.
- G. Rábai and I. Hanazaki, J. Am. Chem. Soc., 1997, 119, 1458–1459 CrossRef.
- A. Kaminaga, G. Rábai and I. Hanazaki, Chem. Phys. Lett., 1998, 284, 109–114 CrossRef CAS.
- Y. Mori and I. Hanazaki, J. Phys. Chem., 1992, 96, 9083–9087 CrossRef CAS.
- Y. Liu, L. Yuan, C. Pan, J. Gao, W. Zhou and Q. Gao, J. Phys. Chem. A, 2018, 122, 1175–1184 CrossRef CAS PubMed.
- M. Reinhard, G. Auboeck, N. A. Besley, I. P. Clark, G. M. Greetham, M. W. D. Hanson-Heine, R. Horvath, T. S. Murphy, T. J. Penfold, M. Towrie, M. W. George and M. Chergui, J. Am. Chem. Soc., 2017, 139, 7335–7347 CrossRef CAS PubMed.
- M. Shirom and G. Stein, J. Chem. Phys., 1971, 55, 3372–3378 CrossRef CAS.
- M. Reinhard, T. J. Penfold, F. A. Lima, J. Rittmann, M. H. Rittmann-Frank, R. Abela, I. Tavernelli, U. Rothlisberger, C. J. Milne and M. Chergui, Struct. Dyn., 2014, 1, 024901 CrossRef CAS PubMed.
- M. Shirom and G. Stein, J. Chem. Phys., 1971, 55, 3379–3382 CrossRef CAS.
- N. Sato, H. H. Hasegawa, R. Kimura, Y. Mori and N. Okazaki, J. Phys. Chem. A, 2010, 114, 10090 CrossRef CAS PubMed.
- G. Rábai, A. Kaminaga and I. Hanazaki, J. Phys. Chem., 1996, 100, 16441 CrossRef.
- T. G. Szanto and G. Rabai, J. Phys. Chem. A, 2005, 109, 5398–5402 CrossRef CAS PubMed.
- A. A. Pogodaev, A. S. Y. Wong and W. T. S. Huck, J. Am. Chem. Soc., 2017, 139, 15296–15299 CrossRef CAS PubMed.
- I. Varga, I. Szalai, R. Mészaros and T. Gilanyi, J. Phys. Chem. B, 2006, 110, 20297–20301 CrossRef CAS PubMed.
- T. Liedl and F. C. Simmel, Nano Lett., 2005, 5, 1894–1898 CrossRef CAS PubMed.
- G. Wang, B. Tang, Y. Liu, Q. Gao, Z. Wang and X. Zhang, Chem. Sci., 2016, 7, 1151–1155 RSC.
- G. Wang, Y. Liu, Y. Liu, N. Xia, W. Zhou, Q. Gao and S. Liu, Colloids Surf., A, 2017, 529, 808–814 CrossRef CAS.
|
This journal is © The Royal Society of Chemistry 2022 |