Simulation-Based Approach to Determining Electron Transfer Rates

Apr 9, 2017 - Philippe Dauphin-Ducharme†‡∥ , Netzahualcóyotl Arroyo-Currás†‡∥ , Martin Kurnik†‡, Gabriel Ortega†‡§, Hui Li†‡, and Kevin W. Plaxco†‡...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Newcastle, Australia

Article

A Simulation-based Approach to Determining Electron Transfer Rates using Square-Wave Voltammetry Philippe Dauphin-Ducharme, Netzahualcoyotl Arroyo-Curras, Martin Kurnik, Gabriel Ortega, Hui Li, and Kevin W. Plaxco Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b00359 • Publication Date (Web): 09 Apr 2017 Downloaded from http://pubs.acs.org on April 11, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Langmuir is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

A Simulation-based Approach to Determining

1

2

Electron Transfer Rates using Square-Wave

3

Voltammetry

4

Philippe Dauphin-Ducharme1,2,‡, Netzahualcóyotl Arroyo-Currás1,2,‡, Martin Kurnik1,2, Gabriel

5

Ortega1,2,3, Hui Li1,2 and Kevin W. Plaxco1,2*

6

1

7

Barbara, CA 93106, USA.

8

2

9

USA.

Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa

Center for Bioengineering, University of California Santa Barbara, Santa Barbara, CA 93106,

10

3

11

‡These authors contributed equally.

12

*

CIC bioGUNE, Bizkaia Technology Park, Building 801 A, 48170 Derio, Spain

Corresponding author.

13

ACS Paragon Plus Environment

1

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 24

14

ABSTRACT

15

The efficiency with which square-wave voltammetry (SWV) differentiates faradaic and charging

16

currents renders it a particularly sensitive electroanalytical approach, as evidenced by its ability

17

to measure nanomolar or even picomolar concentrations of electroactive analytes. Due to the

18

relative complexity of the potential sweep it employs, however, the extraction of detailed kinetic

19

and mechanistic information from square-wave data remains challenging. In response we

20

demonstrate here a numerical approach by which square-wave data can be used to determine

21

electron transfer rates. Specifically, we have developed a numerical approach in which we model

22

the height and shape of voltammograms collected over a range of square-wave frequencies and

23

amplitudes to simulated voltammograms as functions of the heterogeneous rate constant and the

24

electron transfer coefficient. As validation of the approach we have employed it to determine

25

electron transfer kinetics in both freely diffusing and diffusion-less surface-tethered species,

26

obtaining electron transfer kinetics in all cases in good agreement with values derived using non-

27

square-wave methods.

28 29

KEYWORDS. Square-wave voltammetry, numerical simulations, electron transfer rate.

30 31 32

ACS Paragon Plus Environment

2

Page 3 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

33

MAIN TEXT. Square-wave voltammetry (SWV), an electroanalytical technique that involves

34

the application of repeating square-shaped potential pulses superimposed on a staircase potential

35

sweep, enables the effective discrimination of faradaic processes from charging currents.1 It

36

achieves this by imposing two opposite square-wave pulses of the same height at each staircase

37

potential step and measuring the resultant current at the end of each pulse, ifwd and ibwd. Taking

38

the difference between these currents returns a net voltammogram with improved signal-to-noise

39

ratios relative to those of, for example, cyclic voltammetry.1-2 This, in turn, renders SWV one of

40

the most sensitive electroanalytical techniques, enabling the detection of metal ions and organic

41

molecules at nanomolar and picomolar concentrations, respectively.1

42

In this work we build on prior approaches aimed at determining electron transfer kinetics using

43

SWV by modulating the height (amplitude) or width (frequency) of the square-wave pulses.1, 3

44

This earlier work includes that of Lovric4-5 who, for example, developed a mathematical method

45

termed “quasi-reversible maximum” that extracts electron transfer kinetics from the maximum in

46

a plot of SWV peak current versus square-wave frequency. Mirceski and co-workers,6-9 in

47

contrast, have developed a mathematical method that extracts electron transfer kinetics from the

48

relationship between SWV peak current and square-wave amplitude, with the maximum on this

49

plot again correlating with the electron transfer rate. These methods have allowed the

50

determination of the electron transfer rates of a range of freely diffusing8 and surface anchored6

51

redox couples, including linear and duplex DNA strands modified with a redox reporter10-11.

52

Historically, however, these approaches have been limited to the measurement of relatively slow

53

transfer rates of redox couples undergoing quasi-reversible or irreversible electron transfer

54

reactions.3, 8, 12 This leaves uncharacterized many important and fully reversible redox couples,

55

which transfer electrons more rapidly. Here we have expanded upon this prior literature and

