This article is part of the series IFFM&AFM2011 .

Open Access Special Issue Article

Solution reaction design: electroaccepting and electrodonating powers of ions in solution

Keyan Li1, Min Li1 and Dongfeng Xue12*

Author affiliations

1 State Key Laboratory of Fine Chemicals, School of Chemical Engineering, Dalian University of Technology, Dalian 116024, People's Republic of China

2 State Key Laboratory of Rare Earth Resource Utilization, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, People's Republic of China

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:6  doi:10.1186/1556-276X-7-6


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/6


Received:6 September 2011
Accepted:5 January 2012
Published:5 January 2012

© 2012 Li et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

By considering a first-order variation in electroaccepting and electrodonating powers, ω±, induced by a change from gas to aqueous solution phase, the solvent effect on ω± for charged ions is examined. The expression of electroaccepting and electrodonating powers in the solution phase, ω±s, is obtained through establishing the quantitative relationship between the change of the ω± due to the solvation and the hydration free energy. It is shown that cations are poorer electron acceptors and anions are poorer electron donors in solution compared to those in gas phase. We have proven that the scaled aqueous electroaccepting power, ω+s, of cations can act as a good descriptor of the reduction reaction, which is expected to be applied in the design of solution reactions.

Keywords:
electroaccepting; electrodonating; solution phase; solvent effect

Introduction

With the rapid development of functional materials, novel micro/nanostructures of the materials are highly demanded to obtain advanced properties, which can be achieved by the rational design of solution-phase chemical reactions [1-6]. Therefore, it is of significance to thoroughly understand the reactivity of chemical species and the mechanism of chemical reactions to further realize the solution reaction design. Among many chemical reactivity indices, one quantity of importance is the electrophilicity, ω, introduced by Parr and co-workers [7]. They defined ω as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M1">View MathML</a>

(1)

where μ is the chemical potential and η is the chemical hardness of an N-electron system with total energy, E, defined as μ = (∂2 E/∂2N)v(r) and η = (∂2 E/∂2N)v(r). This index has been found to be helpful in analyzing the reactivity behaviors of a variety of compounds as well as the reaction mechanisms of diverse chemical processes [8,9]. As an important contribution to the ω, Gazauez et al. [10] argued that from a chemical perspective, it would make sense to differentiate the response of a system to the electron acceptance from the electron donation grounded on that the left and right derivatives of the total energy, EDFT(N), of an N-electron system with respect to the integer electron number, N, are different. By introducing an electron bath of nonzero chemical potential, μbath, with which the chemical species can exchange electrons, they proposed electroaccepting [ω+] and electrodonating [ω-] powers as the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M2">View MathML</a>

(2)

where the chemical potential, μ±, and the chemical hardness, η±, were defined as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M3">View MathML</a>

(3a)

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M4">View MathML</a>

(3b)

where I and A are the ionization potential and the electron affinity, respectively. A larger value of ω+ corresponds to a larger capability of accepting charges, whereas a smaller value of ω- implies a larger capability of donating charges.

Although some chemical phenomena have been rationalized by establishing the quantitative structure-reactivity relationships using these reactivity indices in the gas phase [11], the presence of solvent is bound to affect the reactivity behaviors of chemical substances. Therefore, studies on the reactivity indices such as ω and ω± in solution are quite necessary to reveal the accurate reactivity of chemical species in solution and further predict and design the solution phase reactions [12]. While several theoretical calculations about the solvent effect on the ω for various chemical species have been performed [12-14], the solvent effect on the ω± which are regarded as better descriptors of the donor-acceptor type interactions [10] has not received much attention to date. In this work, the solvent effect on the ω± is estimated by establishing a linear relationship between the change of the ω± due to the solvation and the hydration free energy, ΔGhyd. The values of aqueous electroaccepting power, ω+s, of 39 metal cations are quantitatively calculated, which are proven to be appropriate descriptors for the reduction reactions.

