Open Access Article
Clark D. Gray
a,
Nigel Gottsa,
Cassio Lima
a,
Kenneth Williamsb,
Sebastien Maussang
b and
Royston Goodacre
*a
aCentre for Metabolomics Research, Department of Biochemistry, Cell and Systems Biology, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, BioSciences Building, Crown St., Liverpool L69 7ZB, UK. E-mail: roy.goodacre@liverpool.ac.uk
bRenishaw plc, Spectroscopy Products Division, New Mills, Kingswood, Wotton-under-Edge GL12 8JR, UK
First published on 27th May 2026
The Briggs–Rauscher reaction is a prototypical oscillating system and is based on chemical mixtures of iodate (IO3−), hydrogen peroxide, malonic acid and Mn2+ ions in an acidic solution. Although the Briggs–Rauscher oscillatory behaviour can be seen visibly, conventional analytical methods lack molecular detail of the exact chemical reactions within this highly dynamic system. Using real-time resonance Raman spectroscopy at 785 nm excitation, we track iodine speciation through distinct triiodide and pentaiodide signatures. Fourier-transform analysis of these time-series Raman data reveals that increasing temperature accelerates oscillation frequency while destabilising oscillatory regularity at elevated temperatures (43 °C); higher iodate concentrations promote regular oscillations at intermediate levels but suppress them through excessive oxidation at high concentrations; and starch actively modulates the reaction by stabilising polyiodide formation and doubling the oscillation period relative to starch-free conditions. This study demonstrates Raman spectroscopy as a precise physicochemical tool for analysing chemical oscillations.
These changes are driven by a series of redox and reduction reactions. Initially, IO3− is oxidised to produce molecular iodine using manganese ions (MnSO4) as the metal catalyst:
![]() | (1) |
I2 is reduced to I− by malonic acid:
| 2I2 + 5CH2(COOH)2 + IO3− + H+ → 5CHI(COOH)2 + 3H2O | (2) |
IO3− is further used in the autocatalytic formation of I2, which forms with starch to generate the striking blue colour.
In this reaction, the excess I2 binds to I− ions in a reversible reaction, forming different polyiodide species. Of these, triiodide (I3−) and pentaiodide (I5−) are well-established in aqueous solution and are the species most readily detectable by Raman spectroscopy under the conditions of this study.2,3 Higher-order polyiodides such as I42−, I62−, I7−, and I82− have been characterised primarily in the solid state or in highly concentrated iodide solutions, and their occurrence in dilute aqueous reaction mixtures under BR conditions is not well established.4
Monitoring these reactions relies on both visual recordings and potentiometer probes, with a strong preference for using a ion-selective electrode5 and less of an emphasis on monitoring the chemical constituents, particularly for fast oscillations where it is currently not possible to quench the reaction so that time consuming chemical analyses can be conducted. In reality, there is an incomplete set of 30 postulated intermediary and sub-elementary processes that induce the oscillation of the BR reaction, some of which have had their reaction coefficients established via computational modelling.6,7
Iodide-selective probes are the most widely used in monitoring the BR reaction.8,9 A key limitation of these probes is their poor specificity: rather than measuring a single species, the recorded potential shifts between different potential-determining species throughout the oscillation cycle, as the probe responds to other anionic contaminants and multiple interconverting iodine species.10 Further drawbacks include the need for regular maintenance and calibration,11 low response sensitivity,12 and susceptibility to fouling—the accumulation of waste residue on the electrode surface can interfere with measurements.13
While UV-visible spectroscopy has also been employed,14,15 in which the direct molecular modes have been probed over-time, broad bands, overlapping absorbance spectra, and long integration times with low signal-to-noise ratio can complicate interpretation of which chemical components are present and changing during fast oscillatory reactions. There is therefore a need to be able to measure reaction dynamics in situ in a non-invasive and non-destructive manner.
The study most directly comparable to the present work is that of Chowdhry et al.,16 who employed kinetic resonance Raman spectroscopy at 632.8 nm excitation to monitor the 160 cm−1 pentaiodide band during BR oscillations at different temperatures, demonstrating that oscillation frequency increases with temperature. The present study extends this work in several important respects: (i) a 785 nm excitation laser is employed, which operates further from the λmax = 620 nm of the starch–iodine complex 17–19 and reduces Franck–Condon enhancement, whereby pre-resonance can simultaneous detect both I3− and I5− as well as additional reagent bands (IO3− at 802 cm−1, H2O2 at 877 cm−1, and MnSO4 at 980 cm−1) in the 700–1000 cm−1 region; (ii) Hue-Saturation-Value (HSV) video analysis is synchronised with Raman acquisition to provide simultaneous visual and molecular-level information; (iii) the effects of iodate concentration, malonic acid concentration, hydrogen peroxide concentration, and starch presence are systematically investigated; and (iv) DCLS decomposition combined with Fast Fourier Transform (FFT) analysis is used to extract quantitative oscillation frequencies from the Raman time series. DCLS (Direct Classical Least Squares) decomposes each measured spectrum as a weighted linear combination of reference spectra (in this case, the initial mixed solution spectrum and the final polyiodide-rich product spectrum), and the coefficient (score) assigned to each reference spectrum reflects its relative contribution to the measured spectrum at that time point. These scores are dimensionless and were extracted at each 0.5 s acquisition interval; oscillation amplitude was assessed as the peak-to-trough difference in DCLS score within each cycle. In the visual representation, Hue encodes the dominant wavelength of colour on a circular scale (0–360°, where 0°/360° corresponds to red, ∼120° to green, and ∼240° to blue), saturation describes the intensity or purity of that colour, and value represents overall luminosity. For monitoring the BR reaction, the Hue channel is particularly informative because it is sensitive to the qualitative colour shifts (colourless → amber → blue-black) that accompany changes in iodine speciation, while being relatively insensitive to changes in illumination intensity.
Spectral analysis employed Renishaw's WiRE™ software package allowing DCLS to be performed, using the initial spectrum post-mixing as an intermediary reference and the final spectrum as the product. DCLS data were further smoothed using a Weighted Moving Average (WMA) with a bin-span of 3 (equivalent to 1.22 s, see SI for details on Python scripts). The Raman spectra of each individual solution and their sequential mixtures are presented in Fig. S1, providing the spectral reference framework for the DCLS decomposition applied throughout this study; mode assignments are summarised in Table S1. Solution A (potassium iodate in sulfuric acid) is dominated by a sharp, intense band at 802 cm−1, assigned to the symmetric stretching mode (ν1) of the IO3− ion,21 which serves as the primary spectral marker for iodate throughout the reaction. Additional bands at 980 and 1050 cm−1 are attributable to the symmetric and asymmetric modes of SO4−2 introduced by sulfuric acid.22,23 Solution B (malonic acid, manganese sulfate monohydrate, and starch) displays a more complex spectral profile, with bands at 422, 773, 919, and 980 cm−1. The bands at 919 cm−1 and 773 cm−1 are consistent C–C chain vibrations and O
C–O deformation of malonic acid,22 respectively, while the 422 cm−1 corresponds to an ion pair HSO4− vibration.24 The band at 980 cm−1 reflects the symmetric SO42− stretching mode of manganese sulfate monohydrate. No strong starch-specific Raman features were observed under 785 nm excitation, consistent with the relatively weak Raman scattering intensity of amylose-rich starch systems at near-infrared wavelengths.25 The combined spectrum of Solutions A and B prior to addition of Solution C confirms that no spontaneous reaction occurs between these components and that their individual spectral features are preserved in the mixture.
Solution C (aqueous hydrogen peroxide) is characterised by a prominent band at 877 cm−1, assigned to the O–O stretching mode of H2O2,26 which provides a direct monitor of peroxide concentration throughout the reaction. Upon addition of Solution C to the A + B mixture, the 877 cm−1 H2O2 band is clearly resolved in the combined intermediary spectrum alongside the 802 cm−1 IO3− band, confirming that all reagents are spectrally distinguishable at the onset of the reaction.
The spectrum recorded at the end of the reaction (after oscillations have ceased and the solution has stabilised to a dark blue colour) shows the emergence of new, intense bands in the low-wavenumber region: at 110 cm−1 (triiodide, I3−) and 160 cm−1 (pentaiodide, I5−).
Once the BR reaction has ended, malonic acid is fully consumed and can no longer reduce I2 into active I− ions. This results in the takeover of the oxidation reaction, resulting in the build-up of I2 and its precipitation of solid I2. Bands at 180 and 189 cm−1 correspond to the in-phase and out-of-phase stretching modes of crystalline iodine, respectively.27 The simultaneous appearance of these solid I2 modes indicates that, once malonic acid is fully consumed and can no longer reduce I2 to I−, molecular iodine precipitates from solution. This is also evident with the increase in temperature by ∼5 °C by the oxidation of I−species to I2, releasing heat in an exothermic reaction. The elapsed spectrum, accumulated after the reaction has reached its maximum polyiodide content, shows these low-wavenumber bands at full intensity, with a marked decrease in the 802 cm−1 IO3− band, confirming near-total consumption of iodate. These reference spectra — particularly the initial A + B + C mixture spectrum and the final polyiodide-rich end spectrum—serve as the two end-member components in the DCLS analysis applied to the oscillatory time series, enabling the relative spectral contribution of polyiodide species to be tracked quantitatively as a function of time (Fig. 4a). The term ‘spectral contribution’ refers to the DCLS score assigned to each reference spectrum component at each time point.
An initial analysis was carried out on the first oscillatory cycle of the BR reaction. Fig. 3a presents the colour evolution extracted from the camera recordings (b) together with the associated spectroscopic data. In aqueous solutions, when I2 interacts with excess I− ions, polyiodides such as pentaiodide (I−5) and triiodide (I−3) are produced. A higher level of enhanced Raman scattering was observed for I−5 than for I−3. Pentaiodides are longer and have weaker bonds compared to triiodides; with a delocalised electronic structure, resulting in greater polarizability changes. These results are consistent with Zhang et al.28
Once the solutions are mixed, the reaction mixture undergoes a rapid colour transition: it changes from colourless to yellow, and then quickly shifts to a deep blue/black within approximately 0.8 s. In contrast, the dominant changes in the spectroscopic signal (DCLS) require around 3.3 s to reach their maximum intensity. A consistent temporal offset is also apparent between the two signals: the Hue channel registers a colour change approximately 2.5 s before the DCLS signal reaches its peak, reflecting the fact that the starch–iodine chromophore is visually detectable at very low polyiodide concentrations, whereas the DCLS signal requires measurable accumulation of I−3 and I−5 to register a significant spectral contribution.
This difference in timing indicates that the appearance of the blue starch–iodine colour occurs before the full chemical accumulation of the relevant iodine species. The early blue colour reflects the initial formation of I− and its interaction with I2 to produce the starch–polyiodide complex, which is highly chromatic even at low concentrations. However, the formation of larger amounts of polyiodides requires both accumulated I2 and an excess of I−, a process that proceeds more slowly. As a result, the DCLS signal continues to rise after the visible colour has already appeared. The progressive increase in cycle period visible in both traces towards the end of the reaction reflects the depletion of reagents, most notably IO−3 and malonic acid, which sustain the autocatalytic feedback loop.
A gradual baseline rise is apparent in the DCLS time series over the course of the reaction. This is attributable to the progressive accumulation of coloured species (I2 and polyiodides) in solution, which increases the broad-band fluorescence background and alters the overall spectral baseline. Increasing turbidity as the reaction approaches termination (due to solid I2 precipitation) also contributes to baseline elevation through enhanced scattering.
To determine the oscillation frequency of the BR reaction, the DCLS Raman time-series was subjected to a Fast Fourier Transform (FFT), and the dominant peak(s) in the resulting frequency spectrum were taken as the characteristic oscillation frequencies.
At lower temperatures (5–13 °C), the FFT reveals a set of single, narrow and dominant frequency peaks, indicating regular, well-defined oscillations. As temperature increases toward 22–29 °C, the dominant frequency increases and additional harmonic content becomes apparent. At 43 °C, the DCLS time series adopts a markedly different profile compared with lower temperatures: rather than displaying well-resolved symmetric oscillations, the trace exhibits a compressed, irregular pattern in which successive peaks are poorly separated, and the reaction rapidly transitions to a terminal state. This behaviour indicates that at 43 °C the accelerated reaction kinetics cause reagent depletion to occur on a timescale comparable to the oscillation period itself, destabilising the autocatalytic feedback loop and disrupting regular oscillatory behaviour. The corresponding FFT at 43 °C consequently shows a broadened, less well-defined frequency peak, consistent with the loss of oscillatory regularity. These observations extend the findings of Chowdhry et al.16 who monitored the 160 cm−1 pentaiodide band at 632.8 nm excitation over a narrower temperature range, by providing additional mechanistic detail across a wider temperature window and revealing the onset of dynamical instability at elevated temperatures.
) from the production of iodous acid (HIO2) which drives multiple oscillatory reactions:4,6
IO3− + 2H+ + I− HOI + HIO2
| (3) |
![]() | (4) |
then oxidises manganese(II) (Mn2+) to form hydroxy manganese (Mn(OH)2+) and additional HIO2 to act as an intermediary in the redox cycles:
![]() | (5) |
The concentration of IO−3 is too low at 47 mM, there is insufficient I− to sustain the autocatalytic feedback loop, triggering no observable reaction (see Fig. S5). Short and fast oscillations are observed at 70 mM: attributed to IO−3 quickly depleting I−, generating rapid fluctuations in I−. With the addition of more KIO3, higher concentrations of IO−3 can stabilise regular oscillations. With ever higher concentrations of IO−3 (200 mM), the fast oscillations quickly slow down and disappear due to excessive oxidation. Higher concentrations of KIO3 still leads to slower oscillations, seen in the decreasing dominant frequencies.
The reaction transitions between yellow and clear phases at an accelerated rate, bypassing the characteristic blue-black complexation observed when starch is present (see Fig. 6a). This is noted by the slight change in relative luminosity, oscillating between 86% and 98%. The luminosity oscillation is not symmetrical, with the peak and trough luminosities increasing by approximately 0.5% each consecutive cycle of the reaction.
The yellow colouring of the solution indicates dissolved I2, which has a minimum in the infrared region of its absorbance spectrum.29 The Raman spectroscopic data reveals minimal chemical and physical changes because of this low absorbance of dissolved I2 in the 785 nm region. At a critical point in the reaction cycle, the 877 cm−1 mode of H2O2 undergoes a marked decrease (Fig. S8). Simultaneously, the 980 and 802 cm−1 modes of manganese sulfate and IO−3 respectively, decrease completely. This depletion coincides with a sharp increase in the 180 and 189 cm−1 Raman modes of solid iodine, indicating the abrupt precipitation of I2. HSV analysis reveals a decline in luminosity and a shift in hue, indicating increased turbidity. As molecular iodine accumulates, a thin solid layer forms at the surface, scattering light and causing a luminosity increase as it reflects ambient light.
DCLS data were obtained by assigning the polyiodide spectra and the I2 spectra for the non-starch and starch mixtures respectively (Fig. S8). Fig. 6 shows (a) a clear set of oscillations may be observed with the starch-present BR reaction, whereby no polyiodide oscillations took place with the starch-absent reaction in the DCLS and (b) FFT was applied to the HSV Hue values to obtain the oscillation frequency.
The starch-absent mixture displays an almost doubling in oscillation frequency compared with its starch-present mixture, parallel to the results of Csepei and Bolla.30 Starch plays a crucial role in modulating the BR reaction, delaying molecular iodine from its reduction into I−, and extends the oscillatory period. In its absence, free I2 is more readily reduced, eliminating a stability and leading to a more rapid cycle of the reaction. The absence of starch also leads to the absence of polyiodide production. Starch's helical amylose structure traps I2, creating a charge–transfer complex that shifts the equilibrium toward polyiodide formation, preventing the immediate reduction of I2 back to I−. Without starch, I2 remains in free solution and is quickly reduced, minimising polyiodide formation. With starch, I2 is retained longer, allowing more interaction with I− ions to generate polyiodides.
Supplementary information (SI): additional experimental methods, Raman spectral assignments, pre-processing details, and supporting figures related to the real-time resonance Raman monitoring of the Briggs–Rauscher reaction. This comprises reference Raman spectra of individual reaction components and intermediates, time-resolved Raman and video correlation analyses, weighted moving average pre-processing procedures, and investigations into the effects of malonic acid, potassium iodate, hydrogen peroxide, temperature, and starch on oscillatory spectral behaviour using DCLS and FFT analyses. See DOI: https://doi.org/10.1039/d6ay00742b.
| This journal is © The Royal Society of Chemistry 2026 |