ACS Paragon Plus Environment

3

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 24

56

developed a method of extracting electron transfer rates from SWV data that overcomes this

57

limitation, enabling access to much more rapid electron transfer kinetics than is possible with the

58

prior approaches.

59

The methods of Lovric and Mirceski rely on the fact that SWV peak height is defined by the

60

interplay between the forward and backward electron transfer rates and the frequency and

61

amplitude of the square-wave pulse.1 Specifically, Lovric’s approach5 extracts rates using the

62

relationship between square-wave frequency and the height of the resulting voltammetric peaks.

63

This works because changing the square-wave frequency optimizes the time scale of the

64

measurement to the time scale of electron transfer and, in turn, influences the peak height.

65

Mirceski’s approach,8 in contrast, extracts kinetic information from the relationship between

66

square-wave amplitude and the height of the resulting voltammetric peaks. This works because

67

changing the square-wave amplitude alters the driving force underlying electron transfer,

68

changing the rates of the forward and backward electron transfer reactions and thus altering peak

69

splitting and, in turn, peak height.

70

The success of the methods of Lovric and Mirceski demonstrates that information regarding

71

electron transfer kinetics can be extracted by observing peak height as a function of either

72

square-wave frequency or square-wave amplitude. Building on this we describe here an approach

73

that captures both effects simultaneously, potentially providing more information regarding rates

74

than methods that extract rates by varying only one or the other. Specifically, we present a

75

numerical approach that, by simultaneously modeling both the heights and shapes of

76

voltammetric peaks, is able to measure electron transfer kinetics for reactions that had proven too

77

rapid to measure using prior SWV methods.

78

ACS Paragon Plus Environment

4

Page 5 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

79

80 81

Figure 1. Description of the numerical model. Our approach starts with one-dimensional numerical models;

82

shown are models for (A) a freely diffusing redox couple and (B) a diffusion-less, surface-bound redox couple. (C)

83

We then simulated the response of these to SWV using a square waveform of frequency, f, and amplitude, ESW. This

84

is superimposed on a staircase that ramps the mean potential (the mean potential per pulse pair, E) by a potential

85

step size per cycle, Estep. In response to this varying potential the system produces an oscillating faradaic current

86

(grey trace in D). We deconvolute this into “forward” and “backward” voltammograms (red and black curves in D)

87

by extracting the current at specific times (ifwd and ibwd). The difference between these two voltammograms then

88

returns a “net” voltammogram (blue curve in D).

89 90 91

Our approach to extracting electron transfer kinetics uses numerical simulations to model the

92

height and shape of experimental SWV scans collected over a range of frequencies and

ACS Paragon Plus Environment

5

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 24

93

amplitudes. To do so we first numerically model the experimental set-up as a simple, one-

94

dimensional electrochemical cell. As examples we explore here two such cases: 1) a freely

95

diffusing redox couple, and 2) a redox couple that is confined to the electrode surface via a short,

96

flexible linker. We define two boundaries in each model, one representing the surface of the

97

working electrode and the other representing either the bulk of the solution (Figure 1A) or the

98

cell wall (Figure 1B). The reaction considered in all cases is reversible and defined as:  +   ⇌ 

Eq. 1

99

where n represents the number of electrons transferred. The electrode boundary condition for

100

both oxidized, O, and reduced, R, species is set in all cases to “flux” as defined by Fick’s second

101

law of diffusion: 

= ∇

   =  ∇  

Eq. 2

Eq. 3

102

where DO and DR are experimentally determined diffusion coefficients (for ferrocenemethanol,

103

for example, which we employ below, both DO and DR are 7.8 x 10-6 cm s-1)13 and CO and CR

104

represent the concentration of the oxidized and reduced species, respectively, and are both

105

functions of time, t, and distance from the electrode surface, x. Per Eq. 1, the initial conditions

106

are CO = CO,initial* and CR = 0 and the right boundary conditions are set as follows:

107

Case 1: diffusive redox couple ∗ lim→  (, ) = , ,

Eq. 4

lim→ " (, ) = 0 108

ACS Paragon Plus Environment

6

Page 7 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

109 110

Langmuir

Case 2: surface-confined redox couple lim→  (, ) =  , lim→ " (, ) = "

Eq. 5

111

Eq. 4 describes the case of the bulk boundary that is positioned far from the electrode surface so

112

that CO = CO, initial* (bulk concentration of oxidized species) and where CR = 0 remains constant at

113

all t and is unaffected by the electrode reaction. Eq. 5, in contrast, describes the case of a thin-

114

layer cell with an inert wall, which mimics the physics of an electroactive monolayer on an

115

electrode surface. It does so by limiting the length of the electrochemical cell to x