Article pubs.acs.org/JPCA
Kinetic Effects of the Complexation Reaction in the Facilitated Ion Transfer at Liquid Membrane Systems of One and Two Polarized Interfaces. Theoretical Insights Á ngela Molina,*,† Encarnación Torralba,† Carmen Serna,† and Joaquín A. Ortuño‡ †
Departamento de Química Física and ‡Departamento de Química Analítica Facultad de Química, Universidad de Murcia, 30100 Murcia, Spain ABSTRACT: An in-depth study of the ion transfer facilitated by complexation in the organic phase (TOC mechanism) in liquid membrane systems of one and two polarized interfaces is carried out by taking into account the kinetic effects associated with the complexation reaction. Explicit analytical equations for the normal pulse voltammetric (I/E) and chronoamperometric (I/t) responses with an explicit dependence on the kinetic parameters of the chemical complexation are presented for both kinds of membrane system, which could be useful for modeling artificial and biological membranes. The equations are compared with those obtained by using the widely used approximation of total equilibrium conditions that leads to the transfer by interfacial complexation mechanism (TIC), which only depends on thermodynamic parameters. Simple methods are proposed that allow quantitative determination of the equilibrium and kinetic constants of the complexation reaction in the organic phase for both kinds of membrane system.
1. INTRODUCTION Ion transfer across liquid−liquid interfaces is a fundamental physicochemical process with important implications in different areas, such as biology, chemistry and technology. It has been the subject of both theoretical and experimental studies.1−5 In cases such as the transfer of highly hydrophilic ions from water to an organic phase, ion transfer can be facilitated by a complexing agent, which normally is present in the organic phase.6−11 This facilitated transfer is widely used in nature, for example, in the uptake of potassium ion in cell membranes by selective transport assisted by the natural antibiotic valynomicin,12 and also in artificial membranes for ion-selective sensors and metal ion separations.6−14 The electrochemical investigation of the facilitated ion transfer reaction at ITIES (interface between two immiscible electrolyte solutions) was pioneered in 1979 by Koryta, who studied the transfer of alkali metal ions facilitated by different natural and synthetic ionophores.15 Since the work of Koryta, extensive studies of facilitated ion transfer have been carried out, many of them focused on the elucidation of the mechanism that this process follows which, even today, is a matter of controversy. Currently, it is accepted that facilitated ion transfer can take place via three mechanisms: through a TOC (transfer followed by complexation in the organic phase) or EC mechanism, where the ion transfer is the electrochemical step and the subsequent complexation the chemical one; through an ACT (aqueous complexation followed by transfer) or CE mechanism, in which the ligand is partitioned from an organic to an aqueous phase, and the complexation occurs before the ion transfer; and through a TIC/TID (transfer by interfacial © 2012 American Chemical Society
complexation/decomplexation), which involves a single ion transfer-ion complexation step occurring at the liquid/liquid interface.7,8 The active mechanism depends on many factors, such as the relative concentration of both the ion and the ligand in the two phases and the stability constant of the complex. In many references on facilitated ion transfer, the TIC mechanism is considered and the mathematical approaches used are based on pure thermodynamic principles, which implies that the kinetic effects of the chemical reaction are not considered. Moreover, in most of these references, the solution for the problem has been carried out with numerical methods in cyclic voltammetry, leading to complex integral solutions. In this work we analyze the facilitated ion transfer followed by organic complexation (TOC or EC mechanism) by taking into account the kinetic of this complexation reaction. We have obtained simple explicit analytical expressions for the flux or current corresponding to a TOC mechanism as a function of the ratio of the interfacial concentrations of the target ion between the aqueous and organic phases or as a function of the applied potential. These equations can be applicable to the electrochemical techniques normal pulse voltammetry (I/E) and chronoamperometry (I/t) in liquid membrane systems of one of two polarized interfaces, be they artificial membrane systems such as supported liquid membranes, solvent polymeric membranes, and bilayers or biological membranes. Both one Special Issue: Herman P. van Leeuwen Festschrift Received: November 14, 2011 Revised: January 4, 2012 Published: January 5, 2012 6452
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
and two polarized interfaces systems consist of an organic membrane that separates two aqueous solutions, but in the former, the polarization phenomena is only effective at one of the two liquid/liquid interfaces involved (the sample solution/ membrane interface or outer interface), whereas in the latter, both interfaces are polarized.1,2,16−22 We have also solved this problem, by following the treatment proposed by Matsuda et al.6 (which does not consider the kinetic of the complexation reaction), in both kinds of membrane system, and a comparison between the responses obtained by following both methodologies has been made. The importance of taking into account the kinetic of the chemical complexation and the fact that the TIC mechanism can be obtained as a limiting case of the TOC or EC mechanism is shown. Simple methods are proposed for the quantitative determination of the homogeneous equilibrium and kinetic constants of the complexation reaction in the membrane based on the evolution with time of the shift of the half-wave potential of the facilitated ion transfer.
Scheme 1. Schematic View of the Facilitated Ion Transfer via TOC or EC Mechanism (Solid Arrows), and via TIC or E Mechanism (Dashed Arrows) in a System of a Single Liquid/ Liquid Polarized Interface
complexation is the rate determining step.8 Note that coupled to the transfer of the target cation X+ through the outer interface (W1/M interface), a common ion present in the supporting electrolyte of the membrane and inner aqueous phase (W2) is transferred through the inner interface (M/W2 interface) to ensure electroneutrality, but as this ion is part of the supporting electrolytes of both phases, the inner interface is non polarizable and the system constitutes a system of a single polarized interface.20 For the sake of simplicity we assume that the bulk concentration of L in the membrane is large enough compared with that of X+ in the same phase for the condition of pseudofirst-order kinetic with respect to X+(M) to be satisfied. When a constant potential, E, is applied to this system and the uptake process establishes a fixed ratio of concentrations of the species involved on both sides of the interfacial region, mass transport can be described by the following diffusive-kinetic equation system:
2. THEORY 2.1. One Polarized Interface System. 2.1.1. Kinetic Effects of the Complexation Reaction (TOC or EC Mechanism). Let us consider the reversible ion transfer of a metal ion X+ from an aqueous solution (W1) to an organic one (M), which is facilitated by a highly lipophilic neutral ionophore, L, present only in the organic phase, according to the process described by the solid arrows in Scheme 1 where k1 and k2 are the forward and backward rate constants of the complexation reaction between the metal X+ and the ligand L to form the complex XL+ in the organic phase, and this
⎫ ⎪ ⎪ ∂t ∂x 2 ⎪ ⎪ ∂c XM+(x ,t ) ∂ 2c XM+(x ,t ) M (x ,t ) + k c M (x ,t ) ⎬ = D XM+ − k c + + 1 2 X XL ⎪ ∂t ∂x 2 ⎪ M (x ,t ) M (x ,t ) ⎪ ∂cXL ∂ 2cXL + + M M (x ,t ) − k c M (x ,t )⎪ = D XL + k c + + + 1 X 2 XL ⎭ ∂t ∂x 2
∂c XW+1(x ,t )
=
D XW+1
∂ 2c XW+1(x ,t )
c XM+(0,t ) = c XW+1(0,t )eη
with the boundary value problem ⎫ ⎬ t ≥ 0, x → −∞⎭ t = 0, x ≤ 0
⎫ t = 0, x ≥ 0 ⎪ ⎬ t ≥ 0, x → ∞⎪ ⎭
c XW+1(x ,0)
=
c XW+1(−∞,t )
= c X*+
η= (3)
⎛ ∂c W1(x ,t ) ⎞ X+ ⎟ ⎟ x ∂ ⎠ ⎝
=
⎛ ∂c M+(x ,t ) ⎞ X ⎟⎟ ∂x ⎠x = 0 ⎝
D XM+⎜⎜
x=0
⎛ ∂c M +(x ,t ) ⎞ M ⎜ XL ⎟⎟ D XL =0 +⎜ ∂x ⎝ ⎠x = 0
F W 0′ ϕ X+) (E − ΔM RT
(7)
where cip(x,t) and Dip are, respectively, the concentration of the species i in phase p (i = X+, XL+; p = W, M) and its diffusion W 0+ coefficient, and cX*+ and ΔM ϕX′ are the initial concentration of + the ion X and its formal ion transfer potential. Other symbols have their usual meaning. For simplicity, we assume that the diffusion coefficients of the free ion X+ and the complex XL+ in the organic phase are equal, i.e.
x = 0, t > 0 D XW+1⎜ ⎜
(6)
with
(2)
c XM+(x ,0) = c XM+(∞,t ) = 0 M (x ,0) = c M (∞,t ) = 0 cXL + XL+
(1)
(4)
M = DM D XM+ = D XL +
(5) 6453
(8)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
with χ being the dimensionless kinetic constant, given by
The current obtained is given by
I=
⎛ ∂c W1(x ,t ) ⎞ X+ ⎟ ⎟ ∂ x ⎠ ⎝
χ = κt
− FAD XW+1⎜ ⎜
κ is the sum of both kinetic constants involved in the chemical complexation, the forward pseudo-first-order kinetic constant (k1) and the backward first-order kinetic one (k2).
(9)
x=0
To solve the problem given by eqs 1−7, we introduce two variables M (x ,t ) ζ(x ,t ) = c XM+(x ,t ) + cXL +
κ = k1 + k 2
(10)
k1′c L*+(M) k2
c * +(M) k = 1 = XL k2 c X*+(M)
δr =
Id = FA
(12)
2 χ
F (x ) =
(1 + γeη(1 + K ))
(14)
c XW+1(x ,t ) c X*+ c XM+(x ,t ) c X*+
=
M (x ,t ) cXL +
c X*+
c X*+
=
(18)
(19)
⎛x⎞ π (x /2)2 xe erfc⎜ ⎟ ⎝2⎠ 2
⎞ ⎛ ⎛ ⎞ eη ⎜erfc⎜ x ⎟ + K δr e−x / δr⎟ ⎟ δ 1 + K δr /δ + γeη(1 + K ) ⎜⎝ ⎜⎝ 2 D Mt ⎟⎠ ⎠ ⎞ ⎛ ⎛ ⎞ K eη ⎜erfc⎜ x ⎟ − δr e−x / δr⎟ ⎟ ⎜ η ⎟ 1 + K δr /δ + γe (1 + K ) ⎜⎝ ⎝ 2 D Mt ⎠ δ ⎠
Idkss γeη(1 + K ) = Id 1 + K δr /δ + γeη(1 + K )
(20)
⎫ ⎪ (a) ⎪ ⎪ ⎪ ⎪ ⎪ ⎬ (b)⎪ ⎪ ⎪ ⎪ (c) ⎪ ⎪ ⎭
(21)
In this treatment, the kinetic effects of the complexation reaction are contained in the ratio δr/δ (=1/(πχ)1/2). This response can be linearized in a formally identical way to that corresponding to a simple ion transfer process, such that we have the following linearized potential/current relation for the facilitated ion transfer
(22)
with δr being the reaction layer thickness, given by eq 17 and δ the diffusion layer thickness, given by
πDMt
πt
⎛ ⎞ ⎜ ⎟ γeη(1 + K ) x =1− erfc⎜ ⎟ η W 1 + K δr /δ + γe (1 + K ) ⎜ 2 D +1t ⎟ X ⎝ ⎠
and
δ=
D XW+
When the kinetic of the coupled chemical reaction is, as usual, even faster, the diffusive kinetic steady state approximation (dkss) applies.23−25 The use of this approximation is highly desirable, because it leads to very useful, intuitive, and simple expressions for the current, the half-wave potential, interfacial concentrations and concentration profiles for all the species involved in the facilitated ion transfer process. Under this approximation, by following the procedure indicated in Appendix A.2, we can derive the following expressions for the concentration profiles and the current corresponding to the facilitated ion transfer in this kind of membrane system
(13)
K 2 DM t = (1 + γeη(1 + K )) K δr
(17)
⎛ M ⎞1/2 D γ=⎜ W⎟ ⎜ D 1⎟ ⎝ X+ ⎠
where
χ EC =
DM /κ
and
cX*+(M), cL*(M), and cXL * (M) being the equilibrium concentrations at phase M. Considering that the kinetic of the homogeneous chemical reaction are fast enough ((k1 + k2)t ≥ 15), the kinetic steady state solution (kss) applies, which relies on the perturbation of the chemical equilibrium being independent of time (∂φ(x,t)/∂t) = 0.23,24 Under these conditions, we can derive the following expression for the current corresponding to the facilitated ion transfer in a system of a single polarized interface (Appendix A.1)
Ikss γeη(1 + K ) = F(χ EC) Id 1 + γeη(1 + K )
(16)
δr is the reaction layer thickness, given by eqs A1.7 and A1.8
M (x ,t ) φ(x ,t ) = Kc XM+(x ,t ) − cXL + (11) where φ(x,t) refers to the distortion of the chemical equilibrium, and K is the stability constant of the complex under the pseudo-first-order assumption
K = K ′c L*+(M) =
(15)
E = E1/2 +
(23) 6454
RT ⎛ Idkss/Id ⎞ ln⎜ ⎟ F ⎝ 1 − Idkss/Id ⎠
(24)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
with
Under total equilibrium conditions the expression for the current associated with the facilitated ion transfer takes the following form (Appendix A.3)
RT ⎛ 1 ⎞ ln⎜ ⎟ F ⎝γ⎠ RT ⎛ 1 + K δr /δ ⎞ ⎟ + ln⎜ F ⎝ 1+K ⎠
W 0′ E1/2 = ΔM ϕ X+ +
Ite γeη(1 + K ) = Id 1 + γeη(1 + K ) (25)
where, as can be seen, the only dependence with the complexation reaction is through the thermodynamic complexation constant, and kinetic effects are not taken into account. This equation can be linearized, obtaining an expression for the potential difference analogous to eq 24, where the half-wave potential is given by
being the half-wave potential corresponding to the process. In conditions of inert equilibrium (K → 0) and very slow kinetic (δ → δr), eqs 21, 22, and 25 for the concentration profiles, current response, and half-wave potential corresponding to the facilitated ion transfer simplify to those corresponding to the simple ion transfer in this kind of membrane system20 c XW+(x ,t ) c X*+ c XM+(x ,t ) c X*+
⎛ γeη x ⎜ erfc⎜ =1− ⎜ 2 D W+t 1 + γeη X ⎝ =
⎛ x ⎞ eη ⎜ ⎟⎟ η erfc⎜ 1 + γe ⎝ 2 DMt ⎠
⎞⎫ ⎟⎪ ⎟⎟ ⎪ ⎠⎪ ⎬ ⎪ ⎪ ⎪ ⎭
W 0′ E1/2 = ΔM ϕ X+ +
RT ⎛ 1 ⎞ ln⎜ ⎟ F ⎝γ⎠ RT ⎛ 1 ⎞ ⎟ + ln⎜ F ⎝1 + K ⎠
W 0′ E1/2 = ΔM ϕ X+ +
(27)
RT ⎛ 1 ⎞ ln⎜ ⎟ F ⎝γ⎠
(28)
with η given by eq 7. 2.1.2. Total Equilibrium Conditions (TIC Mechanism). Equations 22 and 24 are strongly simplified if we suppose that the rates of the complex formation and dissociation processes are sufficiently large in comparison with the diffusion rate, so it may be assumed that the complex formation and dissociation are at equilibrium even when current is flowing, i.e., (k1 + k2)t ≫ 1 or δr → 0. This situation, similar to that described by Matsuda et al. in ref 6, corresponds with that called transfer by interfacial complexation in facilitated ion transfer processes, and it involves a single step occurring at the liquid/liquid interface (see also the ion transfer process described by the dashed arrows in Scheme 1)
X+(W) + L(M) ⇄ XL+(M)
c XW+1(x ,t ) c X*+ c XM+(x ,t ) c X*+ M (x ,t ) cXL +
c X*+ c RW+2(x ,t ) c X*+
⎛ ⎞ ⎜ ⎟ x = 1 − h(η) erfc⎜ ⎟ ⎜ 2 D W+1t ⎟ X ⎝ ⎠ =
⎛ h(η) ⎜ ⎛⎜ x ⎞⎟ + erfc γ(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠
=
⎛ Kh(η) ⎜ ⎛⎜ x ⎞⎟ − erfc γ(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠
=
D XW+1 DRW+2
(30)
which is an equivalent of the Deford−Hume equation deduced for planar diffusion.26 2.2. Two Polarized Interfaces System. 2.2.1. Kinetic Effects of the Complexation Reaction (TOC or EC Mechanism). We consider now the facilitated ion transfer in a membrane system of two polarized interfaces (process described by the solid arrows in Scheme 2). The simple ion transfer process in this kind of membrane system has been quite widely studied.16,17,19 Recently, simple and manageable expressions for the concentration profiles, current response, membrane potential, potential differences at the two interfaces involved, and half-wave potential have been given and experimentally verified by our research group,.22,27−29 When any constant potential, E, is applied to study the facilitated ion transfer in this kind of membrane system, we need to solve two separate problems of mass transport: that corresponding to the facilitated ion transfer of the ion, X+, at the outer interface (W1/M interface) and that of the ion that is simultaneously transferred at the inner one (M/W2 interface) to ensure electroneutrality.22 So, if we assume that coupled to the transfer of the metal ion through the outer interface, the ion of the supporting electrolyte of the membrane, R+, is transferred through the inner one, by means of the procedure indicated in Appendix B, we can derive the following response for the concentration profiles and the current associated to the facilitated ion transfer at systems of two polarized interfaces (TPI)
(26)
Idkss γeη = Id 1 + γeη
(29)
⎛ x−d ⎞ ⎟⎟ h(η) erfc⎜⎜ ⎝ 2 D W2t ⎠
6455
⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎞⎪ K δ r −x / δ r ⎟ e ⎟⎪ δ ⎠⎪ ⎬ ⎞ ⎪ δ r −x / δ r ⎟ ⎪ e ⎟ ⎪ δ ⎠ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭
(31)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
Scheme 2. Schematic View of the Facilitated Ion Transfer via TOC or EC Mechanism (Solid Arrows), and via TIC or E Mechanism (Dashed Arrows) in a System of Two Liquid/ Liquid Polarized Interfaces
and
ITPI,dkss = Idh(η)
(33)
(41)
(42)
(43)
with (34)
h′(η) =
1 {[(λ(1 + K ))2 e2η + 8λ(1 + K )eη]1/2 4 − λ(1 + K )eη}
(35)
(44)
and
with 0′ w 0′ w 0′ EM = ΔM ϕ X+ − ΔM ϕ R+
(36)
E=
W
w +0 DR+ 2, cR*+, and ΔM ϕR ′ are respectively the diffusion coefficient in the inner aqueous solution W2, the initial concentration, and the standard ion transfer potential of the ion R+. Equation 32 can be rewritten as follows
E=
1 {[(λeη)2 + 8λeη]1/2 − λeη} 4
ITPI,te = Idh′(η)
c R*+
F 0′ (E − EM ) RT
RT ⎛⎜ 1 ⎞⎟ ln F ⎝λ⎠
with λ given by eq 34 and η by eq 35. 2.2.2. Total Equilibrium Conditions (TIC Mechanism). If we assume the approximation of total equilibrium (see eq A3.1 and process described by dashed arrows in Scheme 2), by following the same methodology explained in Appendix B, but taking into account the solutions for the outer interface obtained under total equilibrium conditions (eqs 24 and 30), the following expression for the I/E curve and the potential difference in the system of two polarized interfaces are obtained
with K, δr, and δ given by eqs 12, 17, and 23, respectively, and λ and η by
η=
(40)
g(η) =
⎧⎡ 2 ⎤1/2 λ(1 + K ) η⎥ 1 ⎪⎢⎛ λ(1 + K ) ⎞ 2η e h(η) = ⎨ ⎜ ⎟ e +8 4 ⎪⎢⎣⎝ 1 + K δr /δ ⎠ 1 + K δr /δ ⎥⎦ ⎩
c X*+
ITPI = Id g(η)
with g(η) given by
where h(η) is the function that contains the dependence with the applied potential
2 DRW+2DM λ= D XW+1
(39)
1/2 0′ EM = EM +
(32)
⎫ λ(1 + K ) η⎪ − e⎬ 1 + K δr / δ ⎪ ⎭
⎛ ⎞⎫ ⎜ ⎟⎪ x = 1 − g (η) erfc⎜ ⎟⎪ W c X*+ ⎜ 2 D +1t ⎟ ⎪ X ⎝ ⎠⎪ ⎪ ⎛ ⎞ ⎪ c XM+(x ,t ) g (η) ⎜ ⎛ x ⎞⎟ ⎬ ⎟⎟ erfc⎜⎜ = ⎜ ⎟ ⎪ c X*+ γ ⎝ ⎝ 2 D Mt ⎠⎠ ⎪ ⎪ W2 W1 ⎛ x − d ⎞⎪ c R+ (x ,t ) D X+ ⎟⎟ ⎪ g (η) erfc⎜⎜ = c X*+ D RW+2 ⎝ 2 D W2t ⎠ ⎪ ⎭ c XW+1(x , t )
1/2 EM
2⎞ ⎛ RT ⎜ 2(ITPI,dkss/Id) ⎟ + ln⎜ F ⎝ 1 − ITPI,dkss/Id ⎟⎠
2⎞ ⎛ RT ⎜ 2(ITPI,te/Id) ⎟ + ln F ⎜⎝ 1 − ITPI,te/Id ⎟⎠
(45)
with
RT ⎛⎜ 1 ⎞⎟ RT ⎛⎜ 1 ⎞⎟ 1/2 ′ 0′ EM = EM + + ln ln F ⎝λ⎠ F ⎝1 + K ⎠
(46)
which, as can be seen, coincides with eqs 32, 33, 37, and 38 when δr approaches to zero. Equation 46 constitutes an equivalent of the Deford−Hume equation for planar diffusion in a system of two liquid/liquid polarized interfaces. Concentration profiles and interfacial concentrations of the species involved under te conditions are given in Appendix B.
(37)
1/2 with EM being the half-wave membrane potential corresponding to the facilitated ion transfer in this kind of membrane system
RT ⎛⎜ 1 ⎞⎟ ln F ⎝λ⎠ RT ⎛ 1 + K δr /δ ⎞ ⎟ + ln⎜ F ⎝ 1+K ⎠
1/2 ′ EM
1/2 0′ EM = EM +
3. RESULTS AND DISCUSSION Figure 1 shows the I/E curves corresponding to the facilitated ion transfer in a system of a single polarized interface obtained by using the kss solution (solid lines, eq 13), corresponding to a fixed value of the dimensionless kinetic constant χ and different values of the stability constant of the complex K (Figure 1a), and to a fixed value of the stability constant and different values
(38)
When K → 0 and δ → δr, eqs 31, 32, and 38 are transformed into those corresponding to a simple ion transfer in this type of membrane system22 6456
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
Figure 2. I/E curves of the ion transfer facilitated by complexation in the organic phase obtained by using the kss (solid lines, eq 13), dkss (dashed lines, eq 22), and te (dash-dotted lines, eq 29) solutions in a liquid membrane system of a single polarized interface, for a fixed value of χ and w 0+ ϕX′ = 0 mV, cX*+ = different values K (a)−(c), and vice versa (d)−(f). ΔM W+1 −5 2 −1 M+ −8 2 −1 1 mM, DX = 10 cm s , DX = 10 cm s , t = 1 s, A = 0.081 cm2, T = 298.15 K, χ and K values are given on the graphs.
Figure 1. Solid lines: I/E curves corresponding to the ion transfer facilitated by organic complexation in a system of a single polarized interface obtained using the kss solution (eq 13) for a fixed value of χ and different values of K (given on the curves, (a)), and for a fixed value of K and different values of χ (b). Dotted lines: I/E curves of the simple ion transfer (eq 27). ΔMw ϕX0′+ = 0 mV, cX*+ = 1 mM, DXW+1 = 10−5 cm2 s−1, DXM+ = 10−8 cm2 s−1, t = 1 s, A = 0.081 cm2, T = 298.15 K.
(dash-dotted lines), obtained from eqs 13, 22, and 29, respectively, at different values of the stability constant K (Figure 2a−c) and the dimensionless kinetic constant χ (Figure 2d−f). As can be seen from Figure 2a−c, as K increases, the I/E curves corresponding to the dkss and te approximations become less coincident with that corresponding to the more rigorous kss solution, with the dkss curves displaced toward more positive potential values and the te ones toward more negative potentials. This displacement is more important in the case of the te approximation, in which curves show a shift ≥50 mV with respect to the kss curves at potential values near the half-wave potential for a stability constant value of K ≥ 50, whereas the shift shown by the dkss curves is around 8 mV. So, when the chemical equilibrium is displaced toward the complex, the dkss approximation is more suitable than the te one for characterizing the system. From Figure 2d−f it can be noted that, for a very high value of χ (Figure 2f), all the I/E curves corresponding to the facilitated ion transfer overlap, and any of the approximations presented is suitable. However, as χ decreases, the dkss curves and the te ones differ from the kss curves. Thus, for χ = 15 the dkss and te curves are displaced about 8 and 30 mV, respectively, with respect to the kss curve at potential values closer to the half-wave potential of the system. So, from these drafts we can conclude that the dkss approximation can be confidently used for high stability constant values and relatively small values of the dimensionless kinetic constant, the te approximation not being acceptable under these conditions. For other conditions, the kss
of the dimensionless kinetic one (Figure 1b). In the dotted lines the I/E curve corresponding to a simple ion transfer obtained from eq 27 is presented. A value of the diffusion coefficient in the membrane typical of plasticized polymeric membranes has been used.27,29 With respect to the influence of the stability constant K, from Figure 1a it can be seen that, for a fixed value of the dimensionless kinetic constant χ, the I/E curves shift toward less positive potential values with respect to those corresponding to the simple ion transfer as K increases. The more stable the complex, the more displaced the I/E curve. Regarding the influence of χ, in Figure 1b it is shown that for a fixed value of the stability constant, the ion transfer becomes easier, i.e., the I/E curves shift toward less positive potential values with respect to the curve corresponding to the simple ion transfer, as χ (the metal-complex interconversion) increases. It is worth noting that the limiting current is the same for all the curves, independently of the K or χ value (Figure 1a,b), so the facilitated ion transfer cannot be differentiated from the simple ion one at limiting current potentials, because the current depends only on the concentration of the metal in the aqueous phase, as will be shown in Figure 3. In Figure 2 we have plotted the normal pulse voltammograms (I/E curves) corresponding to the facilitated ion transfer in a system of a single polarized interface obtained by using the kss solution (solid lines), the dkss approximation (dashed lines), and the te one 6457
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
system of a single polarized interface in limiting current conditions (Figure 3a), and for a potential applied of 50 mV (Figure 3b), for two different values of the stability constant (shown on the curves). The dotted lines correspond to the concentrations profiles of a simple ion transfer (eq 26) in the same conditions. Figure 3c shows the I/E curves corresponding to the facilitated ion transfer at the two stability constants selected for Figure 3a,b obtained by the dkss solution (solid lines, eq 22), together with those obtained with the te approximation (dashed lines, eq 29), and that corresponding to a simple ion transfer process (dotted lines, eq 27). A vertical dashed line shows the potential chosen for Figure 3b. As can be seen from Figure 3a, when a very high positive potential is applied, the profile corresponding to the facilitated ion transfer of the metal X+ in the aqueous phase coincides with that corresponding to a simple ion transfer in the same conditions. So, the current flux of the facilitated ion transfer in this case is independent of the kinetic of the complexation reaction. This feature is not observed for the metal X+ in the membrane, because part of this ion is lost in its interconversión into the species XL+. The more stable the complex, the less X+ remains free in the membrane. For an intermediate applied potential, the concentration profile corresponding to the metal X+ in the aqueous phase is well differentiated from that corresponding to a simple ion transfer (Figure 3b). Therefore, the current does depend on the kinetic of the chemical complexation in these conditions, and hence, these conditions will be desirable for studying the kinetic of the chemical reaction following the ion transfer. Moreover, from an experimental point of view working at not very high potentials is also preferable, because the background currents of the experiment and the deterioration of the experimental membrane system are expected to be less, and to have a better selectivity. It is necessary to clarify that a slope for the concentration profile of the target ion in the organic phase is 1000 times steeper than in the aqueous phase, and a flat profile for the complex XL+ near the interface is obtained in Figure 3a,b in agreement with that predicted by eqs 4 and 5; however, it cannot be evident at first sight due to the different scales used for the aqueous and organic phases. From Figure 3c it can be observed that, under the conditions selected, the I/E curves corresponding to K = 0.7 obtained by taking into account the kinetic of the coupled chemical reaction overlap with that obtained with the hypotheses of total equilibrium. This is logical, because at these small complexation constants the complexation reaction scarcely takes place and the kinetic do not play an important role. However, as the stability constant increases, the kinetic effects gain importance, and the te curves differ from the dkss ones, showing a displacement of approximately 8 mV for K = 20. The separation shown between the I/E curves of the te approximation and that corresponding to a simple ion transfer is given by (RT/F) ln(1/(1 + K)) (eqs 30 and 28). In Figure 4 we have represented the function perturbation of the chemical equilibrium φ(x,t) (eq 11) for the facilitated ion transfer in conditions similar to those discussed in Figure 3, i.e., for an applied potential corresponding to limiting current φlim(x,t) (dash-dot lines) and to the half-wave potential φ1/2(x,t) (solid lines). In these conditions the expression for the function φ(x,t) simplifies to the following
Figure 3. (a) and (b) Solid lines: Concentration profiles corresponding to the ion transfer facilitated by complexation in the organic phase in a membrane system of a single polarized interface (eq 21) for two different values of K (shown on the curves), obtained under limiting current conditions (a) and for a potential applied of 50 mV (b). Dotted lines: concentration profiles corresponding to a simple ion transfer (eq 26). (c) I/ E curves of the facilitated ion transfer obtained from the dkss and te approximations (solid and dash-dotted lines; eqs 22 and 29, respectively) at the same K values as used in (a) and (b) (shown on the curves), and the I/ E curve corresponding to a simple ion transfer process (dotted line, eq 27). w 0+ ϕX′ = te and dkss curves corresponding to K = 0.7 are overlapped. ΔM W 0 mV, cX*+ = 1 mM, DX+1 = 10−5 cm2s −1, DXM+ = 10−8 cm2 s−1, t = 1s, χ = 500, A = 0.081 cm2, T = 298.15 K.
solution should be used. What is more, by regarding Figure 2 as a whole, it can be seen that for potentials near the foot of the wave, the dkss approximation can be used independently of the kinetic of the complexation reaction, i.e., independently of the K and χ values. Parts a and b of Figure 3 show, in solid lines, the concentration profiles obtained from eq 21 corresponding to the ion transfer facilitated by organic complexation in a membrane
c *+K δr /δ −x / δ r φlim(x ,t ) = 2φ1/2(x ,t ) = X e γ
(47)
As can be observed from Figure 4 and also from eq 47, under limiting current conditions the function of the perturbation of the 6458
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
Figure 4. ϕ(x,t)/x curves of the ion transfer facilitated by complexation in the organic phase in a system of a single polarized interface obtained from eq 11 at different time values (shown on the curves), in limiting current conditions (dash-dotted lines) and when the applied potential coincides with the half-wave potential (solid w 0+ ϕX′ = 0 mV, cX*+ = 1 mM, DXW+1 = 10−5 cm2s −1, DXM+ = lines). ΔM 10−8 cm2 s−1, K = 5, χ = 500, A = 0.081 cm2, T = 298.15 K.
chemical equilibrium takes values double those corresponding to the half-wave potential. For a given time value, φ(x,t) decreases as the distance to the interface increases as a consequence of the inefficiency of the diffusion transport to supply ion X+ to be complexed. The thickness of the linear reaction layer, which is given for the value in the abscissa axis of the cutoff point of the tangent line to the φ curve at the interface (Figure 4), remains constant independently of the duration of the potential pulse and of the potential amplitude, because it only depends on the diffusion coefficient and on the kinetic constants (eq 17). Figure 5a shows the variation of ΔE 1/2 with the dimensionless kinetic parameter χ for different values of the stability constant K (shown on the curves) for a system of a single polarized interface, with ΔE1/2 being the shift of the halfwave potential of the facilitated ion transfer (cL* ≠ 0) with respect to that corresponding to a simple ion transfer process (cL* = 0)25 (eq 48). These curves have been obtained by using the kss solution (solid lines, acquired numerically), the dkss approximation (dashed lines, eq 25), and the te approximation (dash-dotted lines, eq 30). In Figure 5b we present the plots exp(zFΔE1/2/RT) versus 1/√t for a fixed value of K and different values of κ (shown on the curves). For a given value of the stability constant, ΔE1/2 values from the kss and dkss solutions shift to more negative potentials as χ increases (Figure 5a), because the transfer is being favored by the increasingly fast metal-complex interconversion causing a decrease in the Gibbs energy of ion transfer The higher the equilibrium constant, the more stable the complex and the easier the ion transfer. Consequently, a more accentuated displacement to less positive potentials is observed as K increases. On the contrary, the ΔE1/2 obtained by using the te approximation is not affected by the kinetic. So the te treatment cannot be used for small values of χ and higher values K. For example, the deviation observed with the te approximation is approximately 23 mV for K = 10 and nearly 58 mV for K = 50 with respect to the more rigorous values obtained with the kss solution for a χ value of 50, whereas with the dkss aproximation this error is of approximately 8 mV, so the dkss
Figure 5. (a) Evolution of ΔE1/2 with χ (=κt, with t = 1 s) in a system of a single polarized interface for different values of K (shown on the curves), obtained by using the kss solution (solid lines, acquired numerically), the dkss approximation (dashed lines, eq 25), and the te approximation (dash-dotted lines, eq 30). (b) Linear plots of exp(zFΔE1/2/FT) versus 1/√t obtained from the kss approximation for the system of a single polarized interface for K = 50 and different W w 0+ ϕX′ = 0 mV, cX*+ = 1 mM, DX+1 = values of κ (shown on the curves). ΔM 10−5 cm2 s−1, DXM+ = 10−8 cm2 s−1, A = 0.081 cm2, T = 298.15 K.
implies a great improvement in the quality of the measurements with respect to the te aproximation under these conditions. For quantitative determination of the homogeneous kinetic constant of the chemical complexation, it is worth noting that the curves given in Figure 5a can be used as working curves, from which for a given time of the experiment, (κ = χ/t) can be immediately obtained once the stability constant K is known. The variation in the dimensionless kinetic constant can be attained by varying the time of the experiment and also by changing the ligand concentration (eq 12). The plots given in Figure 5b are very useful to the complete characterization of the coupled chemical complexation. The representation of exp(zFΔE1/2/FT) with 1/√t is linear, and the values of the slope and the intercept can be used to obtain κ and K of the complexation reaction, as can be seen in Figure 5b and also in eq 49
ΔE1/2 = E1/2(cL*≠0) − E1/2(cL*=0) RT ⎛ 1 + K / πχ ⎞ = ln⎜ ⎟ F ⎝ 1+K ⎠ 6459
(48)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
See eqs 25 and 28 for the dkss solution,
to more positive potentials than that of the single polarized interface one. However, as we have pointed out in previous papers,28,29 by varying the concentration of the supporting electrolyte of the membrane in the two polarized interfaces system, displacements of the I/E curve can be attained. This fact constitutes an additional advantage of the system of two polarized interfaces besides the well-known onesas the double wide of the potential window, or the huge separation between the signals of cations and anions with similar standard ion transfer potentials17,29 because it confers the signal of this system the possibility of being located in the potential window available. Moreover, as the equilibrium constant increases and the kinetic effects gain importance, the te curves start to differentiate from the dkss ones in both kinds of membrane system, and it is necessary to take into account the kinetic effects.
exp(zF ΔE1/2 /FT ) 1 K 1 = + 1+K (1 + K ) πκ t (49) Note that eq 49 is also applicable to characterize the coupled chemical complexation in liquid membrane systems of two polarized interfaces because, according to our equations, the kinetic dependence of the half-wave potential is the same for both kinds of membrane system (eqs 25 and 38). In Figure 6 we have plotted the I/E curves corresponding to the facilitated ion transfer in systems of one (Figure 6a) and two
4. CONCLUSIONS • Explicit analytical equations for the normal pulse voltammetric (I/E) and chronoamperometic (I/t) responses of the ion transfer facilitated by complexation in the organic phase (TOC mechanism) in liquid membrane systems of one and two polarized interfaces have been obtained by considering the kinetic effects of the chemical complexation. These expressions could be useful for modeling artificial and biological membranes. • The dependence of the ion transfer process on the kinetic parameters of the complexation reaction have been analyzed in-depth for both kinds of membrane system, and the results have been compared with those obtained by using the widely used approximation of total equilibrium conditions, for which the results only depend on thermodynamic parameters. • It is shown that, under certain conditions, the transfer by interfacial complexation mechanism (TIC), which involves a single ion transfer-ion complexation step can be obtained as a limiting case of the two steps comprising TOC or EC mechanism. The conditions under which the treatment of total equilibrium (TIC) is not acceptable and the importance of taking into account the kinetic effects of the complexation reaction have been shown. • Simple methods are proposed for the quantitative determination of the equilibrium and kinetic constants of the complexation reaction in the membrane, based on the evolution with time of the shift of the half-wave potential of the facilitated ion transfer.
Figure 6. I/E curves corresponding to the ion transfer facilitated by complexation in the organic phase in systems of one (a) and two (b) polarized interfaces, obtained by using the dkss solution (solid lines, eqs 22 and 32) and that of total equilibrium (dashed-dotted lines, eqs 29 and 43). Dotted lines show the I/E curves corresponding to the w 0+ w +0 ϕX′ = 0 mV, ΔM ϕR ′ = simple ion transfer process (eqs 27 and 40). ΔM W+1 W2 −5 2 −1 + + −100 mV, cX* = 1 mM, cR* = 50 mM, DX = DR+ = 10 cm s , DXM+ = 10−8 cm2 s−1, χ = 500, t = 1 s, A = 0.081 cm2, T = 298.15 K.
■
(Figure 6b) polarized interfaces, obtained by using the dkss aproximation (solid lines, eqs 22 and 32) and the te one (dasheddotted lines, eqs 29 and 43). In dotted lines, the I/E curves corresponding to a simple ion transfer in both kinds of membrane system obtained for eqs 27 and 40, respectively, have been plotted. Under the experimental conditions selected, the curves corresponding to the two polarized interfaces system are shifted
In the following we consider that X+ is only initially present in the sample aqueous solution because it is very hydrophilic (like the metal ions involved in facilitated ion transfer processes), and L is only initially present in the membrane because it is very hydrophobic, and that all the phases, W1, W2, and M, contain sufficient concentrations of electrolytes such that the different ohmic drops can be neglected. It is worth clarifying that, although Scheme 1 shows a finite thickness
APPENDIX A. ANALYTICAL SOLUTION CORRESPONDING TO THE PROBLEM OF MASS TRANSPORT OF THE ION TRANSFER FACILITATED BY ORGANIC COMPLEXATION IN SYSTEMS OF ONE POLARIZED INTERFACE
Kinetic Steady State (kss)
6460
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
∂ 2c XW+1(x ,t ) ⎫ ⎪ = ⎪ ∂t ∂x 2 ⎬ 2 ⎪ ∂ζ(x ,t ) ∂ ζ(x ,t ) ⎪ = DM 2 ⎭ ∂t ∂x
for the membrane (d), we consider this thickness to be large enough in comparison with that of the diffusion layers of X+ and XL+ in this phase to consider semi-infinite diffusion into the membrane. When the new variables ζ(x,t) and φ(x,t) (given by eqs 10 M+ and 11) were taken into account and DXM+ = DXL = DM was assumed, the diffusive-kinetic equations and boundary value problem given by eqs 1−6 are transformed into the following
(A1.1)
t = 0,x ≤ 0 ⎫ ⎬ t ≥ 0, x → −∞⎭
(A1.2)
D XW+1
t = 0, x ≥ 0 ⎫ ⎬ t ≥ 0, x → ∞⎪ ⎭ ⎪
t = 0, x ≤ 0
∂ 2c XW+1(x ,t )
t = 0, x ≥ 0 ⎫ ⎬ t ≥ 0, x → ∞⎭
(A1.3)
x=0
⎛ ∂ζ(x ,t ) ⎞ = DM⎜ ⎟ ⎝ ∂x ⎠x = 0
(A1.10)
(A1.11)
(A1.12)
DM ⎛ ∂ζ(x ,t ) ⎞ ⎜ ⎟ κ ⎝ ∂x ⎠x = 0
= (1 + K )eηc XW+1(0,t )
x = 0, t > 0
(A1.13)
Using a procedure based on Kouteckýs dimensionless parameter method,30 the expression given by eqs 13−20 is deduced for the current.
⎛ ∂ζ(x ,t ) ⎞ = DM⎜ ⎟ ⎝ ∂x ⎠x = 0
(A1.4)
⎛ ∂ζ(x ,t ) ⎞ ⎛ ∂φ(x ,t ) ⎞ ⎜ ⎟ = K⎜ ⎟ ⎝ ∂x ⎠x = 0 ⎝ ∂x ⎠x = 0
(A1.5)
ζ(0,t ) + φ(0,t ) = (1 + K )eηc XW+1(0,t )
(A1.6)
Diffusive-Kinetic Steady State (dkss)
In this approximation the variable φ retains the form given by eq A1.7 for the kinetic steady state approximation (∂φ(x,t)/∂t = 0) W and a purely diffusive behavior for the species cX1 + and the pseudo-species ζ is assumed; i.e., it is assumed that the solution of equations for cX1W+ and ζ have the same form as that for species that would only suffer planar diffusion
with η and κ being given by eqs 7 and eq 16. If we consider the kinetic steady state approximation, which implies that the perturbation of the chemical equilibrium is independent of time, i.e., ∂φ(x,t)/∂t = 0,23 we can write
⎡ κ φ(x) = φ(0) exp⎢ − ⎣ DM
ζ(x ,0) = ζ(∞,t ) = 0
⎛ ∂c W1(x ,t ) ⎞ X+ ⎟ ⎟ x ∂ ⎠ ⎝
ζ(0,t ) − K
φ(x ,0) = φ(∞,t ) = 0
x=0
c XW+1(x ,0) = c XW+1(−∞,t ) = c X*+
D XW+1⎜ ⎜
ζ(x ,0) = ζ(∞,t ) = 0
⎛ ∂c W1(x ,t ) ⎞ X+ ⎟ ⎟ x ∂ ⎠ ⎝
(A1.9)
x = 0, t > 0
c XW+1(x , 0) = c XW+1(−∞ , t ) = c X*+
D XW+1⎜ ⎜
D XW+
⎫ ⎬ t ≥ 0, x → −∞⎭
⎫ ⎪ = 2 ⎪ ∂t ∂x ⎪ 2 ⎪ ∂ζ(x ,t ) M ∂ ζ(x ,t ) ⎬ =D 2 ⎪ ∂t ∂x ⎪ 2 ⎪ ∂φ(x ,t ) ∂ φ ( x , t ) = DM − κφ(x ,t )⎪ 2 ⎭ ∂t ∂x ∂c XW+1(x ,t )
∂c XW+1(x ,t )
⎤ x⎥ ⎦
⎛ x ⎞ ⎟⎟ ζ(x ,t ) = ζ(0,t ) erfc⎜⎜ ⎝ 2 DMt ⎠
(A2.1)
c XW+1(x ,t ) = c X*+ + (c XW+1(0,t ) − c X*+) ⎛ ⎜ x × erfc⎜ ⎜ 2 D W+1t X ⎝
(A1.7)
and thus,
φ(∞) − φ(0) dφ(x) −φ(0) = = δr dx δr (A1.8) with δr being the reaction layer thickness, given in eq 17 When eq A1.7 is taken into account, the diffusive-kinetic equations and boundary value problem eqs A1.1−A1.6 become c XW+1(0,t ) c X*+ c XM+(0,t ) c X*+ M (0,t ) cXL +
c X*+
=
⎞ ⎟ ⎟ ⎟ ⎠
(A2.2)
By inserting eqs A2.1 and A2.2 into eqs A1.12 and A1.13, we can derive the expressions for the interfacial concentrations and concentration profiles for the species involved in the facilitated ion transfer, given by eqs A2.3 and 21
=
1 + K δr /δ 1 + K δr /δ + γeη(1 + K )
=
eη(1 + K δr /δ) 1 + K δr /δ + γeη(1 + K )
K eη(1 − δr /δ) 1 + K δr /δ + γeη(1 + K )
6461
⎫ (a) ⎪ ⎪ ⎪ ⎪ ⎪ (b)⎬ ⎪ ⎪ ⎪ (c) ⎪ ⎪ ⎭
(A2.3)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
and by substituting eq 21a into eq 9, we deduce the expression for the current for the facilitated ion transfer given by eq 22.
x = 0, t > 0 ⎛ ∂c W1(x ,t ) ⎞ ⎛ ∂c M(x ,t ) ⎞ Xt X+ ⎟ ⎟ = DM⎜⎜ ⎟ ⎟ x x ∂ ∂ ⎝ ⎠x = 0 ⎠x = 0 ⎝
Total Equilibrium Conditions (te)
D XW+1⎜ ⎜
If we assume that the complex formation and dissociation are at equilibrium even when current is flowing, i.e.
c +(x ,t ) K = XL c X+(x ,t )
∀ x, t > 0
Assuming that the interfacial concentrations for cXW+ and ctX M are independent of time, the solution of eqs A3.3 is immediate
(A3.1)
and we define the following variable6 M (x ,t ) c XMt(x ,t ) = c XM+(x ,t ) + cXL + (A3.2) Equations 1−5 corresponding to mass transport for the facilitated ion transfer at systems of a single polarized interface are transformed into
∂ 2c XW+1(x ,t ) ⎫ ⎪ 2 ⎪ ∂t ∂x ⎬ ∂c XMt(x ,t ) ∂ 2c XMt(x ,t ) ⎪ ⎪ = D XM+ ⎭ ∂t ∂x 2
∂c XW+1(x ,t )
c XW+1(x ,t ) = c X*+ + (c XW+1(0) − ⎛ ⎜ x × erfc⎜ ⎜ 2 D W+1t X ⎝
= D XW+
⎫ ⎬ t ≥ 0, x → −∞⎭ t = 0, x ≤ 0
t = 0, x ≥ 0 ⎫ ⎬ t ≥ 0, x → ∞⎭
c X*+) ⎞ ⎟ ⎟ ⎟ ⎠
(A3.7)
⎛ x ⎞ ⎟⎟ c XMt(x ,t ) = c XMt(0) erfc⎜⎜ ⎝ 2 DMt ⎠
(A3.8)
(A3.3)
c XW+1(x ,0) = c XW+1(−∞,t ) = c X*+ c XMt(x ,0) = c XMt(∞,t ) = 0
and by substituting eqs A3.7 and A3.8 in the condition of the equality of fluxes (eq A3.6), taking into account eqs A3.1, A3.2, and 6, the following concentration profiles and interfacial concentrations for the species involved in the facilitated ion transfer are obtained
(A3.4)
(A3.5)
⎛ ⎜ γeη(1 + K ) x =1− erfc⎜ η 1 + γe (1 + K ) ⎜ 2 D W+1t X ⎝
c XW+1(x ,t ) c X*+
c XM+(x ,t ) c X*+ M (x ,t ) cXL +
c X*+
=
=
⎞ ⎟ ⎟ ⎟ ⎠
⎛ ⎛ ⎞⎞ eη ⎜erfc⎜ x ⎟⎟ 1 + γeη(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠⎟⎠
⎛ ⎛ ⎞⎞ K eη ⎜erfc⎜ x ⎟⎟ 1 + γeη(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠⎟⎠
⎫ ⎪ (a) ⎪ ⎪ ⎪ ⎪ ⎪ ⎬ (b)⎪ ⎪ ⎪ ⎪ (c) ⎪ ⎪ ⎭
(A3.9)
As can be seen, these expressions coincide with those obtained by using the dkss approximation when δr → 0. The expression for the current corresponding to the facilitated ion transfer under these conditions is obtained by substituting eqs A3.7 and A3.10a into eq 9 and is given by eq 29.
and
⎫ 1 ⎪ = η c X*+ 1 + γe (1 + K ) ⎪ ⎪ (a) ⎪ ⎪ c XM+(0) ⎪ eη ⎪ = η c X*+ 1 + γe (1 + K ) ⎬ ⎪ (b) ⎪ ⎪ M (0) cXL + ⎪ K eη = ⎪ c X*+ 1 + γeη(1 + K ) ⎪ ⎪ (c) ⎭
(A3.6)
c XW+1(0)
■
APPENDIX B. ANALYTICAL SOLUTION CORRESPONDING TO THE PROBLEM OF MASS TRANSPORT OF THE ION TRANSFER FACILITATED BY ORGANIC COMPLEXATION IN SYSTEMS OF TWO POLARIZED INTERFACES
For simplicity, for the development of the analytical solution to this problem we will consider conditions of diffusive-kinetic steady state. When a constant potential, E, is applied to study the facilitated ion transfer given by Scheme 2 in a system of two polarized liquid−liquid interfaces,22,28 this applied potential is
(A3.10) 6462
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
equality of current fluxes at the outer and inner interfaces (I(outer interface) = I(inner interface)),22 it is easy to deduce the following expression for the potential difference at the inner interface, Einn (eq B.4)
distributed unequally between the outer and inner interfaces. Thus, the differential equations system and boundary value problem describing mass transport at phases W1, M, and the outer interface are identical to those given by eqs 1−7 on changing E for Eout and η for ηout. The differential equations system and boundary value problem describing mass transport at phase W2 and inner interface are given by
∂c RW+2 ∂t
=
∂ 2c RW+2 W2 DR+ 2
(B.1)
∂x
c RM+(x ,t ) = c R*+ t = 0, x > d ⎫ ⎬ t ≥ 0, x → ∞⎭
E
(B.2) c XW+2(x ,0) = c XW+2(∞,t ) = 0
inn
⎛ W2 RT ⎜ DR+ ln = + F ⎜⎜ D W+1 X ⎝ RT ⎛ I ⎞ ln⎜ ⎟ − F ⎝ Id ⎠ W 0′ ΔM ϕ R+
(B.6)
and by subtracting the expressions for the potential differences of the outer and inner interfaces (eqs 24 and B.6, respectively), we can deduce the relations between the potential applied and the measured current in this kind of membrane system, given by eqs 33 and 37. The interfacial concentrations of the ions involved in the facilitated ion transfer as a function of the membrane potential for the system of two polarized interfaces are given by
(B.3)
t > 0, x = d c RM+(d ,t ) = e ηinnc RW+2(d ,t ) (B.4) with F W 0′ ηinn = ϕ R+) (E inm − ΔM (B.5) RT where Einn is the potential difference at the inner interface (membrane/inner solution interface). In this case, besides the considerations made in Appendix A, we assume that neither the diffusion layer nor the reaction layer corresponding to the facilitate ion transfer of X+ overlap with the diffusion layer of R+ and, moreover, that the concentration of R+ is high enough to be considered constant in all the membrane, as this ion is part of the supporting electrolyte of the membrane and the transport is controlled by the diffusion of the ion X+ According to our model, the problem corresponding to the facilitated ion transfer of the ion X+ at the outer interface is identical to that we have already solved for the same ion transfer at the system of a single polarized interface. Therefore, the expression for the current at this interface is equivalent to that given by eq 22 by changing η for ηout and E for Eout. With respect to the problem of mass transport of the ion R+ at the inner interface (eqs B.1−B.5), its solution is relatively easy, because it consists of a simple reversible ion transfer, in such a way that by taking into account the condition of the
⎫ ⎪ = 1 − h(η) ⎪ c X*+ ⎪ ⎪ c XM+(0,t ) K δr ⎞⎪ h(η) ⎛ ⎟ ⎜1 + = δ ⎠⎪ c X*+ (1 + K )γ ⎝ ⎪ ⎬ M (0,t ) cXL + δr ⎞ ⎪ Kh(η) ⎛ ⎜1 − ⎟ ⎪ = δ⎠ ⎪ c X*+ (1 + K )γ ⎝ ⎪ ⎪ c RW+2(d ,t ) D XW+1 ⎪ = η h ( ) c X*+ ⎪ DRW+2 ⎭ c XW+1(0,t )
(B.7)
with h(η) given by eq 33. If we assume that the complex formation and dissociation are at equilibrium even when current is flowing (total equilibrium conditions), the expression for the concentration profiles and the interfacial concentrations of the species involved in the facilitated ion transfer take the following form
⎞⎫ ⎟⎪ ⎟⎪ c X*+ ⎟⎪ ⎠⎪ ⎪ ⎞ ⎛ c XM+(x ,t ) h′(η) ⎜ ⎛⎜ x ⎞⎟⎟ ⎪ = erfc ⎪ γ(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠⎟⎠ ⎪ c X*+ ⎬ M (x ,t ) ⎞ ⎪ ⎛ cXL + Kh′(η) ⎜ ⎛⎜ x ⎞⎟⎟ ⎪ = erfc γ(1 + K ) ⎜⎝ ⎜⎝ 2 DMt ⎟⎠⎟⎠ ⎪ c X*+ ⎪ ⎪ ⎛ x − d ⎞⎪ c RW+2(x ,t ) D XW+1 ⎟⎟⎪ = h′(η) erfc⎜⎜ c X*+ DRW+2 ⎝ 2 D W2t ⎠⎪ ⎭
c XW+1(x ,t )
⎞ c R*+ ⎟ c X*+ ⎟⎟ ⎠
⎛ ⎜ x = 1 − h′(η) erfc⎜ ⎜ 2 D W+1t X ⎝
6463
(B.8)
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464
The Journal of Physical Chemistry A
Article
and
⎫ ⎪ ⎪ c X*+ ⎪ ⎪ M c X+(0,t ) h′(η) ⎪ = (1 + K )γ ⎪ c X*+ ⎪ ⎬ M cXL+(0,t ) ⎪ Kh′(η) = ⎪ (1 + K )γ ⎪ c X*+ ⎪ ⎪ c RW+2(d ,t ) D XW+1 ⎪ = ′ η ( ) h c X*+ ⎪ DRW+2 ⎭ c XW+1(0,t )
(17) Shirai, O.; Kihara, S.; Yoshida, Y.; Matsui, M. J. Electroanal. Chem. 1995, 389, 61−70. (18) Shirai, O.; Yoshida, Y.; Kihara, S. Anal. Bioanal. Chem. 2006, 386, 494−505. (19) Senda, M.; Kakiuchi, T.; Osakai, T. Electrochim. Acta 1991, 36, 253−262. (20) Molina, A.; Serna, C.; Ortuño, J. A.; Torralba, E. Int. J. Electrochem. Sci. 2008, 3, 1081−1091. (21) Langmaier, J.; Samec, Z. Electrochem. Commun. 2007, 9, 2633− 2638. (22) Molina, A.; Serna, C.; González, J.; Ortuño, J. A.; Torralba, E. Phys. Chem. Chem. Phys. 2009, 11, 1159−1166. (23) Molina, A.; Morales, I.; Lopez-Tenes, M. Electrochem. Commun. 2006, 8, 1062−1070. (24) Molina, A.; Morales, I. Electrochem. Commun. 2006, 8, 1453− 1460. (25) Molina, A.; Martínez-Ortíz, F.; Laborda, E.; Puy, J. Phys. Chem. Chem. Phys. 2010, 12, 5396−5404. (26) DeFord, D. D.; Hume, D. N. J. Am. Chem. Soc. 1951, 73, 5321− 5322. (27) Molina, A.; Ortuño, J. A; Serna, C.; Torralba, E.; González, J. Electroanalysis 2010, 22, 1634−1642. (28) Molina, A.; Serna, C.; Ortuño, J. A; González, J.; Torralba, E.; Gil, A. Anal. Chem. 2009, 81, 4220−4225. (29) Molina, A.; Ortuño, J. A.; Serna, C.; Torralba, E. Phys. Chem. Chem. Phys. 2010, 12, 13296−13303. (30) Cízek, J.; Koryta, J.; Koutecký, J. Collect. Czech. Chem. Commun. 1959, 663, 24.
= 1 − h′(η)
(B.9)
with h′(η) being given by eq 44.
■
AUTHOR INFORMATION
Corresponding Author
*Tel: +34 968 367524. Fax: +34 968 364148. E-mail:
[email protected].
■
ACKNOWLEDGMENTS We greatly appreciate the financial support provided by the ́ Dirección General de Investigación Cientifica y Técnica (Project Numbers CTQ2011-27049/BQU and CTQ200913023/BQU), and the Fundación SENECA (Project Number 08813/PI/08). Also, E.T. thanks the Ministerio de Ciencia e Innovación for the grant received. In addition, we thank referees 1 and 2 for their valuable comments.
■
REFERENCES
(1) Samec, Z. Pure Appl. Chem. 2004, 76, 2147−2180. (2) Reymond, F.; Fermín, D.; Lee, H. J.; Girault, H. H. Electrochim. Acta 2000, 45, 2647−2662. (3) Vanýsek, P.; Basáez Ramírez, L. J. J. Chil. Chem. Soc. 2008, 53, 1455−1463. (4) Bender, C. J. Chem. Soc. Rev. 1988, 17, 317−346. (5) Reymond, F.; Chopineaux-Courtois, V.; Steyaert, G.; Bouchard, G.; Carrupt, P.; Testa, B.; Girault, H. H. J. Electroanal. Chem. 1999, 462, 235−250. (6) Matsuda, H.; Kamada, Y.; Kanamori, K.; Kudo, Y.; Takeda, Y. Bull. Chem. Soc. Jpn. 1991, 64, 1497−1508. (7) Shao, Y.; Osborne, M. D.; Girault, H. H. J. Electroanal. Chem. 1991, 318, 101−109. (8) Gulaboski, R.; Ferreira, E. S.; Pereira, C. M.; Cordeiro, M. N.; Garau, A.; Lippolis, V.; Silva, A. J. Phys. Chem. C 2008, 112, 153−161. (9) Stephenson, M. J.; Holmes, S. M.; Dryfe, R. A. W. Electrochem. Commun. 2004, 6, 294−298. (10) Langmaier, J.; Samec, Z. Anal. Chem. 2009, 81, 6382−6389. (11) Kakutani, T.; Nishiwaki, Y.; Osakai, T.; Senda, M. Bull. Chem. Soc. Jpn. 1986, 59, 781−788. (12) Grell, E.; Funck, T.; Eggers, F. In Membranes: A Series of Advances; Eisenman, G., Ed.; Marcel Dekker: New York and Basel, 1975; Vol. 3, pp 1−113. (13) Senda, M.; Katano, H.; Yamada, M. J. Electroanal. Chem. 1999, 475, 90−98. (14) Zhan, D.; Xiao, Y.; Yuan, Y.; He, Y.; Wu, B.; Shao, Y. J. Electroanal. Chem. 2003, 553, 43−48. (15) Koryta, J. Electrochim. Acta 1978, 24, 293−300. (16) Langmaier, J.; Stejskalová, K.; Samec, Z. J. Electroanal. Chem. 2001, 496, 143−147. 6464
dx.doi.org/10.1021/jp2109362 | J. Phys. Chem. A 2012, 116, 6452−6464