Method

In a previous study, Perez et al. [14] examined the solvent effect on the electrophilicity index, ω, by introducing a first-order finite variation in the ω due to the solvation

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M5">View MathML</a>

(4)

where Δμg→s and Δηg→s are the variations in μ and η from the gas to solution phase, respectively.

They rearranged the first contribution as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M6">View MathML</a>

(5)

where ΔEins is the insertion energy of the solute going into the solvent which is suggested as twice the solvation energy.

The second contribution in Equation 4 is rewritten as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M7">View MathML</a>

(6)

Finally, they deduced the expression of Δωg→s

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M8">View MathML</a>

(7)

where Δωg→s showed a linear dependence on the solvation energy, ΔGsolv, with a regression slope, γ. They used 18 well-known electrophilic ligands including hard electrophiles such as Li+ and Na+ to test this linear correlation and obtained good results (R = 0.9925, γ = 1.00765 at B3LYP/6-311G**and R = 0.9918, γ = 0.96843 at HF/6-311G**levels of theory).

Herein, we reconstruct the second contribution in Equation 4 which will directly lead to a quantitative expression for Δωg→s with a definite slope value, γ.

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M9">View MathML</a>

(8)

Substitution of Equations 5 and 8 into Equation 4 leads to the expression of Δωg→s.

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M10">View MathML</a>

(9)

Therefore, the global electrophilicity, ωs, in solution can be calculated by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M11">View MathML</a>

(10)

It should be noted that one key assumption in our approach is ΔNs, max ≈ ΔNg, max which could be justified by the data of Table 1 in Perez's work [14]. Our result, γ = 1, has turned out to be fairly consistent with Perez's regression value, i.e., γ = 1.00765 and γ = 0.96843, which thus approve the reasonableness of our approach to dealing with Δω2, g→s.

Table 1. Calculated electroaccepting power, ω+s, in aqueous solution and the absolute reduction potential, E°abs

Further, we try to extend our approach to examine the solvent effect on the ω±. For the charged ions, we suppose that the chemical potential, μ±bath, of the electron bath equals that of the parent atoms of ions since the charged ions become neutral atoms after accepting or donating the maximum amount of electrons. In addition, as the solvent only has little effect on the chemical potential, μ, of the neutral species [12,13,15], there exists a relationship as μ±bath = μs, atom μg, atom. The ion exchanges electrons from the bath to the point that its chemical potential, μ±, equals the value μ±bath with the maximum amount of electron flow:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M12">View MathML</a>

(11)

The first-order variation in the ω± leads to the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M13">View MathML</a>

(12)

The first part of Equation 12 in terms of the variation in μ± is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M14">View MathML</a>

(13)

where the energy change Δω±1,g→s due to the variation of the chemical potential from the gas to solution phase can be represented by ΔGhyd [16,17].

The second part of Equation 12 in terms of the variation in η± is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M15">View MathML</a>

(14)

Combining Equation 13 with Equation 14 yields

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M16">View MathML</a>

(15)

Therefore, the electroaccepting and electrodonating powers in solution, ω±s, can be calculated by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M17">View MathML</a>

(16)

Results and discussion

According to Pearson's viewpoint that cations are electron acceptors and anions are electron donors [16], we pay attention to the ω+s for cations and ω-s for anions. By using Equation 16, the ω+s values for 39 metal cations with charges from +1e to +3e are calculated and summarized in Table 1. From Table 1, we find that the solvation weakens the capacity of cations to accept electrons due to the negative values of Δω+, in agreement with the previous conclusions [12-16]. Unfortunately, it is impossible to quantitatively calculate the ω-s values for anions so far due to the absence of experimental electron affinities needed in Equation 3b. Herein, these values can be qualitatively estimated:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M18">View MathML</a>

(17)

