Analytical Methods
View Article Online PAPER View Journal | View Issue
Modeling the response of a control-released ion- selective electrode and employing it for the study Cite this: Anal. Methods,2018,10, 467 of permanganate oxidation kinetics†
Dean Song, ab Rongning Liang, *b Xiaohua Jiang,b Huiqing Sun,a Fanyu Kong,*a Bo Lv,a Qiannan Fanga and Wei Qinb
Although polymeric membrane ion-selective electrodes (ISEs) based on outward ion fluxes have been found analytically useful, there is still a lack of a theoretical framework for this detection system. In this study, we attempted to model the response of this kind of permanganate ISE and employed this ISE to analyze the rapid MnO4 /H2O2 reaction. This response is attributed to H2O2 oxidation with MnO4 that is released from the inner solution to the membrane surface layer. The results show that the experimental data can be fitted well to the proposed model that is elucidated mathematically from the viewpoint of chemical kinetics. The second-order rate constant is determined at a near neutral pH and is in agreement with the acid dissociation law to provide the specific value of 370 M 1 s 1. The kinetic mechanism was then investigated by performing DFT calculations. Via analysis of the Mn–O bond length and the HOMO orbital, it has been found that the studied redox system functions similarly as the so- called hydrogen abstraction mechanism with an energy barrier of 24.5 kcal mol 1. This study is considered to be the first report on the simulation of MnO4 attack at the O–H bond. On the basis of Received 26th November 2017 the transition state theory and previous studies on MnO attack at the C]C and C–H bonds, the Accepted 23rd December 2017 4 relationship between the experimental rate constant and computational energy barrier is finally DOI: 10.1039/c7ay02735d constructed. The result indicates the validity of our proposed method and makes the control-released rsc.li/methods ISE a very promising platform to study the kinetics.
Introduction understanding the mechanism of potentiometric response of the ISE.9 Polymeric membrane ion-selective electrodes (ISEs) are one of To date, it has been fully realized that the transmembrane Published on 26 December 2017. Downloaded by Library of Chinese Academy Sciences 9/11/2018 3:33:27 PM. the most important analytical tools widely used for the analysis concentration gradient can cause ion ux across the polymeric of environmental and biological samples.1 They have received interface and thus deteriorate the detection limit of the ISE by signicant attention owing to the attractive features including perturbing the ion concentration at the phase boundary.10,11 high selectivity, ease of handling, and relatively low-cost.2,3 Although the ion ux cannot be entirely eliminated, this char- Recently, these well-established potentiometric sensors have acteristic of ISE has been recently found to be useful for appli- undergone a quiet revolution to reach the detection limits in the cations in polyion sensors,12 pulstrodes,13 switchtrodes,14 and sub-nanomolar levels.4,5 Inspired by this, many innovative ion-channel biosensors.15 sensing schemes have been proposed for a broad range of Recently, a novel ISE system has been exploited using targets over the last few years.6–8 Accompanied by the techno- outward ion ux (efflux) in the direction of the sample solution – logical advances, great emphasis has also been placed on under zero-current conditions.16 18 In this case, the ISE membrane can play the critical roles of both a polymer matrix for reagent release and a sensitive transducer for potentiometric a The State Agriculture Ministry Laboratory of Quality & Safety Risk Assessment for detection. Particularly, our group has developed an ionophore- Tobacco, Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao 266101, China. E-mail: [email protected] free permanganate ISE that utilizes outward ion uxes through bKey Laboratory of Coastal Environmental Processes and Ecological Remediation, the membrane to provide control-released substrates for 16,19 Yantai Institute of Coastal Zone Research (YIC), Chinese Academy of Sciences (CAS), monitoring the analytes of interest. Previous studies have Shandong Provincial Key Laboratory of Coastal Zone Environmental Processes, shown that the design based on permanganate release allows YICCAS, Yantai, Shandong 264003, P. R. China. E-mail: [email protected] a sensitive and reproducible detection of reductants including † Electronic supplementary information (ESI) available: Experimental ascorbate and dopamine. However, the application of this second-order rate constants and computational energy barriers for the selected ffi permanganate oxidation reactions in water and Cartesian coordinates of the technique is still limited since all reductants that can e ciently optimized structures. See DOI: 10.1039/c7ay02735d
This journal is © The Royal Society of Chemistry 2018 Anal. Methods,2018,10, 467–473 | 467 View Article Online Analytical Methods Paper
react with permanganate ions may produce measured potential shown in Fig. S1 in the ESI.† For each ISE, a disk of 7 mm interference. diameter was punched from the obtained membrane and glued This is obviously disadvantageous for the analysis of to a PVC tube with a THF/PVC slurry. A 0.08 M potassium complex samples, but advantageous for the study of perman- permanganate solution containing 0.02 M NaCl was added to ganate oxidation kinetics under controllable conditions. each electrode as the inner lling solution. Before the Because the potential responses are mainly attributed to the measurements, these electrodes were activated by conditioning consumption of permanganate ion on the membrane surface them in 0.1 M NaCl for 3 days. This pretreatment was used to layer, this permits the selective and continuous observation of make permanganate anion (MnO4 ) enter the ionophore-free the rapid permanganate oxidation process. In comparison, the ISE membrane from the inner lling solution and combine spectrophotometric method together with the stopped-ow with tridodecylmethylammonium cation (TDMAC+) via ion- technique would be restricted by the disturbance from the in exchange between MnO4 and Cl .Aer this step, the 20 + situ formed colloidal manganese oxide. Consequently, it is TDMAC sites were considered to be fully saturated by MnO4 , necessary to model the response behavior of the control- which not only could maintain the relatively stable MnO4 released permanganate ISE with respect to the chemical prin- concentration in the membrane phase but could also facilitate ciple of the reaction. the continuous ion ux from the inner lling solution to the Herein, we attempted to investigate the theoretical founda- sample solution. tion of the efflux permanganate ISE under zero-current condi- tions for the rst time and then employed this ISE as a research EMF measurements platform to study the permanganate oxidation kinetics. As Potentiometric measurements were carried out at 25 C using a proof-of-concept experiment, hydrogen peroxide (H2O2) was an LK2006 electrochemical workstation (Tianjin, China) with chosen to be the objective analyte. Note that the permanganate a saturated calomel electrode as a reference electrode in the titration has long been recognized as a reliable procedure for galvanic cell: SCE//sample solution/ISE membrane/inner lling 21 H2O2 determination. Additionally, the corresponding reaction solution/AgCl/Ag. The 20 mL sample solutions were prepared kinetics has been investigated in an acidic aqueous solution by with 0.1 M NaCl as the background electrolyte and adjusted 22 stopped- ow spectrophotometry. By the construction of with NaOH or HCl to the required pH. For each measurement, a theoretical model, the ionophore-free permanganate ISE may the ionophore-free permanganate ISE was rst taken out from have the opportunity to extend its utilization to describe the the condition solution, rinsed with water, and immersed in the related reaction routes at a near neutral pH. To gain a deeper sample solution under vigorous magnetic stirring to obtain insight into the reaction mechanism, the density functional a potential baseline. The potential change was then initiated by theory (DFT) calculation is adopted thereaer to explain the m adding a certain amount (200 L) of H2O2 and was determined 23,24 reaction mechanism. The obtained results were further used continuously with time. to conrm and explain the outcome of the experimental study. Computational details Experimental The density functional theory (DFT) calculation was performed 25 Chemicals and reagents using the Gaussian 16 package program. The B3LYP model together with the LanL2DZ basis set for the manganese atom Published on 26 December 2017. Downloaded by Library of Chinese Academy Sciences 9/11/2018 3:33:27 PM. Tridodecylmethylammonium chloride (TDMAC), 2-nitrophenyl and 6-31+G(d) basis set for other atoms were employed octyl ether (o-NPOE), and high-molecular weight poly(vinyl throughout this study. The solvent effect of water was taken into chloride) (PVC) were obtained from Fluka AG (Buchs, Switzer- account using the IEFPCM model. It is worth mentioning that land). Other chemicals were of analytical grade and purchased the abovementioned method has been utilized successfully in from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). the mechanistic study of permanganate oxidation in an Aqueous solutions were prepared by dissolving the appropriate aqueous solution.20 The geometry was initially optimized, and U salts in freshly deionized water (18.2 M cm speci c resistance) the vibrational frequency was then calculated to determine the obtained using a Pall Cascada laboratory water system. The pH nature of the stationary point and the critical energy parame- measurements were conducted using a Leici pH meter coupled ters. The transition state was further conrmed by intrinsic with a combined electrode (Shanghai Apparatus Co. Ltd.) that reaction coordinate (IRC) analysis. The Cartesian coordinates was calibrated at pH 4.00, 6.86, and 9.18 before use. could be found in the ESI.† Based on the optimized structures, the corresponding wave functions were also analyzed.23,24 The Preparation of membranes and electrodes three-dimensional diagrams of the optimized structure, the The ionophore-free permanganate ISE was prepared according highest occupied molecular orbital (HOMO), the bond critical to the previous procedures.16 Briey, 240 mg of the membrane point (BCP) in the hydrogen bond path, and the spin density component containing 10 wt% TDMAC, 65 wt% o-NPOE, and were constructed with Avogadro. 25 wt% PVC was dissolved in 2.5 mL of tetrahydrofuran (THF). The resulting solution was then transferred into a glass ring of Model development 36 mm diameter xed on a glass plate to allow the solvent to As shown in Fig. 1, MnO4 can diffuse through the polymeric evaporate completely to obtain a homogeneous membrane, as membrane when the ISE is immersed in the sample solution.
468 | Anal. Methods,2018,10, 467–473 This journal is © The Royal Society of Chemistry 2018 View Article Online Paper Analytical Methods
respectively; k represents the second-order rate constant. The integration of eqn (3) provides
aMnO ;t ln 4 ¼ k C Dt gH O H2O2 (4) 2 2 aMnO4 ;0
Considering the Debye–Huckel¨ rule with the Guntelberg¨ approximation, the activity coefficient of the neutral molecule g ¼ can be taken as unity ( H2O2 1), thus indicating the equivalent a 26 of CH2O2 to H2O2. If H2O2 is applied in a large excess (e.g. at least 5 folds) of MnO4 , CH2O2 can be considered as approxi- mately constant. Eqn (4) is thus simplied and transformed Fig. 1 Schematic of the ionophore-free permanganate ISE under the into a pseudo- rst-order kinetic model as follows: fl zero-current condition and the MnO4 ux from the membrane phase aMnO4 ;t ln ¼ kpDt (5) to the membrane surface layer (i.e., Jin) and further into the sample aMnO4 ;0 bulk (i.e., Jout). where kp is the pseudo-rst-order rate constant equaling to
kCH2O2. Upon substituting eqn (5) into (2), we obtained This is attributed to the transmembrane concentration gradient dE DE RT ¼ ¼ k under the zero-current condition.9 It has been well-established dt Dt F p (6) that the ionophore-free permanganate ISE can exhibit a high discrimination ability towards other anions such as thiocya- From eqn (6), the potential change rate (dE/dt) at the initial nate, nitrate, and chloride.16 Therefore, the observed potentio- stage is found to be a function of kp, which is further propor- metric response (E) serves as an indication of the MnO4 activity tional to CH O . a 2 2 ( MnO4 ) on the membrane surface layer:
0 RT E ¼ E ln aMnO (1) F 4 Results and discussion Experimental modeling of the potentiometric response In the absence of a reductant, the change in E is mainly Fig. 2 displays a representative prole of the potentiometric related to the MnO ux from the membrane phase to the 4 response with time obtained in this study. As discussed in the membrane surface layer (i.e., Jin). The E value will approach to model development, the potentiometric response rst reached be stable gradually when J equals to the MnO ux diffusing in 4 a stable value, serving as an indication of a steady-state process further into the sample bulk (i.e., Jout). This is actually the (i.e., Jin ¼ Jout), which meant a constant aMnO at the membrane a 4 steady-state process with a constant MnO4 released at the surface layer. By adding a small aliquot of H2O2 to the sample membrane surface layer. solution, MnO4 on the membrane surface layer was A er the addition of a reductant (i.e.,H2O2) to the sample
Published on 26 December 2017. Downloaded by Library of Chinese Academy Sciences 9/11/2018 3:33:27 PM. consumed; thus, the aMnO value was decreased, and the ffi a 4 solution, an e cient redox reaction occurs, and the MnO4 at the membrane surface layer is decreased; thus, the measured
potential from E0 at time ¼ 0toEt at time ¼ t is increased:
RT aMnO4 ;t DE ¼ Et E0 ¼ ln (2) F aMnO4 ;0
At the initial stage, the reaction rate is signicant, whereas the MnO4 ux remains at a low level as the steady-state process ¼ D D (i.e., Jin Jout). The E variation ( E) with time ( t) should be merely attributed to the consumption of MnO4 by H2O2. This is analytically useful and can be modelled mathematically from the kinetic point of view hereinaer. By speculating that the MnO4 oxidation of H2O2 is a simple a bimolecular reaction, the decay rate of MnO4 can be expressed by the following second-order equation:
d aMnO4 ¼ kaMnO aH O ¼ kgMnO CMnO gH O CH O (3) dt 4 2 2 4 4 2 2 2 2 fi ffi Fig. 2 A representative pro le of the potentiometric response with where a, g, and C represent the activities, the activity coe - C ¼ m fi time by adding H2O2 33 M. Inset: the corresponding linear tting to cients, and the concentrations of the corresponding species, the model proposed by eqn (6).
This journal is © The Royal Society of Chemistry 2018 Anal. Methods,2018,10, 467–473 | 469 View Article Online Analytical Methods Paper
measured potential was increased. It is worthy to note that the permanganate oxidation kinetics. This is mainly attributed to a MnO4 value at the membrane surface layer is in the range from the inherent advantages, including the high selectivity for m 0.1 to 1.0 M, whereas in the sample bulk, it is below 1.4 nM in MnO4 and continuous monitoring without a quenching step, 1 min (i.e. 80 nM in 1 h).16 Therefore, the prerequisite for of the ISE. The similar ISE system has also been successfully maintaining the pseudo-rst-order condition should be satis- employed to detect the enzyme activity.17,18 Hence, this study
ed using at least 5 mMH2O2. This is conrmed by the result may provide a possibility of studying other reactions by shown in the inset of Fig. 2, where the potential change rate (dE/ controlling the release of suitable reagents such as antigens or dt) at the initial stage can be tted well to eqn (6); this illustrates toxins.28,29 the reasonability of the proposed model. The linear relationship The inset of Fig. 3 shows the obtained k values in corre- a ff also con rmed a rst-order dependence on MnO4 and gave the spondence to di erent pH conditions, which are usually – kp value from the tting slope. Note that the initial stage is ascribed to the distribution of the acid base species of the present for a short time. The deviation between the actual and reactants.27,30 On the basis of the acid dissociation constants ideal potential will become larger and can be attributed to the (pKa ¼ 2.25) of HMnO4/MnO4 and (pKa ¼ 11.75) H2O2/ 31,32 increased Jin caused by the greater transmembrane concentra- HO2 , there should be only one reaction pathway involving tion gradient. The increased Jin is apt to inhibit the consump- the reacting species of MnO4 and H2O2 in the near neutral pH tion of MnO4 by H2O2 to reach a new steady-state process. range. This is supported by the relatively stable k values of 381 Consequently, the abovementioned analysis was only per- M 1 s 1, 378 M 1 s 1, and 351 M 1 s 1 at the pH values of 6, 7, formed at the initial stage when the steady-state process was and 8, respectively. The corresponding specic second-order still valid. rate constant was then averaged to be 370 M 1 s 1 for the ff The variation of kp under di erent reaction condition was MnO4 /H2O2 reaction. This level of reaction rate constant is then examined. In Fig. 3, the determined kp was found to higher than that reported for the oxidation of most organic 33 increase with the increasing CH2O2 added to the sample solution. compounds with MnO4 . Therefore, if the MnO4 /H2O2 The positive linear correlation with a zero intercept indicated reaction is involved in the treatment of environmental a rst-order dependence on CH2O2. The reaction between MnO4 contaminants, it would occur rst and may play a critical role in and H2O2 was thus concluded to be of rst-order with respect to the enhancement of the reaction efficiency by forming a highly each reactant and second-order in total, which was in accor- reactive intermediate. However, only the kinetic result was not dance with the results obtained for other permanganate sufficient to make us account for the reaction mechanism. In 20,27 oxidation reactions. The second-order rate constant (k) general, MnO4 prefers to react with the electron-rich moiety
could be calculated by dividing the kp value by CH2O2. Obviously, through oxygen addition, hydrogen abstraction, and electron 20,27 this is the reaction basis for the potentiometric response of the transfer. For the electrophilic attack of H2O2 by MnO4 , the ionophore-free permanganate ISE under the zero-current primary reactive site is hard to be predicted. This is mainly
condition. Although its analytical working range was solely attributed to the structural simplicity of H2O2, thus hindering limited to 5–66 mM, which was not superior to that of the a complete disclosure of the reaction mechanism. titration method,21 the ionophore-free permanganate ISE together with the constructed theoretical model (eqn (6)) Theoretical modeling of the reaction mechanism offered a very promising alternative to study the rapid In addition to the experimental approach, the mechanism of Published on 26 December 2017. Downloaded by Library of Chinese Academy Sciences 9/11/2018 3:33:27 PM. the MnO4 /H2O2 reaction was also studied computationally. The potential energy prole is shown in Fig. 4, where the re- ported energy value is normalized using the reactants (i.e., MnO4 +H2O2) as the zero-point. The transition state (TS) was then obtained with the energy barriers of 24.5 kcal mol 1 from the zero-point and 21.0 kcal mol 1 from the reactant complex
Fig. 3 The plot of pseudo-first-order rate constant (kp) as a function C ¼ k Fig. 4 fi 1 of H2O2 at pH 7. Inset: the second-order rate constant ( ) measured The Gibbs free-energy pro le (kcal mol ) of the reaction C ¼ m at the pH equal to 6, 7, and 8 with H2O2 33 M. between MnO4 and H2O2.
470 | Anal. Methods,2018,10, 467–473 This journal is © The Royal Society of Chemistry 2018 View Article Online Paper Analytical Methods