Since a larger value of ω-s implies a smaller capability of donating charges, we can conclude from Equation 17 that the solvation also weakens the capacity of anions to donate electrons, which agrees with the general viewpoints [12-16].

Many liquid-phase chemical reactions involve the electron-transfer steps, and a key thermodynamic variable that describes the tendency of chemical species in solution to gain or lose electrons is the redox potential. The quantum-chemical computation approach to electrochemistry has become available very recently [18]. However, the estimation of redox potential by the quantum-chemical calculations is a great challenge due to the complexity of the processes involved in a typical electrochemical reaction [19]. For example, the complicated diffusion and adsorption processes on the electrode surface which should be necessarily taken into account in the quantum-chemical modeling of the reduction-oxidation reaction lead to the considerable system size and thus require strong computing power. Therefore, previous studies mainly focus on the one-electron reduction reactions between different oxidation states of transition metals to avoid modeling of an electrode-solution boundary [18-20]. In this work, we try to use the ω+s to describe the many-electron reduction reaction including both main- and sub-group metal cations. According to the reaction formula MZ+ (aq) + z/2H2 (g) = M (aq) + zH+ (aq), the absolute reduction potential, E°abs, can be calculated by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/6/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/6/mathml/M19">View MathML</a>

(18)

where E°abs(SHE) is the conventional reduction potential and E°SHE = 0 is the absolute standard hydrogen electrode potential. Note that the ω+s is the energy lowing associated with a maximum amount of electron flow between two species; it is reasonable to establish a correlation between zE°abs and ω+s. A good relationship shown in Figure 1 approves that our ω+s can act as an appropriate descriptor of the many-electron energy change. Moreover, this method is more simple and convenient compared to the quantum-chemical approach to the estimation of the E°abs.

thumbnailFigure 1. Plot of the ω+s versus zE°abs.

Except for the reduction reaction, the ω±s can also be expected to qualitatively and quantitatively predict other properties of ions in connection with ligand binding, hydrolysis processes, and stability of coordination compounds, etc. In addition, compilation of experimental data on solvation energies in nonaqueous solutions will make it possible to evaluate the corresponding electroaccepting and electrodonating powers, ω±, which will undoubtedly lead to the deeper understanding of the chemical reactivity of ions in these media.

Conclusions

By reconstructing a first-order variation of the ω due to the solvation, the linear relationship between the change in the ω and the solvation energy is reproduced, which suggests that our method is theoretically reasonable. The solvent effect on the electroaccepting and electrodonating powers, ω±, for charged ions is examined, and a definite quantitative expression for the aqueous ω±s is established. It is found that the solvation weakens the capability of both electron-accepting power of cations and electron-donating power of anions. A good relationship between the ω+s and E°abs shows the validity of the electroaccepting powers in determining the chemical reactivity of the ions in aqueous solution. It is expected that our ω±s will be helpful to achieve a better understanding of chemical properties of ions in solution and further be used in many aspects of solution chemistry such as the design of solution-phase reactions according to these indices.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KL participated in the design and coordination of the study and corrected the manuscript. ML assisted in the calculation of the data and prepared the manuscript initially. DX conceived the study, supervised, and corrected the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The financial support of the National Natural Science Foundation of China (Grant Nos. 50872016, 51125009 and 20973033) is greatly acknowledged.

References

  1. Liu J, Xue D: Cation-induced coiling of vanadium pentoxide nanobelts.

    Nanoscale Res Lett 2010, 5:1619. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Wu D, Jiang Y, Liu J, Yuan Y, Wu J, Jiang K, Xue D: Template route to chemically engineering cavities at nanoscale: a case study of Zn(OH)2 template.

    Nanoscale Res Lett 2010, 5:1779. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Wu J, Xue D: In situ precursor-template route to semi-ordered NaNbO3 nanobelt arrays.

    Nanoscale Res Lett 2011, 6:14. OpenURL

  4. Liu J, Xue D, K Li: Single-crystalline nanoporous Nb2O5 nanotubes.

    Nanoscale Res Lett 2011, 6:138. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  5. Yan C, Liu J, Liu F, Wu J, Gao K, Xue D: Tube formation in nanoscale materials.

    Nanoscale Res Lett 2008, 3:473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Liu J, Xue D: Hollow nanostructured anode materials for Li-ion batteries.

    Nanoscale Res Lett 2010, 5:1525. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Parr R, Szentpaly L, Liu S: Electrophilicity index.

    J Am Chem Soc 1999, 121:1922. Publisher Full Text OpenURL

  8. Roy D, Chattaraj P: Reactivity, selectivity, and aromaticity of Be32- and its complexes.

    J Phys Chem A 2008, 112:1612. PubMed Abstract | Publisher Full Text OpenURL

  9. Kostova I, Trendafilova N, Mihaylov T: Theoretical and spectroscopic studies of pyridyl substituted bis-coumarins and their new neodymium (III) complexes.

    Chem Phys 2005, 314:73. Publisher Full Text OpenURL

  10. Gazquez J, Cedillo A, Vela A: Electrodonating and electroaccepting powers.

    J Phys Chem A 2007, 111:1966. PubMed Abstract | Publisher Full Text OpenURL

  11. Chattaraj P, Sarkar U, Roy D: Electrophilicity index.

    Chem Rev 2006, 106:2065. PubMed Abstract | Publisher Full Text OpenURL

  12. Jaramillo P, Perez P, Fuentealba P, Canuto S, Coutinho K: Solvent effects on global reactivity properties for neutral and charged systems using the sequential Monte Carlo quantum mechanics model.

    J Phys Chem B 2009, 113:4314. PubMed Abstract | Publisher Full Text OpenURL

  13. Safi B, Balawender R, Geerlings P: Solvent effect on electronegativity, hardness, condensed Fukui functions, and softness, in a large series of diatomic and small polyatomic molecules: use of the EFP model.

    Phys Chem A 2001, 105:11102. Publisher Full Text OpenURL

  14. Perez P, Toro-Labbe A, Contreras R: Solvent effects on electrophilicity.

    J Am Chem Soc 2001, 123:5527. PubMed Abstract | Publisher Full Text OpenURL

  15. Valiev M, Yang J, Adams J, Taylor S, Weare J: Phosphorylation reaction in cAPK protein kinase-free energy quantum mechanical/molecular mechanics simulations.

    J Phys Chem B 2007, 111:13455. PubMed Abstract | Publisher Full Text OpenURL

  16. Pearson R: Ionization potentials and electron affinities in aqueous solution.

    J Am Chem Soc 1986, 108:6109. Publisher Full Text OpenURL

  17. Contreras R, Andres J, Safont V, Campodonico P, Santos J: A theoretical study on the relationship between nucleophilicity and ionization potentials in solution phase.

    J Phys Chem A 2003, 107:5588. Publisher Full Text OpenURL

  18. Moens J, Roos G, Jaque P, Proft F, Geerlings P: Can electrophilicity act as a measure of the redox potential of first-row transition metal ions?

    Chem Eur J 2007, 13:9331. PubMed Abstract | Publisher Full Text OpenURL

  19. Jaque P, Marenich A, Cramer C, Truhlar D: Computational electrochemistry: the aqueous Ru3+/Ru2+ reduction potential.

    J Phys Chem C 2007, 111:5783. Publisher Full Text OpenURL

  20. Uudsemaa M, Tamm T: Density-functional theory calculations of aqueous redox potentials of fourth-period transition metals.

    J Phys Chem A 2003, 107:9997. Publisher Full Text OpenURL

  21. Marcus Y: Thermodynamics of solvation of ions. Part 5.--Gibbs free energy of hydration at 298.15 K.

    J Chem Soc, Faraday Trans 1991, 87:2995. Publisher Full Text OpenURL